计算机系统建模与仿真论文,计算机建模与仿真-实验指导书-2013.doc

计算机建模与仿真-实验指导书-2013.doc

《计算机建模与仿真》

实验指导书

金陵科技学院机电工程学院电气系

2013年3月

实验一、MATLAB的基本语法

一、实验目的和要求

1熟悉数组和矩阵的运算；

2熟悉二维图的绘制。

二、实验仪器和设备

计算机一台

三、实验过程

记录下面习题的程序和运行后的结果。

计算与的数组乘积。

>> a=[6 9 3;2 7 5];b=[2 4 1;4 6 8];

>> a.*b

ans =

12 36 3

8 42 40

2、对于，如果，，求解X。

>> A=[4 9 2;7 6 4;3 5 7];B=[37;26;28];

>> A\B

ans =

-0.5118

4.0427

1.3318

3、已知：，分别计算a的数组平方和矩阵平方，并观察其结果。

>> a=[1 2 3;4 5 6;7 8 9];

>> a.^2

ans =

1 4 9

16 25 36

49 64 81

>> a^2

ans =

30 36 42

66 81 96

102 126 150

角度，求x的正弦、余弦、正切和余切。

>> x=[30 45 60];

>> y=x/180*pi

y =

0.5236 0.7854 1.0472

>> sin(y)

ans =

0.5000 0.7071 0.8660

>> cos(y)

ans =

0.8660 0.7071 0.5000

>> tan(y)

ans =

0.5774 1.0000 1.7321

>> cot(y)

ans =

1.7321 1.0000 0.5774

5、将矩阵、和组合成两个新矩阵：

(1)组合成一个4(3的矩阵，第一列为按列顺序排列的a矩阵元素，第二列为按列顺序排列的b矩阵元素，第三列为按列顺序排列的c矩阵元素，即

(1)>> a=[4 2;5 7];b=[7 1;8 3];c=[5 9;6 2];

>> d=[a(:) b(:) c(:)]

d =

4 7 5

5 8 6

2 1 9

7 3 2

(2)按照a、b、c的列顺序组合成一个行矢量，即

(2)>> e=[a(:);b(:);c(:)]'

e =

4 5 2 7 7 8 1 3 5 6 9 2

解方程组。

>> a=[2 9 0;3 4 11;2 2 6];b=[13;6;6];

>> a\b

ans =

7.4000

-0.2000

-1.4000

7、绘制曲线，x的取值范围为[-5,5]。

>> x=-5:0.1:5;

>> y=x.^3+x+1;

>> plot(x,y)

8， 有一组测量数据满足，t的变化范围为0~10，用不同的线型和标记点画出a=0.1、a=0.2和a=0.5三种情况下的曲线,在图中添加标题，并标识出各曲线a的取值，并添加图例框。

>> t=0:0.1:10;

>> y1=exp(-0.1.*t+1);y2=exp(-0.2.*t+1);y3=exp(-0.5.*t+1);

>> plot(t,y1,'m:',t,y2,'b+',t,y3,'r-.')

>> title('y=e^{-at+1}')

>> gtext('a=0.1')

>> gtext('a=0.2')

>> gtext('a=0.3')

>>legend('y1=e^{-0.1t+1}','y2=e^{-0.2t+1}','y2=e^{-0.5t+1}')

9，表中列出了4个观测点的6次测量数据，将数据绘制成为条形图。

第1次第2次第3次第4次第5次第6次观测点1367428观测点2673247观测点3972584观测点4643274>> y=[3 6 7 4 2 8;6 7 3 2 4 7;9 7 2 5 8 4;6 4 3 2 7 4];

>> bar(y)

10， x= [66 49 71 56 38]，绘制饼图，并将第五个切块分离出来。

>> x=[66 49 71 56 38];

>> subplot(1,2,1),pie(x)

>> explode=[0 0 0 0 1];

>> subplot(1,2,2)