matlab基礎(chǔ)與應(yīng)用部分習(xí)題答案

matlab基礎(chǔ)與應(yīng)用部分習(xí)題答案matlab基礎(chǔ)與應(yīng)用部分習(xí)題答案matlab基礎(chǔ)與應(yīng)用部分習(xí)題答案xxx公司matlab基礎(chǔ)與應(yīng)用部分習(xí)題答案文件編號：文件日期：修訂次數(shù)：第1.0次更改批準(zhǔn)審核制定方案設(shè)計(jì)，管理制度作業(yè)一4、寫出完成下列操作的命令。（1）將矩陣A第2~5行中第1,3,5列元素賦給矩陣B。>>M=[0:1:48];>>A=reshape(M,7,7)A=0714212835421815222936432916233037443101724313845411182532394651219263340476132027344148>>B=A(2:5,1:2:5)B= 11529216303173141832（2）刪除矩陣A的第七行元素。>>A(7,:)=[]A=071421283542181522293643291623303744310172431384541118253239465121926334047(3)將矩陣A的每個(gè)元素值加30。>>A=A+30A=303744515865723138455259667332394653606774334047546168753441485562697635424956637077(4求矩陣A的大小和維素。sizeA=size(A)dA=ndims(A)sizeA=67dA=2（5）將向量t的0元素用機(jī)器零來代替。>>t=[123405];t=123405>>find(t==0)ans=5>>t(5)=epst=（6）將含有12個(gè)元素的向量x轉(zhuǎn)換成3*4矩陣。>>x=[0:11]x=01234567891011>>y=reshape(x,3,4)y=03691471025811（7）求一個(gè)字符串的ASCII。>>abs('a')ans=97（8）求一個(gè)ASCII所對應(yīng)的字符。>>char(99)ans=c5、已知A[97673410-7875655325-23-590=-12547]寫出完成下列操作的命令。取出A的前3列構(gòu)成矩陣B，前兩列構(gòu)成矩陣C，右下角3*2子矩陣構(gòu)成矩陣D，B與C的乘積構(gòu)成矩陣E。>>B=A(1:3,:)B=97673410-7875655325-23-59>>C=A(:,1:2)C=9767-78753250-12>>D=A(2:4,3:4)D=655-23-59547>>E=B*CE=527111574-1133666419783112輸出【50,100】范圍的全部元素。>>find(A>=50&A<=100)ans=1561012>>output=[A(1)A(5)A(6)A(10)A(12)]output=9767756554二、實(shí)驗(yàn)題>>x=-74*pi/180;y=-27*pi/180;sin((x*x+y*y))/(sqrt(tan((abs(x+y))))+pi)ans=->>a=::,exp.*a).*sin(a+>>x=[24;5];*log(x+sqrt(1+x.*x))ans=>>A=[3542;34-457;879015],B=[1-267;2874;930]>>A*BA.*Bans=12943241977-407-105240259112489ans=3-10813468-3605187832700>>A^3A.^3>>A/BB\A>>[A,B][A([1,3],:);B^2]>>a=1+2i,b=3+4i,c=exp(pi/6*i)c+a*b/(a+b)作業(yè)二從鍵盤輸入一個(gè)4位整數(shù)，按如下規(guī)則加密后輸出。加密規(guī)則：每位數(shù)字都加上7，然后用和除以10的余數(shù)取代該數(shù)字；然后將第1位數(shù)和第3位數(shù)交換，第2位數(shù)和第4位數(shù)交換。M文件名為：程序：x=input('輸入一個(gè)四位整數(shù)：');y=[fix(x/1000),mod(fix(x/100),10),mod(fix(x/10),10),mod(x,10)]z=mod((y+7),10)x=z(3)*1000+z(4)*100+z(1)*10+z(2)結(jié)果：輸入一個(gè)四位整數(shù)：1234y=1234z=8901x=189硅谷公司員工的工資計(jì)算方法如下。工作時(shí)數(shù)超過120h者，超過部分加發(fā)15%。工作時(shí)數(shù)低于60h者，扣發(fā)700元。其余按每小時(shí)84元計(jì)發(fā)。試編程按輸入的工號和該號員工的工時(shí)數(shù)，計(jì)算應(yīng)發(fā)工資。程序：gonghao=input('輸入工號');hour=input('輸入工時(shí)');every=84;ifhour>120gongzi=every*120+*every*(hour-120);elseifhour<60gongzi=every*hour-700;elsegongzi=devery*hour;endformatbank;display([gonghao,gongzi])（1）結(jié)果：輸入工號12輸入工時(shí)130ans=（2）結(jié)果：輸入工號12輸入工時(shí)56ans=（3）結(jié)果：輸入工號12輸入工時(shí)80ans=根據(jù)，求π的近似值。當(dāng)n分別取100、1000、10000時(shí)，結(jié)果是多少要求：分別用循環(huán)結(jié)構(gòu)和向量運(yùn)算（使用sun函數(shù)）來實(shí)現(xiàn)。（1）程序：（循環(huán)結(jié)構(gòu)）：n=input('inputn:');s=0;fori=1:ns=s+1/i^2;pai=sqrt(s*6);enddisplay(pai)inputn:100pai=（2）程序：（向量結(jié)構(gòu)）：n=input('inputn:');i=1:n;pai=sqrt(sum(1./i.^2)*6);display(pai);結(jié)果：inputn:100pai=根據(jù)，求：y﹤3時(shí)的最大n值；與（1）的n值對應(yīng)的y值。y=0;k=0;whiley<3k=k+1;y=y+1/(2*k-1);enddisplay([k-1,y-1/(2*k-1)])結(jié)果：ans=考慮一下迭代公式：其中，a、b為正的常數(shù)。編寫程序求迭代的結(jié)果，迭代的終止條件為||，迭代初值，迭代次數(shù)不超過500次。程序：x0=;a=input('a=');b=input('b=');x=a/(b+x0);m=1;whileabs(x-x0)>1e-5&m<500x0=x;x=a/(b+x0);m=m+1;enddisplay([m,x]);display([(-b+sqrt(b^2+4*a))/2,(-b-sqrt(b^2+4*a))/2]);結(jié)果：a=5b=5ans=ans=如果迭代過程收斂于r，那么r的準(zhǔn)確值是，當(dāng)（a，b）的值?。?，1）、（8,3）、（10,）時(shí)，分別對迭代結(jié)果和準(zhǔn)確值進(jìn)行比較。a=1b=1ans=ans=a=8b=3ans=ans=a=10b=ans=ans=6、已知當(dāng)時(shí)，y的值是多少程序：該題使用了function的調(diào)用。函數(shù)文件functionb=f(n);b=n+10*log(n^2+5);結(jié)果：y=f(40)/(f(30)+f(20))y=當(dāng)時(shí)，y的值是多少程序：函數(shù)文件functionf=fun(n)f=sum([1:n].*([1:n]+1));結(jié)果：y=fun(40)/(fun(30)+fun(20))y=作業(yè)三繪制下列曲線。（1）（2）（3）（4）解：（1）程序：x=-5::5;y=x-x.^3/(3*2*1);plot(x,y);(2)程序：x=-10::10;y=exp(-x.*x/2)/(2*pi);>>plot(x,y);程序：f='x.^2+2*y.^2-64';ezplot(f,[-10,10,-10,10]);ezplot('x.^2+2*y.^2-64',[-10,10,-10,10]);（4）程序：t=0::8*pi;>>x=t.*sin(t);>>y=t.*cos(t);>>plot(x,y); 已知，，完成下列操作：在同一坐標(biāo)系下用不同的顏色和線型繪制兩條曲線，給曲線添加文字說明；解：程序：x1=linspace(0,1,100);>>y1=2*;>>t=linspace(0,pi,100);>>x=sin(3*t).*cos(t);>>y=sin(3*t).*sin(t);>>plot(x1,y1,'k:',x,y,'bp');>>text,1,'y=');>>text,,'x=sin(3t)cos(t)');>>text,,'x=sin(3t)sin(t)'); 以子圖形式，分別用條形圖和散點(diǎn)圖繪制兩條曲線，并為各子圖添加函數(shù)標(biāo)題。程序：subplot(1,2,1);bar(x1,'group');title('y1=');subplot(1,2,2);barh(t,'stack');title('stack');subplot(1,2,1);>>scatter(x1,y1,10);>>title('y=');>>subplot(1,2,2);>>scatter(x,y,10);>>title('x=sin(3t)cos(t),y=sin(3t)cos(t)');分別用plot和fplot函數(shù)繪制函數(shù)的曲線，分別分析兩條曲線的差別。解：程序：subplot(1,2,1);x=1:1:100;y=sin(1./x);plot(x,y);subplot(1,2,2);>>fplot('sin(1/x)',[1,100]);設(shè)，，在同一圖形窗口采用子圖的形式繪制條形圖、階梯圖、桿圖和對數(shù)坐標(biāo)圖。解：程序：t=-pi::pi;>>y=1./(1+exp(-t));>>bar(t,y);>>subplot(2,2,2);>>stairs(t,y);>>subplot(2,2,3);>>stem(t,y);>>subplot(2,2,4);>>semilogy(t,y);結(jié)果：繪制下列極坐標(biāo)圖。（1）（2），解：（1）程序：theta=linspace(-pi,pi,100);ro=5.*cos(theta)+4;polar(theta,ro);（2）程序：fi=linspace(0,2*pi,100);a=1r=a.*(1+cos(fi));polar(fi,r);6、繪制下列三維圖線。（1）（2）解：（1）程序：t=0:pi/20:2*pi;x=exp(-t./20).*cos(t);y=exp(-t./20).*sin(t);z=t;plot3(x,y,z);結(jié)果：程序：t=0::1;x=t;y=t.^2;z=t.^3;plot3(x,y,z);結(jié)果：7、已知，繪制其在范圍內(nèi)的曲面圖和等高線。解：程序：x=-30::0;y=0::30;[x,y]=meshgrid(x,y);z=10.*sin(sqrt(x.^2+y.^2))./sqrt(1+x.^2+y.^2);meshc(x,y,z);結(jié)果：8、已知，繪制其曲面圖，并將部分鏤空。解：程序：x=linspace(-3,3,100);y=linspace(-3,3,100);[xy]=meshgrid(x,y);fxy=-5./(1+x.^2+y.^2);i=find(abs(x)<=&abs(y)<=;fxy(i)=NaN;surf(x,y,fxy)9、繪制曲面圖形，應(yīng)用插值著色處理，并設(shè)置光照效果。解：程序：u=linspace(1,,100);v=linspace(-pi,pi,1000);[u,v]=meshgrid(u,v);x=3.*u.*sin(v);y=2.*u.*cos(v);z=4*u.^2;x=3*u.*sin(v);y=2*u.*cos(v);z=4*u.^2;surf(x,y,z);shadinginterp;light('position',[1，1，1]);10、設(shè)計(jì)一個(gè)藍(lán)色球體沿正弦曲線運(yùn)動(dòng)的動(dòng)畫。解：程序：t=0:pi/250:10*pi;y=sin(t);plot(t,y,'m');h=line('color',[0,0,1],'marker','.','markersize',50);fori=1:length(t);set(h,'xdata',t(i),'ydata',y(i));M(i)=getframe;endmovie(M);第四次作業(yè)P103思考題5、寫出完成下列操作的命令。（1）建立3階單位矩陣A。A=eye(3)A=100010001（2）建立5*6隨機(jī)矩陣A，其元素為[100，200]范圍內(nèi)的隨機(jī)整數(shù)。a=100;b=200;A=fix(a+(b-a)*rand(5,6))A=195176161140105120123145179193135119160101192191181160148182173141100127189144117189113119（3）產(chǎn)生均值為1、方差為的500個(gè)正態(tài)分布的隨機(jī)數(shù)。y=1+sqrt*randn(25,20)（4）產(chǎn)生和A同樣大小的幺矩陣。B=ones(size(A))（5）將矩陣A的對角線的元素加30。A=[1,2,3,4,5,6;5,6,2,1,3,5;7,8,9,4,5,6]A=123456562135789456>>A=A+30*eye(size(A))A=312345653621357839456（6）從矩陣A提取主對角線元素，并以這些元素構(gòu)成對角陣B。C=diag(diag(A))C=310003600039實(shí)驗(yàn)題：4、分別用矩陣求逆、左除運(yùn)算符、矩陣分解等方法求下列對角方程組，并分析程序的執(zhí)行時(shí)間。其中A=2-1000-12-1000-12-1000-12-1000-12>>b=[1;0;0;0;0]b=10000解：程序：%矩陣求逆ticx1=inv(A)*btoc%左除運(yùn)算ticx2=A\btoc%矩陣分解%Q、R分解tic[Q,R]=qr(A);x4=R\(Q\b)；toc%L、U分解tic[L,U]=lu(A);x3=U\(L\b)toc%Cholesky分解ticR=chol(A);x5=R\(R'\b)toc結(jié)果：P121第五章實(shí)驗(yàn)題利用，MATLAB提供的rand函數(shù)生成30000個(gè)符合均勻的隨機(jī)數(shù)，然后檢驗(yàn)隨機(jī)數(shù)的性質(zhì)。均值和標(biāo)準(zhǔn)方差。解：A=rand(1,30000);>>mean(A)std(A)最大元素和最小元素。max(A)min(A)大于的隨機(jī)數(shù)個(gè)數(shù)占總數(shù)的百分比。size(find(A>)/size(A)ans=2、在某處測的海洋不同深部處水溫度數(shù)據(jù)如表。用插值法求出水深500m，900m和1500m處的水溫（℃）。水深（m）4667159504221635水溫（℃）解：程序：m=[466,715,950,1422,1635];w=[,,,,];h=[500,900,1500];wi=interp1(沒,w,h,'spline')3、用5次多項(xiàng)式p（x）在區(qū)間【1,10】區(qū)間內(nèi)逼近函數(shù)lgx，并繪制出lgx和p（x）在【1,101】區(qū)間的函數(shù)曲線。解：x=1::101;y=log(x);f=polyfit(x,y,5);yi=polyval(f,x);plot(x,y,'bp',x,yi,'k-p')4、已知，取N=64，對t從0~5s采樣，用fft作快速傅立葉變換，并繪制相應(yīng)的振幅--頻率圖。解：N=64;T=5;t=linspace(0,T,N);x=exp(-t);dt=t(2)-t(1);f=1/dt;X=fft(x);F=X(1:N/2+1);f=f*(0:N/2)/N;plot(f,abs(F),'r-+')xlabel('Frequency')ylabel('|F(k)|')5、有3個(gè)多項(xiàng)式，，，試進(jìn)行下列操作。（1）求。解、程序：p1=[1,2,0,7];p2=[0,0,1,-2];p3=[1,0,5,1];p12=conv(p1,p2);p=p12+[zeros(1,size(p12,2)-size(p3,2)),p3]結(jié)果：p=0011-412-13（2）求p（x）的根。程序：roots(p)結(jié)果：ans=+-（3）設(shè)A矩陣，當(dāng)以矩陣A的每一元素為自變量時(shí)，求的值，當(dāng)以矩陣A為自變量時(shí)，求的根。程序：A=[-1,4,3;2,1,5;0,5,6];px=polyval(p,A)pxm=polyvalm(p,A)結(jié)果：px=-292919519-3697-136971427pxm=3912084327350226934207720377558926、求非線性方程的數(shù)值解。（1）、求在附近的根。解：程序：z=fzero('3*x-sin(x)+1',0)（2）、，初值。建立函數(shù)文件：functionF=fun(X)x=X(1);y=X(2);F(1)=x*x+y*y-9;F(2)=x+y-1;在命令窗口輸入：x=fsolve(@fun,[3,0]',optimset('Display','off'))第7章符號與計(jì)算課后實(shí)驗(yàn)題1、已知A=求矩陣的秩、逆矩陣和行列式的值。（這里的A是符號矩陣）解：程序及結(jié)果：A=sym('[1,2,3;x,y,z;3,2,1]');>>rank(A)%矩陣秩的值ans=3>>inv(A)%逆矩陣的值ans=[1/4*(y-2*z)/(-2*y+z+x),1/(-2*y+z+x),-1/4*(-2*z+3*y)/(-2*y+z+x)][-1/4*(x-3*z)/(-2*y+z+x),-2/(-2*y+z+x),1/4*(-z+3*x)/(-2*y+z+x)][1/4*(2*x-3*y)/(-2*y+z+x),1/(-2*y+z+x),-1/4*(-y+2*x)/(-2*y+z+x)]>>det(A)%矩陣的行列值ans=-8*y+4*z+4*x2、求函數(shù)的符號導(dǎo)數(shù)。（1）求。解：程序及結(jié)果：symsx;>>y=sqrt(‘x+sqrt(x+sqrt(x’)));>>diff(y)%一階導(dǎo)數(shù)ans=1/2/(x+(x+x^(1/2))^(1/2))^(1/2)*(1+1/2/(x+x^(1/2))^(1/2)*(1+1/2/x^(1/2)))>>diff(y,x,2)%二階導(dǎo)數(shù)ans=-1/4/(x+(x+x^(1/2))^(1/2))^(3/2)*(1+1/2/(x+x^(1/2))^(1/2)*(1+1/2/x^(1/2)))^2+1/2/(x+(x+x^(1/2))^(1/2))^(1/2)*(-1/4/(x+x^(1/2))^(3/2)*(1+1/2/x^(1/2))^2-1/8/(x+x^(1/2))^(1/2)/x^(3/2))（2）已知，求。解：程序：symsxy;>>fxy=sin(x^2*y)*exp(-x^2-y);>>c=diff(diff(fxy,x),y)結(jié)果：c=-2*sin(x^2*y)*x^3*y*exp(-x^2-y)+2*cos(x^2*y)*x*exp(-x^2-y)-2*cos(x^2*y)*x*y*exp(-x^2-y)-2*cos(x^2*y)*x^3*exp(-x^2-y)+2*sin(x^2*y)*x*exp(-x^2-y)3、求積分。（1）（2）（1）解：程序：f=sym('1/(x^4+1)');>>s1=int(f)（2）解：程序：symsxt;>>f=sym('(-2*x^2+1)/(2*x^2-3*x+1)^2');>>int(f,cos(t),exp(2*t))結(jié)果：ans=-(-2*exp(2*t)*cos(t)^2-exp(2*t)+cos(t)+2*exp(2*t)^2*cos(t))/(2*exp(2*t)-1)/(exp(2*t)-1)/(2*cos(t)-1)/(cos(t)-1)4、求級數(shù)的和函數(shù)，并求之和。解：程序：symsxn;>>s1=symsum(2/((2*n+1)*(2*x+1)^(2*n+1)),n,0,inf)結(jié)果：s1=1/(2*x+1)*(4*x^2+4*x+1)^(1/2)*log((1+1/(4*x^2+4*x+1)^(1/2))/(1-1/(4*x^2+4*x+1)^(1/2)))s2=symsum(2/((2*n+1)*(2*x+1)^(2*n+1)),n,1,5)s2=2/3/(2*x+1)^3+2/5/(2*x+1)^5+2/7/(2*x+1)^7+2/9/(2*x+1)^9+2/11/(2*x+1)^11>>s=s1+s2s=1/(2*x+1)*(4*x^2+4*x+1)^(1/2)*log((1+1/(4*x^2+4*x+1)^(1/2))/(1-1/(4*x^2+4*x+1)^(1/2)))+2/3/(2*x+1)^3+2/5/(2*x+1)^5+2/7/(2*x+1)^7+2/9/(2*x+1)^9+2/11/(2*x+1)^115、求函數(shù)在的泰勒展開式。(1)(2)（1）解：程序：sym('x');>>y=(exp(x)+exp(-x))/2;>>taylor(y,x,5,0)結(jié)果：ans=1+1/2*x^2+1/24*x^4（2）解：程序：>>symsax>>s=taylor(exp(-5*x)*sin(3*x+pi/3),5,a)結(jié)果：s=exp(-5*a)*sin(3*a+1/3*pi)+(3*exp(-5*a)*cos(3*a+1/3*pi)-5*exp(-5*a)*sin(3*a+1/3*pi))*(x-a)+(8*exp(-5*a)*sin(3*a+1/3*pi)-15*exp(-5*a)*cos(3*a+1/3*pi))*(x-a)^2+(33*exp(-5*a)*cos(3*a+1/3*pi)+5/3*exp(-5*a)*sin(3*a+1/3*pi))*(x-a)^3+(-161/6*exp(-5*a)*sin(3*a+1/3*pi)-40*exp(-5*a)*cos(3*a+1/3*pi))*(x-a)^46、求非線性方程的符號解。（1）（2）（1）解：程序：s=solve(sym('x^3+a*x+1=0'))（2）解：程序：三種方法：[xy]=solve(sym('sqrt(x^2+y^2)-100=0'),sym('3*x+5*y-8=0'))結(jié)果：x=12/17-10/17*21246^(1/2)12/17+10/17*21246^(1/2)y=20/17+6/17*21246^(1/2)20/17-6/17*21246^(1/2)>>%或者下面的方式也可以>>s1=sym('sqrt(x^2+y^2)-100=0');>>s2=sym('3*x+5