MATLAB基础与控制系统CAD仿真练习1(9)

控制系统计算机辅助设计,matlab仿真的一些小练习

于是，我们可知：

（1） k值为1.8447e+019；

（2） y=x*sin(x) 在0<x<100条件下共有16个峰值，分别为：1.8197，7.9167，

14.1724，20.4448，26.7222，33.0019，39.2826，45.5641，51.8459，58.1281，64.4104，70.6929，76.9755，83.2582，89.5410，95.8238。

3-3．绘制下列函数的图形。

（1）sin(1/t), -1<t<1 (2) 1 cos3(7t) -1<t<1

解答：

编写M-file并命名为zuoye23.m

t=-1:0.01:1; y=sin(1./t);

subplot(1,2,1),plot(t,y),grid on;

xlabel('t'),ylabel('sin(1/t)'); %标横纵坐标

z=1-(cos(7.*t)).^3;

subplot(1,2,2),plot(t,z),grid on;

xlabel('t'),ylabel('1-cos^3(7t)'); %标横纵坐标

随后在Command Window中运行zuoye23.m，回车得到结果（如图1）：

图1

控制系统计算机辅助设计,matlab仿真的一些小练习

其中，左边一个图表为（1）sin(1/t), -1<t<1的函数图形，右边一个图表为(2)

1 cos3(7t) -1<t<1的函数图形。

3-4．已知元件的实验数据如下，拟合这一数据，给出其特性方程并绘出图形。

解答：

从表中可知X是以1为间距的等差数列。 编写M-file并命名为zuoye24.m

x=0.01:9.01; y=[2.5437 7.8884 9.6242 11.6071 11.9727 13.2189 14.2679 14.6134 15.4045 15.0805]; plot(x,y,'b'),grid on; hold on;

p=polyfit(x,y,5) %特性方程的各阶系数 y1=polyval(p,x); plot(x,y1,'r');

随后在Command Window中运行zuoye24.m，回车得到结果：

p =

0.0015 -0.0450 0.4962 -2.5636 7.2232 2.5195

以及拟合图形（如图2）：