MATLAB軟件實(shí)習(xí)課程實(shí)習(xí)報(bào)告

精選優(yōu)質(zhì)文檔-----傾情為你奉上精選優(yōu)質(zhì)文檔-----傾情為你奉上專心---專注---專業(yè)專心---專注---專業(yè)精選優(yōu)質(zhì)文檔-----傾情為你奉上專心---專注---專業(yè)《MATLAB軟件實(shí)習(xí)》課程實(shí)習(xí)報(bào)告姓名：專業(yè)：數(shù)學(xué)與應(yīng)用數(shù)學(xué)班級：學(xué)號：《MATLAB軟件實(shí)習(xí)》課程實(shí)習(xí)報(bào)告姓名專業(yè)數(shù)學(xué)與應(yīng)用數(shù)學(xué)班級09級一班學(xué)號實(shí)習(xí)地點(diǎn)305、107實(shí)習(xí)時(shí)間第二周實(shí)習(xí)章節(jié)2.1MATLAB程序設(shè)計(jì)語言基礎(chǔ)、2.2基本數(shù)學(xué)運(yùn)算實(shí)習(xí)內(nèi)容：1，tic,A=rand(500);>>B=inv(A);>>norm(A*B-eye(500)),tocans=9.8984e-0122．symsxs;>>f=x^5+3*x^4+4*x^3+2*x^2+3*x+6;>>f1=subs(f,x,(s-1)/(s+1))f1=(s-1)^5/(s+1)^5+3*(s-1)^4/(s+1)^4+4*(s-1)^3/(s+1)^3+2*(s-1)^2/(s+1)^2+3*(s-1)/(s+1)+63.>>A=[1,2,3,4;4,3,2,1;2,3,4,1;3,2,4,1]A=1234432123413241>>B=[1+4j,2+3j,3+2j,4+1j;4+1j,3+2j,2+3j,1+4j;2+3j,3+2j,4+1j,1+4j;3+2j,2+3j,4+1j,1+4j]B=1.0000+4.0000i2.0000+3.0000i3.0000+2.0000i4.0000+1.0000i4.0000+1.0000i3.0000+2.0000i2.0000+3.0000i1.0000+4.0000i2.0000+3.0000i3.0000+2.0000i4.0000+1.0000i1.0000+4.0000i3.0000+2.0000i2.0000+3.0000i4.0000+1.0000i1.0000+4.0000i>>A(5,6)=5A=1234004321002341003241000000054.>>A=magic(8),B=A(2:2:end,:)A=64236160675795554121351501617474620214342244026273736303133323435292838392541232244451918484915145253111056858595462631B=95554121351501640262737363031334123224445191848858595462631>>A=magic(8);B=A(2:2:end,:)B=95554121351501640262737363031334123224445191848858595462631實(shí)習(xí)總結(jié)：通過本次課的學(xué)習(xí)初步掌握了matlab中一些最基本語句的命令，以及一些特殊矩陣。例如，冒號表達(dá)式，子矩陣提取等。在第二節(jié)還學(xué)到Matlab中矩陣的代數(shù)運(yùn)算，比較運(yùn)算等的命令，最后還學(xué)了化簡函數(shù)、替換函數(shù)、多項(xiàng)式展開、因式分解等函數(shù)。這兩節(jié)的學(xué)習(xí)使matlab有了最初步的認(rèn)識。學(xué)生簽名：《MATLAB軟件實(shí)習(xí)》課程實(shí)習(xí)報(bào)告姓名專業(yè)數(shù)學(xué)與應(yīng)用數(shù)學(xué)班級09級一班學(xué)號實(shí)習(xí)地點(diǎn)305、107實(shí)習(xí)時(shí)間第三周實(shí)習(xí)章節(jié)2.3MATLAB語言的流程結(jié)構(gòu)、2.4函數(shù)編寫與調(diào)試實(shí)習(xí)內(nèi)容：5、functiony=xuexiaopei5(x,D,h)ifx>D,y=h;elseif-D<=x<=D,y=h/(D*x);elsex<-D,y=-h;end6、sum(sym(2).^[1:63])ans=8、v=[123];A=hankel(v)A=123230300、>>H=myhankel(v)H=12239、functiony=xuexiaopei9(n)ifround(n)==n&n>=1ifn>=3y=xuexiaopei9(n-1)+xuexiaopei9(n-2);elsey=1;endelseerror('nmustbepositiveinteger.')endans=55>>tic,xuexiaopei9(20),tocans=6765Elapsedtimeis0.seconds.>>tic,a=[11];>>fori=3:30,a(i)=a(i-1)+a(i-2);end,a(20),tocans=6765Elapsedtimeis69.seconds.10、M=[15192018;24222938;44283610;12151628];>>iM=inv(M)iM=-0.1547-0.08510.07010.19000.0748-0.22170.02930.24230.14690.2774-0.0808-0.4420-0.0577-0.00330.00040.0771>>A=diag([1234]);B=hankel([1234]);C=diag([4321]);>>iM1=xuexiaopei10(A,B,C)iM1=0.0553-0.03890.00170.0041-0.03890.0555-0.0210-0.00210.0017-0.02100.0328-0.01370.0041-0.0021-0.01370.0244>>M1=sym(M);iM0=inv(M1)iM0=[-1613/10429,-888/10429,731/10429,1981/10429][780/10429,-2312/10429,306/10429,2527/10429][1532/10429,2893/10429,-843/10429,-4610/10429][-602/10429,-34/10429,9/20858,804/10429]>>norm(double(iM0)-iM)ans=1.4082e-016>>norm(double(iM0)-iM1)ans=0.6859實(shí)習(xí)總結(jié)：了解了MATLAB中的三大結(jié)構(gòu)即循環(huán)結(jié)構(gòu)、轉(zhuǎn)移結(jié)構(gòu)、開關(guān)結(jié)構(gòu)。學(xué)生簽名：《MATLAB軟件實(shí)習(xí)》課程實(shí)習(xí)報(bào)告姓名專業(yè)數(shù)學(xué)與應(yīng)用數(shù)學(xué)班級09級一班學(xué)號實(shí)習(xí)地點(diǎn)305、107實(shí)習(xí)時(shí)間第4周實(shí)習(xí)章節(jié)2.5二維圖形繪制、2.6三維圖形表示實(shí)習(xí)內(nèi)容：12、t=[0,120,240,0]*pi/180;x=[];y=[];fori=0;5:360tt=i*pi/180;x=[x;cos(tt+t)];y=[y;sin(tt+t)];endplot(x,y,'b'),axis('square')13、%sin1/t等間距圖形t=-1:0.03:1;y=sin(1./t);plot(t,y)%sin1/t的不等間距圖形t=[-1:0.01:-0.2,-0.21:0.003:0.41,0.4:0.1:1];y=sin(1./t);plot(t,y)15、ezplot('x^2+y^2-3*x*y^2');holdonezplot('x^3-x^2-y^2+y')16、[x,y]=meshgrid(-pi:0.1:2*pi);z=sin(x.*y);surf(x,y,z),figure;contour3(x,y,z,30)14、%(1)的極坐標(biāo)圖t=0:0.01:2*pi;y1=1.0013.*t.^2;polar(t,y1)%(2)的極坐標(biāo)圖t=0:0.01:2*pi;y2=cos(7.*t./2);polar(t,y2)%(3)的極坐標(biāo)圖t=0:0.01:2*pi;y3=sin(t./t);polar(t,y3)17、[x,y]=meshgrid(-1:.1:1);z=sin(x.*y);i=find(x.^2+y.^2<=0.5^2);z(i)=NaN;surf(x,y,z)實(shí)習(xí)總結(jié)：本次實(shí)習(xí)我學(xué)會(huì)了基本的圖形的繪制及三維圖形的表示。學(xué)生簽名：《MATLAB軟件實(shí)習(xí)》課程實(shí)習(xí)報(bào)告姓名專業(yè)數(shù)學(xué)與應(yīng)用數(shù)學(xué)班級09級一班學(xué)號實(shí)習(xí)地點(diǎn)D樓305、107機(jī)房實(shí)習(xí)時(shí)間第5周實(shí)習(xí)章節(jié)3.1極限問題的解析解、3.2函數(shù)導(dǎo)數(shù)的解析解、3.3積分問題的解析解實(shí)習(xí)內(nèi)容：1、%(1)symsx;f=(3^x+9^x)^(1/x);L=limit(f,x,inf)%(2)symsx;f=(x+2)^(x+2)*(x+3)^(x+3)/(x+5)^(2*x+5);L=limit(f,x,inf)L=92、%(1)symsxy;f=(x^2*y+x*y^3)/(x+y)^3;L1=limit(limit(f,x,-1),y,2)%(2)symsxy;f=x*y/(sqrt(x*y+1)-1);L2=limit(limit(f,x,0),y,0)%(3)symsxy;f=(1-cos(x^2+y^2))/(x^2+y^2)*exp(x^2+y^2);L3=limit(limit(f,x,0),y,0)L1=-6L3、%(1)symsx;f=sqrt(x*sin(x)*sqrt(1-exp(x)));y1=diff(f,x)%(2)symsxa;f=(1-sqrt(cos(a*x)))/(x*(1-cos(sqrt(a*x))));y2=diff(f,x)%(3)symsx;f=atan(y/x)/log(x^2+y^2);y3=diff(f,x)%(4)symsxan;f=-1/(n*a)*log((x^n+a)/x^n);y4=diff(f,x)y1=1/2/(x*sin(x)*(1-exp(x))^(1/2))^(1/2)*(sin(x)*(1-exp(x))^(1/2)+x*cos(x)*(1-exp(x))^(1/2)-1/2*x*sin(x)/(1-exp(x))^(1/2)*exp(x))y2=1/2/cos(a*x)^(1/2)*sin(a*x)*a/x/(1-cos((a*x)^(1/2)))-(1-cos(a*x)^(1/2))/x^2/(1-cos((a*x)^(1/2)))-1/2*(1-cos(a*x)^(1/2))/x/(1-cos((a*x)^(1/2)))^2*sin((a*x)^(1/2))/(a*x)^(1/2)*ay3=-y/x^2/(1+y^2/x^2)/log(x^2+y^2)-2*atan(y/x)/log(x^2+y^2)^2*x/(x^2+y^2)y4=-1/n/a*(n/x-(x^n+a)/(x^n)*n/x)/(x^n+a)*x^n4、symsx;f=sqrt((x-1)*(x-2)/(x-3)/(x-4));y=diff(f,x,4)y=-15/16/((x-1)*(x-2)/(x-3)/(x-4))^(7/2)*((x-2)/(x-3)/(x-4)+(x-1)/(x-3)/(x-4)-(x-1)*(x-2)/(x-3)^2/(x-4)-(x-1)*(x-2)/(x-3)/(x-4)^2)^4+9/4/((x-1)*(x-2)/(x-3)/(x-4))^(5/2)*((x-2)/(x-3)/(x-4)+(x-1)/(x-3)/(x-4)-(x-1)*(x-2)/(x-3)^2/(x-4)-(x-1)*(x-2)/(x-3)/(x-4)^2)^2*(2/(x-3)/(x-4)-2*(x-2)/(x-3)^2/(x-4)-2*(x-2)/(x-3)/(x-4)^2-2*(x-1)/(x-3)^2/(x-4)-2*(x-1)/(x-3)/(x-4)^2+2*(x-1)*(x-2)/(x-3)^3/(x-4)+2*(x-1)*(x-2)/(x-3)^2/(x-4)^2+2*(x-1)*(x-2)/(x-3)/(x-4)^3)-3/4/((x-1)*(x-2)/(x-3)/(x-4))^(3/2)*(2/(x-3)/(x-4)-2*(x-2)/(x-3)^2/(x-4)-2*(x-2)/(x-3)/(x-4)^2-2*(x-1)/(x-3)^2/(x-4)-2*(x-1)/(x-3)/(x-4)^2+2*(x-1)*(x-2)/(x-3)^3/(x-4)+2*(x-1)*(x-2)/(x-3)^2/(x-4)^2+2*(x-1)*(x-2)/(x-3)/(x-4)^3)^2-1/((x-1)*(x-2)/(x-3)/(x-4))^(3/2)*((x-2)/(x-3)/(x-4)+(x-1)/(x-3)/(x-4)-(x-1)*(x-2)/(x-3)^2/(x-4)-(x-1)*(x-2)/(x-3)/(x-4)^2)*(-6*(x-1)*(x-2)/(x-3)/(x-4)^4-6/(x-3)/(x-4)^2-6/(x-3)^2/(x-4)+6*(x-2)/(x-3)^3/(x-4)+6*(x-2)/(x-3)^2/(x-4)^2+6*(x-2)/(x-3)/(x-4)^3+6*(x-1)/(x-3)^3/(x-4)+6*(x-1)/(x-3)^2/(x-4)^2-6*(x-1)*(x-2)/(x-3)^4/(x-4)-6*(x-1)*(x-2)/(x-3)^3/(x-4)^2+6*(x-1)/(x-3)/(x-4)^3-6*(x-1)*(x-2)/(x-3)^2/(x-4)^3)+1/2/((x-1)*(x-2)/(x-3)/(x-4))^(1/2)*(-24*(x-2)/(x-3)/(x-4)^4-24*(x-1)/(x-3)/(x-4)^4-24*(x-2)/(x-3)^4/(x-4)+24/(x-3)^3/(x-4)+/(x-3)/(x-4)^3+24/(x-3)^2/(x-4)^2+24*(x-1)*(x-2)/(x-3)^2/(x-4)^4+24*(x-1)*(x-2)/(x-3)/(x-4)^5-24*(x-2)/(x-3)^2/(x-4)^3-24*(x-1)/(x-3)^4/(x-4)-24*(x-1)/(x-3)^3/(x-4)^2-24*(x-2)/(x-3)^3/(x-4)^2+24*(x-1)*(x-2)/(x-3)^5/(x-4)+24*(x-1)*(x-2)/(x-3)^4/(x-4)^2-24*(x-1)/(x-3)^2/(x-4)^3+24*(x-1)*(x-2)/(x-3)^3/(x-4)^3)5、%對分子，分母分別求極限symsx;f1=log(1+x)*log(1-x)-log(1-x^2);f2=x^4;y1=limit(f1,x,0);y2=limit(f2,x,0);y0=y1/y2%用L’Hopital法則symsx;f=(log(1+x)*log(1-x)-log(1-x^2))/x^4;y=limit(f,x,0)y0=NaNy=1/126、%一階導(dǎo)symst;x=log(cos(t));y=cos(t)-t*sin(t);f=xuexiaopei(y,x,t,1);[n,d]=numden(f);F1=simple(n)/simple(d)%二階導(dǎo)及t=pi/3f=xuexiaopei(y,x,t,3);[n,d]=numden(f);F2=simple(n)/simple(d)F1=(2*sin(t)+t*cos(t))*cos(t)/sin(t)F2=cos(t)*(4*cos(t)^4*sin(t)+cos(t)^5*t-9*cos(t)^2*sin(t)+2*sin(t)-2*cos(t)^3*t+4*t*cos(t))/sin(t)^5f1=1.91467、%uxysymsxy;u=(cos(x/y)^1/2)^(-1);u1=diff(diff(u,x),y)%uyxu2=diff(diff(u,y),x)u1=-4/cos(x/y)^3*sin(x/y)^2/y^3*x-2/cos(x/y)*x/y^3-2/cos(x/y)^2*sin(x/y)/y^2u2=-4/cos(x/y)^3*sin(x/y)^2/y^3*x-2/cos(x/y)*x/y^3-2/cos(x/y)^2*sin(x/y)/y^28、symsxy;u=1/(y^2-x^2);u1=diff(diff(u,x),y)u1=-8/(y^2-x^2)^3*x*y9、symst;f=int(exp(-t^2),t,0,x*y);f1=diff(diff(f,x),x);f2=diff(diff(f,x),y);f3=diff(diff(f,y),y);f0=x/y*f1-2*f2+f3f0=2*x^2*y^2*exp(-x^2*y^2)-2*exp(-x^2*y^2)-2*x^3*y*exp(-x^2*y^2)10、symsxyz;f=[3*x+exp(y)*z;x^3+y^2*sin(z)];H=jacobian(jacobian(f,[x,y,z]),[x,y,z])H=[0,0,0][6*x,0,0][0,exp(y)*z,exp(y)][0,2*sin(z),2*y*cos(z)][0,exp(y),0][0,2*y*cos(z),y^2*sin(z)]11、%(1)symsxa;I1=int((3*x^2+a)/x^2/(x^2+a)^2,x)%(2)symsx;I2=int(sqrt((x*(x+1)))/(sqrt(x)+sqrt(1+x)),x)%(3)symsxab;I3=int(x*exp(a*x)*cos(b*x),x)%(4)symsxabc;I4=int(exp(a*x)*sin(b*x)*sin(c*x),x)I1=x/a/(x^2+a)-1/a/xI2=2/15*(x*(1+x))^(1/2)/(1+x)^(1/2)*x*(3*x+5)-2/15*(x*(1+x))^(1/2)/x^(1/2)*(1+x)*(-2+3*x)I3=(a/(a^2+b^2)*x-(a^2-b^2)/(a^2+b^2)^2)*exp(a*x)*cos(b*x)-(-b/(a^2+b^2)*x+2*a*b/(a^2+b^2)^2)*exp(a*x)*sin(b*x)I4=1/2*a/(a^2+(b-c)^2)*exp(a*x)*cos((b-c)*x)-1/2*(-b+c)/(a^2+(b-c)^2)*exp(a*x)*sin((b-c)*x)-1/2*a/(a^2+(b+c)^2)*exp(a*x)*cos((b+c)*x)+1/2*(-b-c)/(a^2+(b+c)^2)*exp(a*x)*sin((b+c)*x)12、%(1)symsx;I1=int(cos(x)/sqrt(x),x,0,inf)%(2)symsx;I2=int((1+x^2)/(1+x^4),x,0,1)I1=1/2*2^(1/2)*pi^(1/2)I2=1/4*2^(1/2)*pi13、symsx;f=exp(-5*x)*sin(3*x+pi/3);f1=subs(f,x,t+x);R=int(f*f1,x,0,t)R=1/17680*(-31200*sin(t)^2*3^(1/2)*cos(t)^5-1170*sin(t)^2*3^(1/2)*cos(t)+24960*sin(t)^2*3^(1/2)*cos(t)^7+100*sin(t)*3^(1/2)+10875*sin(t)*3^(1/2)*cos(t)^2+402*cos(t)+195*sin(t)+49920*cos(t)^8*sin(t)-425*sin(t)^3*3^(1/2)+56160*3^(1/2)*cos(t)^7-87360*cos(t)^6*sin(t)+46800*cos(t)^4*sin(t)-42120*3^(1/2)*cos(t)^5-24960*3^(1/2)*cos(t)^9+18289*cos(t)^3-71900*cos(t)^5+1700*sin(t)^3*3^(1/2)*cos(t)^2-41600*cos(t)^9-83200*sin(t)*3^(1/2)*cos(t)^8+10920*sin(t)^2*3^(1/2)*cos(t)^3-585*cos(t)*3^(1/2)+93600*cos(t)^7-52000*sin(t)^2*cos(t)^5+16500*cos(t)^3*sin(t)^2-675*cos(t)*sin(t)^2+*cos(t)^7*sin(t)^2-7800*sin(t)*cos(t)^2+11310*cos(t)^3*3^(1/2)+*sin(t)*3^(1/2)*cos(t)^6-76300*sin(t)*3^(1/2)*cos(t)^4+1300*exp(t)^10*sin(t)*3^(1/2)*cos(t)^2+195*exp(t)^10*sin(t)-3627*exp(t)^10*cos(t)+4836*exp(t)^10*cos(t)^3-325*exp(t)^10*sin(t)*3^(1/2)+780*exp(t)^10*cos(t)^3*3^(1/2)-585*exp(t)^10*cos(t)*3^(1/2)-780*exp(t)^10*sin(t)*cos(t)^2)/exp(t)^1514、symsxa;f=cos(a*x)/(1+x^2);t=[0:0.1;pi];y=[];forn=tf1=subs(f,a,n);I=int(f1,x,0,inf);y=[y,double(I)]endplot(t,y)15、symstab;f=sin(t);I1=int(f,t,a,b)I2=int(f,t,b,a)I1=-cos(b)+cos(a)I2=-cos(a)+cos(b)16、%(1)symsxy;f=sqrt(4-x^2-y^2);I1=int(int(f,y,0,sqrt(4-x^2)),x,0,2)%(2)symsxyz;f=x*y*z;I2=int(int(int(f,z,0,3-x-y),y,0,3-x),x,0,3)%(3)symsxyz;f=z*(x^2+y^2);I3=int(int(int(f,z,0,sqrt(4-x^2-y^2)),y,0,sqrt(4-x^2)),x,0,2)%(4)symsxyzu;f=x*y*z*u*exp(6-x^2-y^2-z^2-u^2);I4=int(int(int(int(f,u,0,z),z,0,y),y,0,x),x,0,1)I1=4/3*piI2=81/80I3=4/3*piI4=1/384*exp(2)-1/96*exp(3)+1/64*exp(4)-1/96*exp(5)+1/384*exp(6)實(shí)習(xí)總結(jié)：在這節(jié)課上讓我了解到極限問題的解析解、函數(shù)導(dǎo)數(shù)的解析解、積分問題的解析解的相關(guān)問題，并且學(xué)會(huì)了初步的運(yùn)用。學(xué)生簽名：《MATLAB軟件實(shí)習(xí)》課程實(shí)習(xí)報(bào)告姓名專業(yè)數(shù)學(xué)與應(yīng)用數(shù)學(xué)班級09級一班學(xué)號實(shí)習(xí)地點(diǎn)305、107實(shí)習(xí)時(shí)間第七周實(shí)習(xí)章節(jié)3.4函數(shù)的級數(shù)展開與級數(shù)求和問題求解、3.5曲線積分與曲面積分的計(jì)算內(nèi)容：16>symsxy;f=sqrt(4-x^2-y^2);I1=int(int(f,y,0,sqrt(4-x^2)),x,0,2)I1=4/3*pi>>symsxyz;f=x*y*z;I2=int(int(int(f,z,0,3-x-y),y,0,3-x),x,0,3)I2=81/80>>symsxyz;f=z*(x^2+y^2);I3=int(int(int(f,z,0,sqrt(4-x^2-y^2)),y,0,sqrt(4-x^2)),x,0,2)I3=4/3*pi>>symsxyzu;f=x*y*z*u*exp(6-x^2-y^2-z^2-u^2);I4=int(int(int(int(f,u,0,z),z,0,y),y,0,x),x,0,1)I4=1/384*exp(2)-1/96*exp(3)+1/64*exp(4)-1/96*exp(5)+1/384*exp(6)function[dy,dx]=diff_ctr(y,Dt,n)yx1=[y00000];yx2=[0y0000];yx3=[00y000];yx4=[000y00];yx5=[0000y0];yx6=[00000y];switchncase1dy=(-diff(yx1)+7*diff(yx2)+7*diff(yx3)-diff(yx4))/(12*Dt);L0=3;case2dy=(-diff(yx1)+15*diff(yx2)-15*diff(yx3)+diff(yx4))/(12*Dt^2);L0=3;case3dy=(-diff(yx1)+7*diff(yx2)-6*diff(yx3)-6*diff(yx4)+7*diff(yx5)-diff(yx6))/(8*Dt^3);L0=5;case4dy=(-diff(yx1)+11*diff(yx2)-28*diff(yx3)+28*diff(yx4)-11*diff(yx5)+diff(yx6))/(6*Dt^4);L0=5;enddy=dy(L0+1:end-L0);dx=([1:length(dy)]+L0-2-(n>2))*Dt;實(shí)習(xí)總結(jié)：這次上級實(shí)習(xí)我了解了Taylor冪級數(shù)展開、Fourier級數(shù)展開、級數(shù)求和的計(jì)算、序列求積問題以及曲線積分及MATLA求解還有曲面積分問題。學(xué)生簽名：《MATLAB軟件實(shí)習(xí)》課程實(shí)習(xí)報(bào)告姓名專業(yè)數(shù)學(xué)與統(tǒng)計(jì)學(xué)院班級09級一班學(xué)號實(shí)習(xí)地點(diǎn)107、305實(shí)習(xí)時(shí)間第八周實(shí)習(xí)章節(jié)3.6數(shù)值微分問題、3.7數(shù)值積分問題實(shí)習(xí)內(nèi)容：22、>>symskntheta;s=symsum(cos(k*theta),k,1,n);s1=sin(n*theta/2)*cos((n+1)*theta/2)/sin(theta/2);s2=simple(symsum((s-s1),n,1,inf))s2=025、>>f=@(t)(pi-t).^(1/4).*exp(-t).*sin(3.*t+1);y=quad(f,0,pi)y=0.341426、x=[0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.,1.1,1.2];y=[0,2.2077,3.2058,3.4435,3.241,2.8164,2.311,1.8101,1.3602,0.98172,0.6797,0.4473,0.27684];[dy1,dx1]=diff_ctr(y,x(2)-x(1),1);[dy2,dx2]=diff_ctr(y,x(2)-x(1),2);[dy3,dx3]=diff_ctr(y,x(2)-x(1),3);[dy4,dx4]=diff_ctr(y,x(2)-x(1),4);plot(dx1+x(1),dy1,'+',dx2+x(1),dy2,':',dx3+x(1),dy3,'*',dx4+x(1),dy4,'x')281、[x,y]=meshgrid(0:.1:2,0:.1:2);z=4-x.^2-y.^2;27>>(-x.^2-y.^2-z.^2-x.*z);I2=triplequad(f,0,2,0,2,0,2,1e-7,@quadl)I2=0.2078[fx,fy]=gradient(z);fx=fx/0.2;fy=fy/0.2;contour(x,y,z,30);holdon;quiver(x,y,fx,fy)[x,y]=meshgrid(0:.1:2,0:.1:2);z=4-x.^2-y.^2;[fx,fy]=gradient(z);fx=fx/0.2;fy=fy/0.2;zx=-2*x;zy=-2*y;surf(x,y,fx-zx);figure;surf(x,y,fy-zy);實(shí)習(xí)總結(jié)：這次的上機(jī)實(shí)習(xí)讓我學(xué)會(huì)了數(shù)值微分法、中心差分法等的數(shù)值微分問題，以及，由給定數(shù)據(jù)進(jìn)行梯形求積、單變量數(shù)值積分問題的求解等等的數(shù)值積分問題。學(xué)生簽名：《MATLAB軟件實(shí)習(xí)》課程實(shí)習(xí)報(bào)告姓名專業(yè)數(shù)學(xué)與統(tǒng)計(jì)學(xué)院班級09級一班學(xué)號實(shí)習(xí)地點(diǎn)實(shí)習(xí)時(shí)間第九周實(shí)習(xí)章節(jié)4.1特殊矩陣的輸入、4.2矩陣基本分析實(shí)習(xí)內(nèi)容：1、>>J=diag([-5-5-5-5-5])+diag([1111],1)J=-510000-510000-510000-510000-53、>>A=ones(6);B=sym(A);C=A*BC=[6,6,6,6,6,6][6,6,6,6,6,6][6,6,6,6,6,6][6,6,6,6,6,6][6,6,6,6,6,6][6,6,6,6,6,6]4>>A=[5765165;2310014;6420644;3963662;10760077;7244077;4867217];B=[3550123;3254625;1211346;3515212;4101201;-3-4-737812;1-107-6815];A1=sym(A)B1=sym(B)A1=[5,7,6,5,1,6,5][2,3,1,0,0,1,4][6,4,2,0,6,4,4][3,9,6,3,6,6,2][10,7,6,0,0,7,7][7,2,4,4,0,7,7][4,8,6,7,2,1,7]B1=[3,5,5,0,1,2,3][3,2,5,4,6,2,5][1,2,1,1,3,4,6][3,5,1,5,2,1,2]Columns1through50.5186-0.26120.3128-0.38590.03910.72980.17350.2343-0.3007-0.1059-1.3574-0.2625-0.58260.68760.34760.6226-0.11030.2362-0.3215-0.1793-0.3479-0.05010.01310.1821-0.01430.21340.14480.00620.0213-0.1040-0.52070.2799-0.18810.2413-0.0010Columns6through7-0.2886-0.0402-0.3969-0.16560.54490.3635-0.1901-0.05540.14350.08850.0484-0.18930.29360.10087、>>A=[5765165;2310014;6420644;3963662;10760077;7244077;4867217];B=[3550123;3254625;1211346;3515212;4101201;-3-4-737812;1-107-6815];p1=poly(A);p2=poly(B);[V1,D1]=eig(A)[V2,D2]=eig(B)polyvalm(p1,A)polyvalm(p2,B)V1=Columns1through30.41660.0385-0.1349i0.0385+0.1349i28.6796000-1.6337+6.9741i000-1.6337-6.9741i000000000000Columns4through6000000000-3.47660000.0275+1.1755i0000.0275-1.1755i000Column70000005.0094V2=D2=Columns1through312.3669+6.6610i00012.3669-6.6610i0001.8504+4.6859i000000000000Columns4through60000000001.8504-4.6859i000-2.4345000-0.0000000Column7000000-0.00000.09840.07850.13300.04550.13240.04490.0946-0.39480.1692-0.22950.43190.0406-0.33550.1091-0.21600.2930Columns6through70.10960.21270.09820.18390.15860.33050.0066-0.01450.0005-0.01130.25670.56620.17700.4002>>A=[27577;74933;39838;59636;26854];B=[703795980137660;54795727112284;445523252894469;69588087619965;621173737299988];d1=eig(A)d2=eig(B)[V1,D1]=eig(A)[V2,D2]=eig(B)s1=svd(A)s2=svd(B)d1=27.86292.6062-2.2306+1.8926i-2.2306-1.8926i-5.0078d2=1.0e+003*Columns4through5000000-2.2306-1.8926i00-5.0078V2=Columns1through30.5371-0.52450.86000.32330.0727-0.28910.42330.7294-0.37880.4293-0.42880.00120.4934-0.0613-0.1826Columns4through5-0.1142+0.0175i-0.1142-0.0175i0.1919+0.3393i0.1919-0.3393i0.2293-0.2825i0.2293+0.2825i0.4703-0.2277i0.4703+0.2277i-0.6551-0.6551D2=1.0e+003*Columns1through3[4,1,0,1,2,0,1][-3,-4,-7,3,7,8,12][1,-10,7,-6,8,1,5]6>>A=[5765165;2310014;6420644;3963662;10760077;7244077;4867217];B=[3550123;3254625;1211346;3515212;4101201;-3-4-737812;1-107-6815];a1=sym(A)b1=sym(A)a2=rank(A)b2=rank(A)%b為奇異矩陣，所以不能求逆%a3=trace(A)b3=trace(A)a4=inv(A)a1=[5,7,6,5,1,6,5][2,3,1,0,0,1,4][6,4,2,0,6,4,4][3,9,6,3,6,6,2][10,7,6,0,0,7,7][7,2,4,4,0,7,7][4,8,6,7,2,1,7]b1=[5,7,6,5,1,6,5][2,3,1,0,0,1,4][6,4,2,0,6,4,4][3,9,6,3,6,6,2][10,7,6,0,0,7,7][7,2,4,4,0,7,7][4,8,6,7,2,1,7]a2=7b2=7a3=27b3=27a4=0.1249-0.1426+0.1548i-0.1426-0.1548i0.33570.2001+0.2826i0.2001-0.2826i0.40340.60210.60210.4462-0.1147+0.3618i-0.1147-0.3618i0.4137-0.3247+0.0374i-0.3247-0.0374i0.4064-0.0929-0.4405i-0.0929+0.4405iColumns4through60.3888-0.2542+0.3398i-0.2542-0.3398i-0.0123-0.3717-0.0529i-0.3717+0.0529i-0.78010.62410.62410.2996-0.2558+0.1058i-0.2558-0.1058i0.37040.1728+0.0865i0.1728-0.0865i-0.1053-0.1023-0.2178i-0.1023+0.2178i0.04690.2885-0.1854i0.2885+0.1854iColumn7-0.25430.4548-0.1545-0.2834-0.1041-0.57030.5368D1=Columns1through3Columns1through30.0306-0.3831i0.0306+0.3831i-0.1431-0.2811i-0.0138-0.4299i-0.0138+0.4299i-0.0257+0.0616i0.1817-0.2554i0.1817+0.2554i0.1802-0.1874i-0.1582-0.3661i-0.1582+0.3661i0.0014+0.0928i-0.0772-0.1682i-0.0772+0.1682i-0.1943+0.1869i0.48380.48380.82680.3697+0.0701i0.3697-0.0701i-0.2500+0.0316iColumns4through6-0.1431+0.2811i-0.02050.1932-0.0257-0.0616i0.1918-0.25700.1802+0.1874i-0.24320.12870.0014-0.0928i-0.12730.1923-0.1943-0.1869i-0.1376-0.12360.8268-0.68870.7843-0.2500-0.0316i0.6279-0.4612Column7-0.25750.4087-0.1835-0.34670.67740.0649-0.3869ans=1.0e-005*Columns1through5-0.3010-0.3411-0.2478-0.1875-0.1224-0.0903-0.1023-0.0743-0.0560-0.0368-0.2423-0.2750-0.1996-0.1510-0.0991-0.2910-0.3304-0.2400-0.1819-0.1189-0.3223-0.3654-0.2653-0.2004-0.1313-0.2990-0.3386-0.2459-0.1858-0.1215-0.2936-0.3325-0.2416-0.1830-0.1192Columns6through7-0.2610-0.3209-0.0783-0.0962-0.2102-0.2584-0.2525-0.3105-0.2796-0.3435-0.2592-0.3186-0.2545-0.3132ans=1.0e-007*Columns1through50.22260.13760.18380.19720.27830.23630.19270.18460.24800.27170.1799-0.03990.17890.05490.32950.17010.28460.10190.28152.6697-0.5809-0.13580.5730+0.3156i0.5730-0.3156iV1=Columns1through30.44520.57560.3203-0.3633i0.42560.1987-0.56970.4933-0.5629-0.0168+0.0797i0.46070.5010-0.0594+0.3503i0.4062-0.24780.5457-0.1010iColumns4through50.3203+0.3633i-0.0706-0.56970.6063-0.0168-0.0797i-0.4901-0.0594-0.3503i-0.52290.5457+0.1010i0.3373D1=Columns1through327.86290002.6062000-2.2306+1.8926i0000002.6697000-0.5809000-0.1358000000Columns4through50000000.5730+0.3156i000.5730-0.3156is1=28.78196.58304.23452.69491.4393s2=1.0e+003*2.80470.95650.75910.43550.1016實(shí)習(xí)總結(jié)：學(xué)生簽名：《MATLAB軟件實(shí)習(xí)》課程實(shí)習(xí)報(bào)告姓名專業(yè)數(shù)學(xué)與統(tǒng)計(jì)學(xué)院班級09級一班學(xué)號實(shí)習(xí)地點(diǎn)實(shí)習(xí)時(shí)間第十周實(shí)習(xí)章節(jié)4.3矩陣的基本變換與分解、4.4矩陣方程的計(jì)算機(jī)求解實(shí)習(xí)內(nèi)容：實(shí)習(xí)總結(jié)：學(xué)生簽名：《MATLAB軟件實(shí)習(xí)》課程實(shí)習(xí)報(bào)告姓名專業(yè)數(shù)學(xué)與統(tǒng)計(jì)學(xué)院班級學(xué)號實(shí)習(xí)地點(diǎn)實(shí)習(xí)時(shí)間第十一周實(shí)習(xí)章節(jié)7.1常系數(shù)線性微分方程的解析解方法、7.2微分方程問題的數(shù)值解法實(shí)習(xí)內(nèi)容：1、symst;u=exp(-2*t)*(sin(2*t)+pi/3)+cos(3*t);y=dsolve(['D5y+13*D4y+64*D3y+152*D2y+176*Dy+80*y=',char(u)],'y(0)=1','y(1)=3','y(pi)=2','Dy(0)=1','Dy(1)=2')2、symst;[x,y]=dsolve('D2x+5*Dx+4*x+3*y=exp(-6*t)*sin(4*t)','2*Dy+y+4*Dx+6*x=exp(-6*t)*cos(4*t)','x(0)=1','x(pi)=2','y(0)=0')3、x=dsolve('(1-t^2)*D2x-2*t*Dx+n*(n+1)*x=0')y=dsolve('t^2*D2x+t*Dx+(t^2-n^2)*x=0')x=C1*LegendreP(n,t)+C2*LegendreQ(n,t)y=C1*besselj(n,t)+C2*bessely(n,t)x=C1*LegendreP(n,t)+C2*LegendreQ(n,t)y=C1*besselj(n,t)+C2*bessely(n,t)實(shí)習(xí)總結(jié)：學(xué)生簽名：《MATLAB軟件實(shí)習(xí)》課程實(shí)習(xí)報(bào)告姓名專業(yè)數(shù)學(xué)與統(tǒng)計(jì)學(xué)院班級學(xué)號實(shí)習(xí)地點(diǎn)實(shí)習(xí)時(shí)間第十二周實(shí)習(xí)章節(jié)8.1插值與數(shù)據(jù)擬合、8.3由已知數(shù)據(jù)擬合數(shù)學(xué)模型實(shí)習(xí)內(nèi)容：實(shí)習(xí)內(nèi)容：18x0=0:.1:1;y0=(x0.^2-3*x0+5).*exp(-5*x0).*sin(x0);p3=polyfit(x0,y0,3);x=0:.01:1;ya=(x.^2-3*x+5).*exp(-5*x).*sin(x);y1=polyval(p3,x);plot(x,y1,x,ya,x0,y0,'o')p4=polyfit(x0,y0,4);y2=polyval(p4,x);p5=polyfit(x0,y0,5);y3=polyval(p5,x);p8=polyfit(x0,y0,8);y4=polyval(p8,x);xplot(x,ya,x0,y0,'o',x,y2,x,y3,x,y4)19x0=-1+2*[0:10]/10;y0=1./(1+25*x0.^2);x=-1:.01:1;ya=1./(1+25*x.^2);p3=polyfit(x0,y0,3);y1=polyval(p3,x);p5=polyfit(x0,y0,5);y2=polyval(p5,x);p8=polyfit(x0,y0,8);y3=polyval(p8,x);p10=polyfit(x0,y0,10);y4=polyval(p10,x);plot(x,ya,x,y1,x,y2,'-',x,y3,'--',x,y4,':')>>maple('with(numtheory):');f=maple(['cfe:=cfrac(pi,20)'])n=maple('nthnumer','cfe',4);d=maple('nthdenom','cfe',4);vpa(n)/vpa(d)Warning:Warning,theprotectednameorderhasbeenredefinedandunprotectedf=cfe:=3+1/(7+1/(15+1/(1+1/(292+1/(1+1/(1+1/(1+1/(2+1/(1+1/(3+1/(1+1/(14+1/(2+1/(1+1/(1+1/(2+1/(2+1/(2+1/(2+1/(1+`...`))))))))))))))))))))ans=3.>>symsx;fun='sin(x)*exp(-x)/(x+1)^3';maple('with(numtheory):');f=maple(['cfe:=cfrac('fun',x,10)'])n=collect(maple('nthnumer','cfe',8),x);d=collect(maple('nthdenom','cfe',8),x);[n,d]=numden(n/d);G=n/d;n=collect(maple('nthnumer','cfe',10),x);d=collect(maple('nthdenom','cfe',10),x);[n,d]=numden(n/d);G1=n/d;ezplot(fun,[0,2])holdonezplot(G,[0,2]);ezplot(G1,[0,2])figure;ezplot(fun,[0,5])holdonezplot(G,[0,5]);ezplot(G1,[0,5])f=cfe:=x/(1+4*x/(1-5*x/(3+43*x/(20-337*x/(43+28274*x/(1685-*x/(-*x/(+8444*x/(-*x/(+`...`)))