2025年MATLAB考试必考题 .pdfVIP

穷则独善其身，达则兼善天下。——《孟子》

一．线性方程组求解

(1)求方程的解（三种方法）。

(2)将方程右边向量元素b为[0.5:0.01:0.6]，再求解，并比较b的变化和解的相对变化。

33

解（1）○A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6];

1

b=[0.95,0.67,0.52];

x1=A\b;

x2=inv(A)*b;

x3=pinv(A)*b

○A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6];

2

b=[0.95,0.67,0.52];

[L,U]=lu(A);

x1=U\(L\b);

x2=inv(U)*(inv(L)*b);

x3=pinv(U)*(pinv(L)*b)

（2）A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6];

b1=0.95;

b2=0.67;

b3=[0.5:0.01:0.6];

fori=1:length(b3);

b(1)=b1;

b(2)=b2;

b(3)=b3(i);

x0=A\b;

x1(i)=x0(1);

x2(i)=x0(2);

x3(i)=x0(3);

end

subplot(2,2,1)

plot(b3,x1)

holdon

plot(b3,x2,r-);

holdon

plot(b3,x3,g--)

holdoff

二．非线性方程的数值解

去留无意，闲看庭前花开花落；宠辱不惊，漫随天外云卷云舒。——《幽窗小记》

2

sinxylnz70

y3

3x2z10,初值x1,y1,z1。

000

xyz50

解：functionf=fxx(x)

f=[sin(x(1))+x(2).^2+log(x(3))-7;3*x(1)+2^x(2)-x(3).^3+1;x(1)+x(2)+x(3)-5];

end

x=fsolve(@fxx,[1,1,1])

三．Montecarlosimilaion

例10：MonteCarlosimulation

求曲线y=x^2与曲线x=y^2所围成的面积。

解：%%MentoCarlosimulation

x3=rand(1,100000);

y3=rand(1,100000);

Y3=sqrt(x3);