matlab结课论文.doc.docx

北京林业大学全校公共选修课结课作业
MATLAB程序应用
姓名:
课程序号:
51
班级:
电气10-1
学号:
101054118
2011年12月
1、(指选题1)有一天小猴摘下了若干个桃子,当即吃掉了一半,还觉得不过隐,又多吃了2个。第二天接着吃下剩下的一半,又多吃了2个。以后每天都是吃掉尚存的一半零两个。到第十天早上,小猴准备吃桃子时,看到只剩下n个桃子了。问小猴第一天共摘下多少个桃子?
写出求解小猴第一天共摘下多少个桃子的数学表达式
根据确定的数学表达式,利用matlab编写出一个函数文件用于求解
编写一个命令文件,完成输入一个n值,然后调用上述函数求解,并将结果输出显示出来
【思路和方法】
假设猴子第一天吃掉……
因此本题可以转化利用matlab求解级数的问题……(正文格式:中文宋体,西文times New Roman,)
【源程序】
input n %输入n
【运行结果】
>> display('Please input the value of n');
n=input('n=');
for i=1:9;
n=(n+2)*2;
end
m=n
Please input the value of n
n=
2、(指选题2)
1、分别用矩阵求逆、矩阵除法、矩阵分解求下面线性方程组的解;
矩阵除法:
A =
1 -2 3 -4
0 1 -1 1
1 3 0 1
0 -7 3 1
>> b=[4,-3,1,-3]
b =
4 -3 1 -3
>> x=A/b
x =

-
-

矩阵求逆
A=[1,-2,3,-4;0,1,-1,1;1,3,0,1;0,-7,3,1]
A =
1 -2 3 -4
0 1 -1 1
1 3 0 1
0 -7 3 1
>> b=[4,-3,1,-3]
b =
4 -3 1 -3
x=inv(A)*b
x =

-
-

矩阵分解法(LU分解):
[L,U]=lu(A);
>> x=U/(L/b)
x =

-
-

4、某乡镇企业某产品年销售额与所获利润从2000年到2010年的数据如下:
年度
2000
2001
2002
2003
2004
2005
2006
2007
2008
2009
2010
销售额x











纯利润











试求对的经验回归直线方程,并做回归分析,对回归的参数进行区间估计;
x=[,,,,,,,,,,];
>> y=[,,,,,,,,,,];
>> p=polyfit(x,y,1)
p =

xi=:;
>> yi=polyval(p,xi)
yi =
Columns 1 through 8

Columns 9 through 16

Columns 17 through 24

Columns 25 through 29

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、(自选题)
我们是电气工程及其自动化专业,对电路分析有时会需