5、、min函数123456890maxi二9mini二0程序:◎勻・:♦回厂目目b目~1.0+■■1.1X熾1-a=[1,2,3$4$5,6,7,8,9_>0]2一maxi二max(a)3一mini二min(a)结果:1.求Fibonacci数列⑴大于4000的最小项。(2)5000之内的项数。解：(1)大于4000的最小项文件：網「昌I・二0
6、+
7、1.1IX12-3-4-□functiont=fibonacclf=[l,1]:n=2：□while1f=[f,f(n-1)+f(n)]:n=length(f):iff(n)>4000break:10--en
8、d11-t=f(n):12一return9-end结果:arts=4181(2)5()00之内的项数文件：-1.0+1.1X+回「回回b回Efunctionn=fibonacc2f=[l,1]:n=2;Ewhile110f=[f,f(n-l)+f(n)];n=length(f);iff(n)>4000break:endend1112n=n-1;return结果:ans=18第四讲习题1.解方程组Ax=b,分别用求逆解法与直接解法求其解。_2-1-120-1000■0~10A=0-12-10,b=000-12-10000-12__0解：(1)逆解法文件：
9、馆謂
10、・1.0[+
11、专i・i[x
12、癖疵
13、Q1一A=[2,-l,0,0,0;-1,2,-1,0,0;0,-1,2,-1,0:0,0,-1,2,-1:0,0,0,-1,2]:2-b=[l;0：0;0;0]：3一x=irLy(A)*b结果:0.83330.66670.50000.33330.1667A»(2)直接法文件：•昌「昌
14、-
15、1・0]+
16、专
17、1・1[x
18、癖疫
19、Q1-A二[2,-1,0,0,0：-1,2,-1,0,0;0,-1,2,-1,0;0,0,-1,2,-1;0,0,0,-1,2]:2-b=[l;0;0：0;0];3—x=Ab结果:0.8333
20、0.66670.50000.33330.16672.编一个m程序，求N阶方阵A的行列式的值。解:文件:1.0+
21、^-1.1X
22、癖疾Q结果:♦目「回回卜回1□functions=mydet(A)2-n=size(A):3-if(n(l)^=n(2))5_6—u8—9-return;endn=n(1);if(n==1)s=A;return;ends=0;fori=l:ns=s-(-1)"i*A(1,i)*mydet(A(2:m[1:i-1i+1:n]));end>>A=mydet([1,3;5,7])1.今有多项式P1(x)=x4-2x+1,P2(x)=x
23、2+4x-0.5,要求先求得P(x)二Pl(x)+P2(x),然后计算xi=0.2*i各点上的P(xi)(i=0,l,2,…,5)值。解：文件：[增頃
24、・1.0+
25、号
26、1・1x
27、囁1一pl=[l00-21]:2-p2=[l4-0.5]:3一plx=poly2sym(pl);4一p2x=poly2sym(p2);5一p=plx+p2x6一x=0:57-x.*4+2*x+l/2+x.*2
28、结果:p=x-4+x*2+2林+1/2012345ans=0.50004.500024.5000
29、9&5000280.5000660.50002.试编一个m程序，将一维数
30、组x中的N个数按颠倒的次序重新存储。如N=5,原来x为:x=[13579]而经过颠倒处理后x中数据的次序应该为：X=l97531]解：文件：:♦回「目•目b目<%1-x=[13579]：2一[%n]=size(x):3一y=zeros(size(x)):4-Efori=l:n5-y(i)=x(n-i+l):6--end结果:97531A».•I•••第五章习题1・在[02兀]范围内绘制二维曲线图y=sin(x)*cos(5x)o解：文件：["昌瑁
31、-]o[+
32、壬1・1
33、癖疾
34、Q1-x=0:pi/200:2*pi;2-y=sin(x).*cos(5*x)
35、3-plot(x^y)结果:2.把自己姓名中的各个拼音字母（小写）所对应的ASCII值转换为二
