控制系统数字仿真和CAD第三章习题.docVIP

3-1． 求解下列线性方程，并进行解得验证： (1),(2) 由A*X=B得：X=A\B 解： a=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] a = 7 2 1 -2 9 15 3 -2 -2 -2 11 5 1 3 2 13 b=[4 7 -1 0] b = 4 7 -1 0 x=a\b x = 0.4979 0.1445 0.0629 -0.0813 (2)解： a=[5 7 6 5 1 7 10 8 7 2 6 8 10 9 3 5 7 9 10 4 1 2 3 4 5] a = 5 7 6 5 1 7 10 8 7 2 6 8 10 9 3 5 7 9 10 4 1 2 3 4 5 b=[24 96 34 136 36 144 35 140 15 60] b = 24 96 34 136 36 144 35 140 15 60 x=a\b x = 1.0000 4.0000 1.0000 4.0000 1.0000 4.0000 1.0000 4.0000 1.0000 4.0000 3-2．进行下列计算，给出不使用for和while等循环语句的计算方法。 （1） 解:根据等比数列求和方法，在利用matlab中的m文件，编写程序求解。 M文件为 n=64; q=2; k=(1-q^n)/(1-q); disp(k的值为); disp(k); 保存文件q1.m 在matlab命令框中输入 q1 k的值为 1.8447e+019 (2)求出y=x*sin(x) 在0x100条件下的每个峰值 解：画出图形 x=0:0.01:100; y=x.*sin(x); plot(x,y); grid on title(y=x*sin(x)) xlabel(x) ylabel(y) 方法1。从图形中不难看出峰值点取决于函数sin(x),即在sin(x)为峰值时，y就得到峰值。所以求取函数的峰值转化为求取正弦函数波峰问题。而sin(x)在x= +2k(k为整数)，所以求取y在上述x时刻的数值就是峰值。 在matlab命令行里键入 x=pi/2:pi*2:100; y=x.*sin(x) %注意是。*不是*% 得到结果y=1.5708 7.8598 14.1481 20.4350 26.7198 33.0019 39.2804 45.5549 51.8245 58.0887 64.3467 70.5978 76.8414 83.0769 89.303 95.5204 方法2. a=size(y) a=1 1001 b=([y(2:1000)][y(1:999)])([y(2:1000)][y(3:1001)]); at=find(b==1); disp(y(at)) 就可以找到最大值点 3-3．绘制下面的图形。 （1）sin(1/t),-1t1 (2) -1t1 (1)解: t=-1:0.01:1; y=sin(1./t); %注意是./不是/% Warning: Divide by zero. plot(t,y) grid on xlabel(t) ylabel(y) title(y=sin(1/t)) （2）解： t=-1:0.01:1; y=1-(cos(7.*t)).^3; %注意是.*与.^% plot(t,y) grid on xlabel(t) ylabel(y) title( y=1-cos(7t)^3) 3-4．已知元件的实验数据如下，拟合这一数据，并尝试给出其特性方程。 X 0.0100 1.0100 2.0100 3.0100 4.0100 Y 2.5437 7.8884 9.6242 11.6071 11.9727 X 5.0100