matlab经典例题及答案,matLAB经典例题及答案.pdf

.

一 对以下数据分别作二次,三次多项式拟合,并画出图形.

x 1:16;

y [4, 6.4, 8, 8.4, 9.28, 9.5, 9.7, 9.86, 10, 10.2, 10.32, 10.42, 10.5,

10.55, 10.58, 10.6];

答：程序如下

(1)x (1:16);

y erf(x);

p polyfit (x,y,2);

f polyval (p,x);

plot (x,y,x,f);

结果 p

-0.0010 0.0202 0.9096

图形

1.04

1.02

1

0.98

0.96

0.94

0.92

0.9

0.88

0.86

0.84

0 2 4 6 8 10 12 14 16

(2)y [4, 6.4, 8, 8.4, 9.28, 9.5, 9.7, 9.86, 10, 10.2, 10.32, 10.42, 10.5,

10.55, 10.58, 10.6];

y erf(x);

p polyfit (x,y,3)

f polyval (p,x);

plot (x,y,x,f)

结果

P

0.0002 -0.0071 0.0628 0.8404

图形

1.04

1.02

1

0.98

0.96

0.94

0.92

0.9

0.88

0.86

0.84

0 2 4 6 8 10 12 14 16

.

二 在[0,4pi]画sin(x),cos(x)(在同一个图象中); 其中cos(x)图象用红色

小圆圈画.并在函数图上标注 “y sin(x)”, “y cos(x)” ,x轴,y轴,标题

为 “正弦余弦函数图象”.

答：程序如下

x [0:720]*pi/180;

plot (x,sin (x),x,cos(x),'ro');

x [2.5;7];

y [0;0];

s ['y sin (x)';'y cos(x)'];

text (x,y,s);

xlabel ('正弦余弦函数图象'),ylabel ('正弦余弦函数图象')

图形如下

1

0.8

0.6

0.4

象

图 0.2

数

函 0 y sin(x ) y cos(x )

弦

余

-0.2

弦

正

-0.4

-0.6

-0.8

-1

0 2 4 6 8 10 12 14

正 弦 余 弦 函 数 图 象

三.选择一个单自由度线性振动系统模型，自定质量、弹簧刚度、阻尼、激振力

等一组参数，分别编程 (m文件)计算自由和强迫振动时的响应，并画出振动曲

线图。

(要求画出该单自由度线性振动系统模型图)

其中质量为m 1000kg,弹性刚度k 48020N/m,阻尼c 1960N.s/m,激振力f(t) 0.

阻尼比ζ的程序 p 1960/(2*sqrt(48020*1000))

求得 p 0.1414 而p为阻尼比ζ

强迫振动时的响应程序

g tf([-1 0 1],[48020 0 48020*1.98 48020]);

bode (g)

图形