求一个MATLAB自定义函数并调用的实例程序

Question

如题所述

Answer 1

没其他条件吗？那例子太多了。

给你一个刚才帮别人调的Jacobi迭代求解线性方程组的程序吧。

 

函数（保存成Jacobi.m）：

function x=Jacobi(A,b,x0)
epsilon=0.001;
n=length(b);
x1=zeros(n,1);
M=100;
for k=1:1:M
    for i=1:1:n
        s=sum(A(i,:)*x0(:))-A(i,i)*x0(i);
        x1(i)=(b(i)-s)/A(i,i);
    end
    fprintf('%.4f ',x1);
    fprintf('\n ');% 按格式输出数据
    if norm(x1-x0,'inf')<epsilon
        break
    end
    x0=x1;
end
x=x1;

 

调用：

A = [10., -1., 2., 0.;
    -1., 11., -1., 3.;
    2., -1., 10., -1.;
    0.0, 3., -1., 8.];
b = [6., 25., -11., 15.]';
x0 = zeros(4,1);
x=Jacobi(A,b,x0)