资源预览内容
第1页 / 共24页
第2页 / 共24页
第3页 / 共24页
第4页 / 共24页
第5页 / 共24页
第6页 / 共24页
第7页 / 共24页
第8页 / 共24页
第9页 / 共24页
第10页 / 共24页
亲,该文档总共24页,到这儿已超出免费预览范围,如果喜欢就下载吧!
资源描述
1.求f(x)=4x6-x+x3-95在(0,100)上最大值、最小值,根。(使用函数fminbnd、roots)。程序: syms x f=4*x.6+x.3-x-95;x1=fminbnd( f ,0,100)x2=fminbnd(-4*x.6-x.3+x+95,0,100)x=x1;minf = eval(f)x=x2;maxf = eval(f)p=4 0 0 1 0 -1 -95;x=roots(p)结果:x1 = 0.4432x2 =99.9999minf =-95.3258maxf = 4.0000e+012x =1.6860 0.8525 + 1.4852i 0.8525 - 1.4852i -1.7050 -0.8431 + 1.4514i -0.8431 - 1.4514i2.求解常微分方程x=-x+x+1,x(0)=1,x(0)=0(使用函数dsolve)。程序:y=dsolve(D2x=-Dx+x+1,Dx(0)=0,x(0)=0)结果:y =exp(1/2*(5(1/2)-1)*t)*(1/10*5(1/2)+1/2)+exp(-1/2*(5(1/2)+1)*t)*(-1/10*5(1/2)+1/2)-1 3已知t=an2+bn,测得对应数据如下:(多项式插值interp1)t=0,20,40,60,80,100,120,140,160,183.5;n=0,1153,2045,2800,3466,4068,4621,5135,5619,6152;试求a和b的值。程序:n=0,1153,2045,2800,3466,4068,4621,5135,5619,6152;t=0,20,40,60,80,100,120,140,160,183.5;A=polyfit(n,t,2);a=A(1),b=A(2)结果:a =2.5083e-006b =0.01444请用梯形法、辛普森法分别计算积分值 (trapz、quad)程序:F = (x)sqrt(x.2+x+1);Q = quad(F,0,1)x=0:0.01:1;y=sqrt(x.2+x+1);q=trapz(x,y)结果:Q =1.3369q = 1.33695计算二重积分 (使用函数dblquad)程序:F = (x,y)x.2+y.2+x*y+2*x+y+1;Q = dblquad(F,0,2,0,1)结果:Q =11.33336矩阵M=1,2,6; 4,2,7; 8,9,3,求M的LU分解,QR分解,对角阵,特征值分解。(使用函数lu、qr、shol、eig)程序:A=1 2 3;4 5 6;7 8 9;L U=lu(a)Q R=qr(a)X D=eig(a)结果:L =1U =2.5083e-006Q =1R =2.5083e-006X =1D =2.5083e-0067 .a=3,A=4,b=a2,B=b2-1,c=a+A-2B,C=a+B+2c,求C。(使用函数solve)程序:D=solve(a=3,A=4,b=a2,B=b2-1,c=a+A-2*B,C=a+B+2*c);C=D.C结果:C =-2238用两种方法求解Ax=b的解。(A为四阶随机矩阵,b为四阶向量,自己创建)。程序:A=fix(50*rand(4) B=fix(50*rand(4,1) if det(A)=0 x=AB X=inv(A)*B else C=B; D=A C rref(B) end结果:A = 47 44 41 46 11 38 22 36 30 22 30 8 24 0 39 20B =46 45 20 44x =-0.8583 0.0430 1.3169 0.6621X =-0.8583 0.0430 1.3169 0.6621 9,用两种方法求函数的根,并求其极值与零点。程序:y1=x.3+(x-0.8).2/(x+1.25).3-5*(x+1/x);ezplot(y1,-100,100)xmin,ymin=fminbnd(y1,0,60);yminy2=-x.3-(x-0.8).2/(x+1.25).3+5*(x+1/x);xmax,y=fminbnd(y2,-50,-20);ymax=-y结果:ymin =-1.1489e+005ymax =-7.8999e+00310、f(x)的定义如下:,写一个matlab函数func1实现该函数。程序:function y=hs(x)x=input(请输入x值:x=);if x=0&x realmax。程序:n=1;while prod(1:n)1000, 在(xk, yk) 处亮一点(注意不要连线)可得所谓Henon引力线图。程序:x(1)=0; y(1)=0;for k=1:3000 x(k+1)=1+y(k)-1.4*x(k)2; y(k+1)=0.3*x(k);endplot(x(1001:3001),y(1001:3001),.);结果: 20某公司投资2000万元建成一条生产线。投产后,在时刻t 的追加成本和追加
网站客服QQ:2055934822
金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号