湖 北 汽 车 工 业 学 院 《MATLAB及应用》实验指导书

第一题

Clear,clc

z=@(x)x(1)^2-(x(2)-1)^2;

[x,fvalue,flag,output]=fminsearch(z,[0,0]) y=x(2);x=x(1);z=x^2-(y-1)^2; x,y,z

第二题

Clear,clc

fun=@(x,y)sin(2.*y).*(x+y)./(x.^3+5); q=dblquad(fun,1,20,1,10)

第三题

function dydt=vdp1(t,y)

dydt=[y(2);100*(1-y(1)^2)*y(2)-y(1)-5];

新建脚本

[t,y]=ode45(@vdp1,[0 16],[5;10]); 第四题 Clc,clear a=10000;t=0 while a<1000000 a=a*1.06+10000 t=t+1; end t

第五题

Clc,clear

24

湖 北 汽 车 工 业 学 院 《MATLAB及应用》实验指导书

[X,Y]=meshgrid([-16:0.4:16],[-4:0.4:4]) Z=X.^2/16-Y.^2/4

subplot(2,2,1),mesh(X,Y,Z) subplot(2,2,2),surf(X,Y,Z)

第六题

25