matlab習(xí)題(課程設(shè)計(jì))

可編輯版6/61.已知矩陣,矩陣,解線性方程X,使得A.X=B,并將矩陣A的右下角2×3子矩陣賦給矩陣C。程序：A=sym<[5,7,6,5,1;7,10,8,7,2;6,8,10,9,3;5,7,9,10,4;1,2,3,4,5]>;B=sym<[24,96;34,136;36,144;35,140;15,60]>;X=A\BC=A<[4,5],[3,4,5]>運(yùn)行結(jié)果：X=[1,4][1,4][1,4][1,4][1,4]C=[9,10,4][3,4,5]2.根據(jù)下面兩個(gè)矩陣執(zhí)行下列的矩陣運(yùn)算命令,并回答有關(guān)的問(wèn)題A+5*B和A-B+I分別是多少〔其中I為單位矩陣？A.*B和A*B將分別給出什么結(jié)果,它們是否相同,為什么？得出A.^B、A/B及A\B的結(jié)果,并分別解釋它們的物理意義。程序：A=[4,12,20;12,45,78;20,78,136];B=[1,2,3;4,5,6;7,8,9];I=[1,0,0;0,1,0;0,0,1];C=A+5*BD=A-B+IE=A.*BF=A*BG=A.^BH=A/BJ=A\BA=sym<[41220;124578;2078136]>；B=sym<[123;456;789]>；I=sym<[100;010;001]>；A+5*BA-B+IA.*B%A.為A的共軛矩陣,本式為A的共軛矩陣和B相乘A*B%A直接和矩陣B相乘,結(jié)果與上面顯然不同A.^B%A的共軛矩陣的B次方A/B%右除解XB=AA\B%左除解AX=B運(yùn)行結(jié)果：ans=[9,22,35][32,70,108][55,118,181]ans=[4,10,17][8,41,72][13,70,128]ans=[4,24,60][48,225,468][140,624,1224]ans=[192,228,264][738,873,1008][1284,1518,1752]ans=[4,144,8000][20736,184528125,225199600704]Warning:Systemisrankdeficient.Solutionisnotunique.>InC:\MATLAB6p5\toolbox\symbolic\@sym\mldivide.matline38InC:\MATLAB6p5\toolbox\symbolic\@sym\mrdivide.matline24InC:\MATLAB6p5\work\Untitled2.matline9ans=[28/3,-4/3,0][40,-7,0][212/3,-38/3,0]Warning:Systemisrankdeficient.Solutionisnotunique.>InC:\MATLAB6p5\toolbox\symbolic\@sym\mldivide.matline38InC:\MATLAB6p5\work\Untitled2.matline10ans=[-1/12,5/6,7/4][1/9,-1/9,-1/3][0,0,0]分別用for和while循環(huán)結(jié)構(gòu)編寫(xiě)程序,求出并考慮一種避免循環(huán)的簡(jiǎn)潔方法來(lái)求和,并比較各種算法的運(yùn)行時(shí)間?！沧ⅲ阂还彩侨N方法for的方法：程序：A=63;n=0;k=0;fori=0:An=2^i;k=k+n;endk運(yùn)行結(jié)果：k=1.8447e+019While的方法：程序：i=0;n=0;k=0;whilei<64n=2^i;i=i+1;k=k+n;endk運(yùn)行結(jié)果：k=1.8447e+019k=0;k=1*<2^64-1>/<2-1>;kk=1.8447e+0194.畫(huà)圖題要求：1、兩根曲線畫(huà)在一幅圖中,不采用子圖的方式。2、橫坐標(biāo)〔0~5π,縱坐標(biāo)為〔-1.5~+1.53、要求給橫坐標(biāo)軸加說(shuō)明"t<deg>",縱坐標(biāo)軸加"magnitude",并且加網(wǎng)格4、在〔π/2,1.2處加創(chuàng)建說(shuō)明性的文字"這是我的程序運(yùn)行結(jié)果！"5、圖中,sint用紅色*線繪制,cost用籃色實(shí)線繪制。6、給圖形加標(biāo)題"sinewavefromzeroto5/pi"程序：clf;t=0:pi/50:5*pi;y1=sin<t>;y2=cos<t>;plot<t,y1,'r*',t,y2,'b'>,gridon,axis<[0,5*pi,-1.5,1.5]>text<pi/2,1.2,'\fontsize{16}\fontname{隸書(shū)}這是我的程序運(yùn)行結(jié)果！'>title<'sinewavefromzeroto5/pi'>xlabel<'t<deg>'>ylabel<'magnitude'>運(yùn)行結(jié)果：5.體會(huì)各種繪圖命令及效果。假設(shè)用戶有下面的繪圖數(shù)據(jù)向量t=0:0.1:10;y=sin<tan<t>>-tan<sin<t>>;且繪圖函數(shù)采用下面各個(gè)函數(shù),如polar<>,bar<>,stem<>,stairs<>,那么試得出并解釋所得出的結(jié)果。程序：clf;t=0:.1:10;y=sin<tan<t>>-tan<sin<t>>;subplot<2,2,1>,polar<t,y,'b-.'>;subplot<2,2,2>,bar<t,y,'y--'>;subplot<2,2,3>,stem<t,y,'g'>;subplot<2,2,4>,stairs<t,y,':r'>;運(yùn)行結(jié)果：6.證明：函數(shù)z=xy的圖形是雙曲拋物面?！蔡崾荆涸趨^(qū)域-2≤x≤2,-2≤y≤2上作出它的圖形。clfx=-2:2/10:2;y=x; [X,Y]=meshgrid<x,y>; Z=X.*Y; surf<X,Y,Z>; colormap<hot> xlabel<'x'>,ylabel<'y'>,zlabel<'z'> 7.求解微分方程y=dsolve<'Dy-2*y/<x+1>-<x+1>^2.5=0','x'>y=2/3*<x+1>^<3/2>*x^2+4/3*<x+1>^<3/2>*x+2/3*<x+1>^<3/2>+C1*x^2+2*C1*x+C18.用subplot分別在不同的坐標(biāo)系下作出下列四條曲線,為每幅圖形加上標(biāo)題,①概率曲線〔-5≤x≤5；②四葉玫瑰線〔0≤x≤2pi；③葉形線〔-2≤x≤2；④曳物線〔-1≤x≤1。clfx=-5:.1:5;y=exp<-<x.^2>>;subplot<2,2,1>,plot<x,y,'r'>;title<'概率曲線'>xlabel<'x'>ylabel<'y'>gridonphi=0:.1:2*pi;r=sin<2*phi>;subplot<2,2,2>,polar<phi,r,'r'>;title<'四葉玫瑰線'>xlabel<'phi'>ylabel<'r'>gridonsymstx=-2:0.1:2;x=3*t/<1+t.^3>;y=3*<t.^2>/<1+t.^3>;subplot<2,2,3>,ezplot<'3*t/<1+t.^3>','3*<t.^2>/<1+t.^3>',[020]>title<'葉形線'>xlabel<'x'>ylabel<'y'>gridonsymsx,y;subplot<2,2,4>,f1=inline<'log<<1+sqrt<1-y.^2>>./y>-sqrt<1-y.^2>-x'>;f2=inline<'log<<1-sqrt<1-y.^2>>./y>+sqrt<1-y.^2>-x'>;ezplot<f1>;holdon;ezplot<f2>;title<'曳物線'>gridon;求,并簡(jiǎn)化結(jié)果。symsabcdD=[a^2+<1/a^2>,a,1/a,1;b^2+<1/b^2>,b,1/b,1;c^2+<1/c^2>,c,1/c,1;d^2+<1/d^2>,d,1/d,1];DA=det<D>DB=simple<simple<DA>>DA=<c*d^2*b^3*a^4-c^2*d*b^3*a^4-c^3*d^2*b*a^4+c^3*d*b^2*a^4+c^2*d^3*b*a^4-c*d^3*b^2*a^4-c*d^2*a^3*b^4+c^2*d*a^3*b^4+c^3*d^2*a*b^4-c^3*d*a^2*b^4-c^2*d^3*a*b^4+c*d^3*a^2*b^4+b*d^2*a^3*c^4-b^2*d*a^3*c^4-b^3*d^2*a*c^4+b^3*d*a^2*c^4+b^2*d^3*a*c^4-b*d^3*a^2*c^4-b*c^2*a^3*d^4+b^2*c*a^3*d^4+b^3*c^2*a*d^4-b^3*c*a^2*d^4-b^2*c^3*a*d^4+b*c^3*a^2*d^4+c*d^2*b^3-c^2*d*b^3-c^3*d^2*b+c^3*d*b^2+c^2*d^3*b-c*d^3*b^2-c*d^2*a^3+c^2*d*a^3+c^3*d^2*a-c^3*d*a^2-c^2*d^3*a+c*d^3*a^2+b*d^2*a^3-b^2*d*a^3-b^3*d^2*a+b^3*d*a^2+b^2*d^3*a-b*d^3*a^2-b*c^2*a^3+b^2*c*a^3+b^3*c^2*a-b^3*c*a^2-b^2*c^3*a+b*c^3*a^2>/a^2/c^2/d^2/b^2DB=<d-c>*<-c+b>*<-d+b>*<-c+a>*<-d+