资源描述：
《系统建模与仿真-第9次课-第四章ppt课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、[例4.3]试用循环语句和条件语句编程实现下述方波信号，要求：仿真步长取为h=0.1s；建立两个向量t和U；用MATLAB画出此方波信号。1程序h=0.1;fork=1:1:101t1=h*(k-1);t(k)=t1;ift1>=0&t1<2u1=5;elseift1>=2&t1<4u1=0;elseift1>=4&t1<6u1=5;elseift1>=6&t1<8u1=0;elseu1=5;endu(k)=u1;endplot(t,u,'r-');axis([01105.5]);xlabel('t/s');ylabel('U/V');24.2.6MATLAB的M文件
2、及M函数的编写与调用MATLAB提供了两种文件格式，其中一种是普通的ASCII码构成的文件，在这样的文件中只有MATLAB语言所支持的语句，它类似于DOS下的批处理文件，这种文件称为M文件，它的执行方式很简单，用户只需在MATLAB的提示符>>下键入该M文件的文件名，这样MATLAB就会自动执行该M文件中的各条语句。MATLAB的另一种，也是最常用的特殊M文件称为MATLAB函数，这样的函数是由function语句引导的，其基本格式如下：function返回变量列表=函数名（输入变量列表）注释说明语句段函数体语句3如果返回变量多于一个，则用方括号给括起来。[例4.3]
3、用户想要建立函数，则可编写M函数如下functiony=exam(x)y=cos(x)-exp(3*x);其中exam既是该函数名，也是该文件名，即该文件名为exam.m。[例4.4]用户想要建立一个函数形式的一阶微分方程组4则可编写M函数如下functiondx=exam1(t,x)dx(1)=x(2)-4*x(3);dx(2)=cos(x(1));dx(3)=-2*x(2);dx=dx';当MATLAB想要调用函数时，只需在程序中引用该函数对应的文件名即可。5functiondx=exam1(t,x)dx=[x(2)-4*x(3);cos(x(1));-2*x(2
4、)];上述函数还可以写成6[t,X]=ode45(‘exam1’,[0,10],[-2,0.2,0.4]);X代表三个列向量组成的长方矩阵，这三个列向量分别为：现在欲求解该微分方程组，则只要给定其初值条件及仿真的时间区间，即可调用函数。而在函数ode45()中再调用函数exam1。ode45()7如果模型是线性的，则M函数还可进一步简写。举例线性系统functiondx=f(t,x)A=[1,-2,0;-3,1,-1;-4,0,1];dx=A*x;84.3MATLAB的绘图功能MATLAB允许用户同时打开若干个图形窗口，打开一个图形窗口可以由命令figure()函数来
5、完成，每产生一个图形窗口，伴随着都有一个图形窗口的句柄。当然用户还可以用close()函数来关闭该图形窗口。4.3.1MATLAB下二维图形的绘制简单的说，绘制二维图形的函数为plot()函数，其基本调用格式为[例4.5]用户想要绘制一个周期内的正弦曲线，则首先应当用t=0:0.1:2*pi命令来产生自变量t，然后键入命令y=sin(t)，再调用plot(t,y)。绘制的曲线如下plot(x,y)901234567-1-0.8-0.6-0.4-0.200.20.40.60.81图4.1用MATLAB绘制的正弦曲线10MATLAB还允许在一个窗口同时绘制多条曲线，例如下
6、面的命令>>t=0:0.1:2*pi;y=[sin(t);cos(t)];plot(t,y)可以同时产生一组曲线如下01234567-1-0.8-0.6-0.4-0.200.20.40.60.81图4.2用MATLAB绘制的正弦曲线及余弦曲线11为了区分不同曲线的不同颜色，还可以在plot函数中加入选项，参见下例。t=0:0.1:2*pi;y1=sin(t);y2=cos(t);y3=sin(t).*cos(t);plot(t,y1,'-',t,y2,':',t,y3,'x');[例4.6]在同一张图中绘制曲线，，。则可以先编写一个M文件如下或者使用命令，该文件还可写
7、为：holdon12设这个文件的名为exa.m，则在MATLAB的工作空间中键入文件名exa，回车后就可得到曲线如下t=0:0.1:2*pi;y1=sin(t);y2=cos(t);y3=sin(t).*cos(t);plot(t,y1,'-’);holdon;plot(t,y2,':’);plot(t,y3,'x');.*含义？1301234567-1-0.8-0.6-0.4-0.200.20.40.60.81图4.3MATLAB用不同线型绘制的三条曲线14在上面的文件中，plot函数中就加入了附加选项：'-'、':'及'x'，他们分别表示实线、虚
