Question

如何用matlab解这个方程组

40=x^(7/4)*1.81371+y^(7/4)*2.39720,
60=x^(7/4)*1.81371+z^(7/4)*4.79441,
0.007854*x=0.0050265*y+0.0050265*z

Answer 1

你好，这个方程组使用赛德尔迭代法，matalb没有自带这个函数，所以自己写个。
function [x,n,data]=Seidelgroudfun(f,x0,tol,max1)
N=length(x0);
for k=1:max1
x=x0;
for j=1:N
a=feval(f,x);
x(j)=a(j);
end
err=abs(norm(x-x0));
x0=x;
n=k;
data(:,n)=x;
if(err<tol)
break
end
end
然后，对你的方程组进行编写，也写为函数文件。
function out=myfun(X);
x=X(1);
y=X(2);
z=X(3);
out=zeros(3,1);
out(1)=(0.0050265*y+0.0050265*z)/0.007854;
out(2)=((40-x^(7/4)*1.81371)/2.39720)^(4/7);
out(3)=((60-x^(7/4)*1.81371)/4.79441)^(4/7);
最后，可以直接在命令窗口调用了
clc
clear all
x0=[1 1 1];
tol=1e-6;
max1=20;
[x,n,data]=Seidelgroudfun('myfun',x0,tol,max1);
祝好。

Answer 2

不是道