matlab结课论文doc.docx
- 文档编号:4865681
- 上传时间:2022-12-11
- 格式:DOCX
- 页数:9
- 大小:151.57KB
matlab结课论文doc.docx
《matlab结课论文doc.docx》由会员分享,可在线阅读,更多相关《matlab结课论文doc.docx(9页珍藏版)》请在冰豆网上搜索。
matlab结课论文doc
北京林业大学全校公共选修课结课作业
MATLAB程序应用
姓名:
课程序号:
51
班级:
电气10-1
学号:
101054118
2011年12月
1、(指选题1)有一天小猴摘下了若干个桃子,当即吃掉了一半,还觉得不过隐,又多吃了2个。
第二天接着吃下剩下的一半,又多吃了2个。
以后每天都是吃掉尚存的一半零两个。
到第十天早上,小猴准备吃桃子时,看到只剩下n个桃子了。
问小猴第一天共摘下多少个桃子?
1)写出求解小猴第一天共摘下多少个桃子的数学表达式
2)根据确定的数学表达式,利用matlab编写出一个函数文件用于求解
3)编写一个命令文件,完成输入一个n值,然后调用上述函数求解,并将结果输出显示出来
【思路和方法】
假设猴子第一天吃掉……
因此本题可以转化利用matlab求解级数的问题……(正文格式:
中文宋体,西文timesNewRoman,1.25倍行距)
【源程序】
inputn%输入n
【运行结果】
>>display('Pleaseinputthevalueofn');
n=input('n=');
fori=1:
9;
n=(n+2)*2;
end
m=n
Pleaseinputthevalueofn
n=
2、(指选题2)
1、分别用矩阵求逆、矩阵除法、矩阵分解求下面线性方程组的解;
矩阵除法:
A=
1-23-4
01-11
1301
0-731
>>b=[4,-3,1,-3]
b=
4-31-3
>>x=A/b
x=
0.7143
-0.2000
-0.2286
0.6000
矩阵求逆
A=[1,-2,3,-4;0,1,-1,1;1,3,0,1;0,-7,3,1]
A=
1-23-4
01-11
1301
0-731
>>b=[4,-3,1,-3]
b=
4-31-3
x=inv(A)*b
x=
0.7143
-0.2000
-0.2286
0.6000
矩阵分解法(LU分解):
[L,U]=lu(A);
>>x=U/(L/b)
x=
0.7143
-0.2000
-0.2286
0.6000
4、某乡镇企业某产品年销售额
与所获利润
从2000年到2010年的数据如下:
年度
2000
2001
2002
2003
2004
2005
2006
2007
2008
2009
2010
销售额x
6.1
7.5
9.4
10.7
14.6
17.4
21.1
24.4
29.8
32.9
34.3
纯利润
4.5
6.4
8.3
8.4
9.7
11.5
13.7
15.4
17.7
20.5
22.3
试求
对
的经验回归直线方程,并做回归分析,对回归的参数进行区间估计;
x=[6.1,7.5,9.4,10.7,14.6,17.4,21.1,24.4,29.8,32.9,34.3];
>>y=[4.5,6.4,8.3,8.4,9.7,11.5,13.7,15.4,17.7,20.5,22.3];
>>p=polyfit(x,y,1)
p=
0.56391.9096
xi=6.1:
34.3;
>>yi=polyval(p,xi)
yi=
Columns1through8
5.34915.91306.47687.04077.60458.16848.73229.2961
Columns9through16
9.860010.423810.987711.551512.115412.679213.243113.8069
Columns17through24
14.370814.934615.498516.062316.626217.190017.753918.3177
Columns25through29
18.881619.445420.009320.573121.1370
plot(x,y,':
o',xi,yi,'-*')
3、(自选题)
>>surf(x,y,z)
title('surf(x,y,z)')
t=0:
pi/20:
2*pi;
[x,y,z]=cylinder(2+sin(t),30);
subplot(2,1,2);
surf(x,y,z);
4、(自选题)
我们是电气工程及其自动化专业,对电路分析有时会需要进行傅里叶变换和拉普拉斯变换,通过学习这门课,我们发现matlab可以实现对一个函数进行快速傅氏变换(fft),以下为具体操作:
对函数x(t)=12sin(10πt+π/2)+10(cos20πt)进行fft:
N=128
N=
128
>>T=1
T=
1
>>t=linspace(0,T,N);
>>x=12*sin(2*pi*5*t+pi/2)+10*cos(2*pi*10*t)
x=
Columns1through8
22.000020.435616.05249.70932.6203-3.9179-8.8128-11.4085
Columns9through16
-11.6121-9.8832-7.0931-4.2877-2.4143-2.0769-3.3828-5.9146
Columns17through24
-8.8369-11.1107-11.7627-10.1454-6.1215-0.12976.889213.6876
Columns25through32
18.983221.746121.430718.099512.40905.4608-1.4463-7.1059
Columns33through40
-10.6652-11.7996-10.7577-8.2703-5.3457-3.0031-2.0086-2.6789
Columns41through48
-4.8010-7.6861-10.3460-11.7479-11.0862-8.0045-2.71194.0343
Columns49through56
11.079417.120020.992721.936419.764014.90528.30901.2290
Columns57through64
-5.0568-9.5273-11.6334-11.3977-9.3774-6.4998-3.8129-2.2127
Columns65through72
-2.2127-3.8129-6.4998-9.3774-11.3977-11.6334-9.5273-5.0568
Columns73through80
1.22908.309014.905219.764021.936420.992717.120011.0794
Columns81through88
4.0343-2.7119-8.0045-11.0862-11.7479-10.3460-7.6861-4.8010
Columns89through96
-2.6789-2.0086-3.0031-5.3457-8.2703-10.7577-11.7996-10.6652
Columns97through104
-7.1059-1.44635.460812.409018.099521.430721.746118.9832
Columns105through112
13.68766.8892-0.1297-6.1215-10.1454-11.7627-11.1107-8.8369
Columns113through120
-5.9146-3.3828-2.0769-2.4143-4.2877-7.0931-9.8832-11.6121
Columns121through128
-11.4085-8.8128-3.91792.62039.709316.052420.435622.0000
>>dt=t
(2)-t
(1)
dt=
0.0079
>>f=1/dt
f=
127
>>x=fft(x);
>>F=x(1:
N/2+1);
>>f=f*(0:
N/2)/N;
>>plot(f,abs(F),'-*')
>>xlabel('Frequency');
>>ylabel('|F(k)|')
5、意见和建议
>>notebook-setup
WelcometotheutilityforsettinguptheMATLABNotebook
forinterfacingMATLABtoMicrosoftWord
Setupcomplete
>>notebook
Warning:
MATLABisnowanautomationserver
尊敬的老师,您好,一个学期过去,我们所选修的matlab课程接近尾声。
通过一个学期的学习,我觉得我初步掌握了如何运用matlab这一功能强大的数学软件,在做三次作业的过程中也发现了一些问题,
主要有两个:
1.matlab没有中文汉化,用起来不是十分方便,但是我相信,当我们的英语水平十分强大时,这也不是什么问题了。
2.输入错了一点它就会报错,所以工作的效率很难达到非常高的水平。
为此我想提出以下建议:
matlab本身可以学习微软的很对东西,例如自动纠错的功能,我觉得有了错误的提示,又有正确的提示,效率会更高,另外,就是期待汉化版本的出现,这样一来,它的表现一定会更出色。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- matlab 论文 doc
![提示](https://static.bdocx.com/images/bang_tan.gif)