x=0:0.1:20;
y1=1-1/sqrt(1-0.3^2)*exp(-0.3*x).*sin(sqrt(1-0.3^2)*x+acos(0.3))
plot(x,y1,'r') %画阻尼系数为 0.3 的曲线
hold on
y2=1-1/sqrt(1-0.707^2)*exp(-0.707*x).*sin(sqrt(1-0.707^2)*x+acos(0.707))
plot(x,y2,'g') %画阻尼系数为 0.707 的曲线
y3=1-exp(-x).*(1+x)
plot(x,y3,'b') %画阻尼系数为 1 的曲线
(3) 选择 M 文件编辑器菜单“Debug”——“Run”,就可以在图形窗中看到如图 5.2
所示的曲线。
查看工作空间的变量:
whos
Name Size Bytes Class
x 1x201 1608 double array
y1 1x201 1608 double array
y2 1x201 1608 double array
y3 1x201 1608 double array
Grand total is 804 elements using 6432 bytes
5.1.4 M 函数文件
函数文件的特点:
(1) 第一行总是以“function”引导的函数声明行;
函数声明行的格式:
function [输出变量列表] = 函数名(输入变量列表)
(2) 函数文件在运行过程中产生的变量都存放在函数本身的工作空间;
评论0
最新资源