MATLAB数值运算实验报告Word文件下载.docx
《MATLAB数值运算实验报告Word文件下载.docx》由会员分享，可在线阅读，更多相关《MATLAB数值运算实验报告Word文件下载.docx（10页珍藏版）》请在冰豆网上搜索。
（3）多项式的乘、除法分别用函数conv和deconv实现。
（4）多项式求根用函数roots
（5）多项式求值用函数polyval
练习1：
求（s^2+1）（s+3）（s+1）/（s^3+2*s+1）的“商〞及“余〞多项式
2.多项式插值和拟合
有一组实验数据如表所示
X
1
2
3
4
5
6
7
8
9
10
Y
16
32
70
142
260
436
682
1010
1432
1960
请分别用拟合〔二阶至三阶〕和插值〔线性和三次样条〕的方法来估测X=9.5时Y的值。
3.数值微积分
（1）差分使用diff函数的实现
（2）可以用因变量和自变量差分的结果相除得到数值微分
（3）Cumsum函数求累计积分，trapz函数用梯形法求定积分，即曲线的面积
练习：
如图瑞士地图，为了算出其国土面积，首先对地图作如下测量：
以由西向向为X轴，由南到北方为Y轴，选择方便的原点，并将从最西边点到最东边界点在X轴的区间适当划分假设干级，在每个分点的Y方向测出南边界点和北边界点的Y坐标Y1和Y2,这样就得到了下表，根据地图比例知道18mm相当于40km，试有测量数据计算瑞士国土近似面积，与其准确值41228km^2比拟。
10.5
13
17.5
34
40.5
44.5
48
56
61
68.5
76.5
80.5
91
Y1
44
45
47
50
38
30
36
41
46
Y2
59
72
93
100
110
117
118
116
96
101
104
106.5
111.5
123.5
136.5
146
150
157
158
43
37
33
28
65
55
54
52
66
68
121
124
122
83
81
82
86
85
4.构造数组与细胞数组
（1）构造数组的创立
（2）构造数组的操作
创立一构造数组stusorce，其域为：
No,Name,English，Math，Chinese，Total，Average。
构造数组的大小为2*2。
（3）细胞数组的创立
（4）细胞数组的操作
创立一大小为2*2细胞数组stucell，其元素的类型分别为：
构造类型、字符串、矩阵和细胞类型。
四.实验具体过程及数据分析
1.
S1=[242]
S1=
242
roots（S1）
ans=
-1
S=[13-509]
S=
13-509
S1=[2311]
2311
S2=[13-547]
S2=
13-547
S3=conv（S1,S2）
S3=
291026-296577
S4=deconv（S3,S1）
S4=
polyval（S1,3）
32
x=1:
x=
12345678910
y=polyval（S1,x）
y=
81832507298128162200242
clearall
s1=[101]
s2=[13]
s3=[11]
s4=conv（s1,s2）
Y=conv（s4,s3）
X=[1021]
[Q,R]=deconv（Y,X）
s1=
101
s2=
13
s3=
11
s4=
1313
Y=
14443
X=
1021
Q=
14
R=
002-5-1
poly2sym（Q）
x+4
poly2sym（R）
2*x^2-5*x-1
conv（Q,X）+R-Y
00000
101
poly2sym（Q）
poly2sym（R）
ans=
conv（Q,X）+R-Y
2.
x=1:
y=[163270142260436682101014321960]
p1=polyfit（x,y,1）
y1=polyval（p1,9.5）
3.
2:
diff（x）
x=linspace（0,2*pi,100）;
y=sin（x）;
plot（x,y）
y1=diff（y）./diff（x）;
plot（x（1:
end-1）,y1）
x=ones（1,10）
cumsum（x）
x=linspace（0,pi,100）;
s=trapz（x,y）
练习2：
x=[710.51317.53440.544.548566168.576.580.59196101104106.5111.5118123.5136.5142146150157158];
y1=[444547505038303034363441454643373328326555545250666668];
y2=[4459707293100110110110117118116118118121124121121121116122838182868568];
X=x./18*40;
Y1=y1./18*40;
Y2=y2./18*40;
t1=trapz（X,Y1）,t2=trapz（X,Y2）,t=t2-t1
4.
〔1〕>
student.number='
20050731001'
;
student.name='
jack'
student
（2）.number='
20050731002'
student
（2）.name='
lucy'
；
〔2〕
student
（1）.subject=[]
student
（1）.sorce=[]
student
fieldnames（student）
getfield（student,{2},'
name'
）
student=rmfield（student,'
subject'
student=setfield（student,{1},'
sorce'
90）;
stuent
（2）.sorce=88;
〔3〕
A={'
How
are
you!
'
ones（3）;
[1
2;
3
4],{'
cell'
}};
B（1,1）={'
Hello
world'
};
B（2,1）={[
1
2
4]};
〔4〕
ans1=A（1,1）
ans2=A{1,1}
whos
ans1
ans2
celldisp（A）
a1=A{2,1}（1,2）
[a2
a3]=deal（A{1:
2}）
五.实验原始记录
六.实验心得、体会及思考
通过本次实验，我更熟练的掌握了MATLAB的数值运算及其运算中所用到的函数，更好的明白了一些MATLAB数值运算的根本操作。
做实验不仅仅要我们有清晰的思路，还要有一丝不苟的态度，认真严谨的实验才能得出更准确的实验结果。
升级会员 