资源描述：
《matlab结课论文.doc》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、北京林业大学全校公共选修课结课作业MATLAB程序应用姓名：课程序号：51班级：电气10-1学号：1010541182011年12月2011-2012学年MATLAB程序应用结课作业-51-乔林翰1、（指选题1）有一天小猴摘下了若干个桃子，当即吃掉了一半，还觉得不过隐，又多吃了2个。第二天接着吃下剩下的一半，又多吃了2个。以后每天都是吃掉尚存的一半零两个。到第十天早上，小猴准备吃桃子时，看到只剩下n个桃子了。问小猴第一天共摘下多少个桃子？1)写出求解小猴第一天共摘下多少个桃子的数学表达式2)根据确定的数学表达式，利用mat
2、lab编写出一个函数文件用于求解3)编写一个命令文件，完成输入一个n值，然后调用上述函数求解，并将结果输出显示出来【思路和方法】假设猴子第一天吃掉……因此本题可以转化利用matlab求解级数的问题……（正文格式：中文宋体，西文timesNewRoman，1.25倍行距）【源程序】inputn%输入n【运行结果】>>display('Pleaseinputthevalueofn');n=input('n=');fori=1:9;n=(n+2)*2;endm=nPleaseinputthevalueofnn=112011-20
3、12学年MATLAB程序应用结课作业-51-乔林翰2、（指选题2）1、分别用矩阵求逆、矩阵除法、矩阵分解求下面线性方程组的解；矩阵除法：A=1-23-401-1113010-731>>b=[4,-3,1,-3]b=4-31-3>>x=A/bx=0.7143-0.2000-0.22860.6000矩阵求逆A=[1,-2,3,-4;0,1,-1,1;1,3,0,1;0,-7,3,1]A=1-23-401-11112011-2012学年MATLAB程序应用结课作业-51-乔林翰13010-731>>b=[4,-3,1,-3]b=
4、4-31-3x=inv(A)*bx=0.7143-0.2000-0.22860.6000矩阵分解法（LU分解）：[L,U]=lu(A);>>x=U/(L/b)x=0.7143-0.2000-0.22860.60004、某乡镇企业某产品年销售额与所获利润从2000年到2010年的数据如下：年度20002001200220032004200520062007200820092010销售额x6.17.59.410.714.617.421.124.429.832.934.3纯利润4.56.48.38.49.711.513.715.
5、417.720.522.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)112011-2012学年MATLAB程序应用结课作业-51-乔林翰p=0.56391.9096xi=6.1:34.3;>>yi=polyval(p,xi)yi=Columns
6、1through85.34915.91306.47687.04077.60458.16848.73229.2961Columns9through169.860010.423810.987711.551512.115412.679213.243113.8069Columns17through2414.370814.934615.498516.062316.626217.190017.753918.3177Columns25through2918.881619.445420.009320.573121.1370plot(x,y
7、,':o',xi,yi,'-*')112011-2012学年MATLAB程序应用结课作业-51-乔林翰112011-2012学年MATLAB程序应用结课作业-51-乔林翰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);112011-2012学年MATLAB程序应用结课作业-51-乔林翰4、（自选题）我们是电气工程及其自动化专业，对电路分析有时会需要进
8、行傅里叶变换和拉普拉斯变换，通过学习这门课，我们发现matlab可以实现对一个函数进行快速傅氏变换（fft），以下为具体操作：对函数x(t)=12sin(10πt+π/2)+10(cos20πt)进行fft：N=128N=128>>T=1T=1>>t=linspace(0,T,N);>>x=12*sin
