第一章 MATLAB应用.pptVIP

* MATLAB图形功能 常用命令： plot(x); plot(x,y,’s’) plot(x1,y1, ’s1’,x2, y2,’s2’) plot3 (x,y,z,’s’) plot3(x1,y1,z1, ’s1’,x2, y2,z2,’s2’) fplot(f,[a,b],N) ezplot(f,[a,b],N) ezplot(f,g,[a,b],N) t=0:0.1:2*pi; fplot(‘1/sin(x)’,[]) 新疆大学电气工程学院 图形修饰： grid on grid off xlabel(‘text’) ylabel(‘text’) title((‘text’) text(x,y,’string’) axis([xmin xmax ymin ymax]) hold on hold off subplot() t=0:0.01:10; y1=sin(t); y2=cos(t); y3=sin(t^.2); y4=sin(t^.0.5) figure(1) plot(t,y1,’r’,t,y2,’b’) grid on axis([0 10 -1.2 1.2]) xlabel(‘时间’) ylabel(‘幅值’) title(‘三角函数’) 新疆大学电气工程学院 MATLAB典型应用 1求解线性方程组 例：一曲线y=a+bx+cx2+dx3通过四点(1,3)、(2,4)、(3,3)、(4,-3)。试求系数a，b，c，d。 A=[1 1 1 1;1 2 4 8;1 3 9 27;1 4 16 64]; b=[3 4 3 -3]; d=det(A); if d~=0 x=A\b else disp(矩阵行列式为0) end 新疆大学电气工程学院 2. 函数极值 x=fmins(‘funname’,x0,options) x=fminsearch(‘funname’,x0,options) x=fminunc(‘funname’,x0,options) % 一元函数求极值 f=3*x-x^3; fplot(f,[-3 3]) x0=0; x=fminsearch(f,x0) ymin=3*x-x^3 % 多元函数求极值 x0=[1 0] x1=fminsearch(zhl1_fun,x0) ymin1=zhl1_fun(x1) function y=zhl1_fun(x) y=100*(x(2)-x(1)^2)^2+(1-x(1))^2; 新疆大学电气工程学院 3 代数方程求解 x=fzero(funname,x0) x=fsolve(funname,x0) 一元方程求解 多元元方程求解 f=2*sin(3*x)+0.5*x^2-1; fplot(f,[1.5 3]); grid on y=fzero(f,[1.5,3]) function y=zhl2_fun(x) y(1)=sin(x(1))+x(2)^2+log(x(2))-7; y(2)=3*x(1)+2*x(2)-x(3)^3+1; y(3)=x(1)+x(2)+x(3)-5; y=fsolve(zhl2_fun,[1 1 1]) 新疆大学电气工程学院 4 微分方程求解 x=ode45(funname,[t0 tf],x0,tol) 求 例如已知： 的解。 function y=zhl3_fun(t,x) y=[x(2);(1-x(1)^2)*x(2)-x(1)]; [t,x]=ode45(‘zhl3_fun’,[0 30],[1;0]); figure(1) plot(t,x(:,1),t,x(:,2)) xlabel(‘t’) ylabel(‘x(t)’) figure(2) plot(x(:,1),x(:,2)) 新疆大学电气工程学院 新疆大学电气工程学院 混沌简介： 是确定非线性系统的特有复杂运动形态 lorenz方程： 主要特点：初值敏感性、有界性、各态遍历性， 内随机性、分数维等 新疆大学电气工程学院 %lorenz方程求解 tf=150; x0=[0;0;0.0000001]; [t,x]=ode45(zhl_lorenz_fun,[0,tf],x0); figure(1);plot(t,x) figure(2);plot3(x(:,1),x(:,2),x(:,3));grid on figure(3) comet3(x(:,1),x(:,2),x(:,3)) function y=zhl_lorenz_fun(t,x) sigma=28;rho=10;beta=8/3; y=[-beta*x(1)+x(2)*x(3);-rho