Matlab復(fù)習(xí)題

1、84ye4、編程實現(xiàn)當(dāng)t-,間隔為1°時求解正弦和余弦值。解：t=(-1*pi:1/180:pi);y1=sin(t)y2=cos(t)5、利用rand函數(shù)產(chǎn)生（0,1）間的均勻分布的10*10隨機矩陣A，然后統(tǒng)計A中大于等于0.6的元素的個數(shù)。解：A=rand(10);B=A >= 0.6;C=sum(B);count=sum(C)運行結(jié)果（每次運行結(jié)果是不同的，僅作參考）：count=326、利用randn函數(shù)產(chǎn)生均值為0，方差為1的10*10隨機矩陣A，然后統(tǒng)計A中大于-0.5且小于0.5的元素的個數(shù)。解：A=randn(10);B=(A<0.5)&(A&g

2、t;-0.5);C=sum(sum(B)運行結(jié)果（每次運行結(jié)果是不同的，僅作參考）：C=488、有一矩陣A，找出矩陣中值等于1的元素，并將它們重新排列成列向量B。解：A=2*rand(4);k=find(A<=1);A(k)=;%刪除下標(biāo)為k的元素B=A'9、在一測量矩陣A（100*3）中，存在有奇異值（假設(shè)大于100的置認(rèn)為是奇異值），編程實現(xiàn)刪去奇異值所在的行。解：A=120*randn(10,3);i,j=find(A>100);A(i,:)= %刪去存在奇異值的行運行結(jié)果（每次運行結(jié)果是不同的，僅作參考）：A =49.5355 -23.7550 -73.0269-1

3、18.4354 39.3214 -88.447291.1482 -28.5962 -209.985521.2336 -74.0239 -9.5871-15.8184 72.1322 22.044471.4429 11.0770 34.89484、在給定的100*100矩陣中，刪去整行為0的行，刪去整列為0的列。解：A=diag(1 2 3 4,1)B=any(A)i,j=find(B=0)A(:,i)= %刪除全為0的列B=any(A')i,j=find(B=0)A(j,:)= %刪除全為0的行運行結(jié)果：初始值：A =0 1 0 0 00 0 2 0 00 0 0 3 00 0 0 0

4、 40 0 0 0 0操作后：A =1 0 0 00 2 0 00 0 3 00 0 0 4164ye2、將圖形窗口分成兩個，分別繪制正割和余割曲線，并加上標(biāo)注。程序為：x1=0:pi10:2*pi;figure(1)subplot(2,1,1)plot(x,sec(x),'k-'),grid onlegend('ity=sec(x)')title('y=sec(x)')xlabel('x'), ylabel('y')subplot(2,1,2)plot(x,csc(x),'k-'),grid on

5、legend('ity=csc(x)')title('y=csc(x)')xlabel('x'), ylabel('y')3、設(shè)有函數(shù)y=exp(x+5)+x.3，在半對數(shù)坐標(biāo)系中繪制曲線。程序為：x=1:0.01:10;y=exp(x+5)+x.3;figure(1)subplot(2,1,1)plot(x,y,'r-'),grid onlegend('ity=exp(x+5)+x.3')title('平面坐標(biāo)')xlabel('x'), ylabel('y

6、')subplot(2,1,2)semilogx(x,y,'k-'),grid on %半對數(shù)坐標(biāo)軸legend('ity=exp(x+5)+x.3')title('半對數(shù)坐標(biāo)')xlabel('x'), ylabel('y')5、將圖形窗口分成兩個，分別繪制以下函數(shù)在-3,3區(qū)間上的曲線，并利用axis 調(diào)整軸刻度，使他們具有相同縮放尺度。y1=2x+5； y2=x2-3x+1。程序為：x=-3:0.1:3;y1=2*x+5;y2=x.2-3*x+1;figure(1)subplot(2,2,1)plot

7、(x,y1,'r-'),grid onlegend('ity1=2*x+5')title('y1=2x+5')xlabel('x'), ylabel('y1')subplot(2,2,2)plot(x,y2,'k-'),grid onlegend('ity2=x.2-3*x+1')title('y2=x2-3x+1')xlabel('x'), ylabel('y2')subplot(2,2,3)plot(x,y1,'r-'

8、;),grid onlegend('ity1=2*x+5')title('調(diào)整后的y1=2x+5')axis(-3 3 -10 10)xlabel('x'), ylabel('y1')subplot(2,2,4)plot(x,y2,'k-'),grid onlegend('ity2=x.2-3*x+1')title('調(diào)整后的y2=x2-3x+1')axis(-3 3 -10 10) %調(diào)整坐標(biāo)軸xlabel('x'), ylabel('y2')7、繪

9、制餅圖。程序為：x=190 33 45 42 45;explode=0 1 0 0 0;figure(1)subplot(2,1,1)colormap hsvpie(x,explode)gtext('生活費')gtext('資料費')gtext('電話費')gtext('衣服')gtext('其它')title('二維餅圖')subplot(2,1,2)colormap hsvpie3(x,explode)title('三維餅圖')8、畫出函數(shù)z=（x-2）2+（y-1.2）2+si

10、n(xy)的三維曲線和網(wǎng)格曲線。程序為：x,y=meshgrid(0:0.5:10); %為三維繪圖產(chǎn)生x，y數(shù)據(jù)矩陣z=(x-2).2+(y-1.2).2;figure(1)subplot(2,1,1)mesh(x,y,z),grid on %繪制網(wǎng)格曲線title('網(wǎng)格曲線')subplot(2,1,2)plot3(x,y,z),grid ontitle('三維曲線')9、畫出下列函數(shù)的曲面及等高線圖z=x2+y2+sin(xy)。程序為：x,y=meshgrid(0:pi/10:2*pi);z=x.2+y.2+sin(x*y);figure(1)subp

11、lot(2,1,1)surfc(x,y,z), grid ontitle('曲面和等高線')subplot(2,1,2)c,h=contour(x,y,z);set(h,'showtext','on','textstep',get(h,'levelstep')*2);title('等高線')184ye1、編寫M 函數(shù)實現(xiàn)：求一個數(shù)是否為素數(shù)，在編寫一主程序，要求通過鍵盤輸入一個整數(shù)，然后完成判斷其是否為素數(shù)。解：function prime(x)n=fix(sqrt(x);for i=2:nif r

12、em(x,i)=0a='fasle'returnelse a='true'endend4、function x1,x2=jfc(a,b,c)d=b2-4*a*c;if d>0; x1=(-b-sqrt(d)/(2*a); x2=(-b+sqrt(d)/(2*a);elseif d=0; x1=-b/(2*a); x2=x1;enda=input('a=');b=input('b=');c=input('c=');x1,x2=jfc(a,b,c);if x1=x2; disp('原方程存在兩個不同的根&

13、#39;) disp('x1=',num2str(x1) disp('x2=',num2str(x2);elseif x1=x2; disp('原方程存在兩個相同的根') disp('x1=x2=',num2str(x1)else disp('原方程的根不存在')end%MATLAB上機作業(yè)，184頁習(xí)題4，求二次方程的實根function s1,s2=solve(a,b,c)d=b2-4*a*c;if d>0 s1=(-b-sqrt(d)/(2*a) s1=(-b-sqrt(d)/(2*a) %disp(&

14、#39;原方程有兩個不同的根'num2str(s1)'和'nums2tr(s2)elseif d=0 s1=(-b-sqrt(d)/(2*a); s1=(-b-sqrt(d)/(2*a); %disp('原方程有兩個相同的根'num2str(s1)else disp('原方程無實根')endend5、編寫程序計算x 在（-3,3）上，并畫出曲線。解：function y=func2(x)n=length(x);for i=1:n;if (x(i)>=-3)&&(x(i)<-1)y(i)=-x(i).2-4*x(

15、i)-3/2;elseif (x(i)>=-1)&&(x(i)<1)y(i)=-x(i).2+1;else (x(i)>=1)&&(x(i)<3)y(i)=-x(i).2+4*x(i)-3/2;endend腳本為：x=-3:.01:3;y=func2(x);figure(1)plot(x,y),grid ontitle('y=func2(x)')xlabel('x'), ylabel('y')第二種方法：5、x=-3:0.01:3;if x>=-3 & x<-1; y=(

16、-x.2-4*x-3)/2;elseif x>=-1 & x<1; y=-x.2+1;else y=(-x.2+4*x-3)/2;endplot(x,y)6、利用menu 函數(shù)輸入選擇參數(shù)ch。當(dāng)ch=1 時，產(chǎn)生-10,10之間均勻分布的隨機數(shù)；當(dāng)ch=2 時，產(chǎn)生-5,5之間均勻分布的隨機數(shù)；當(dāng)ch=3 時，產(chǎn)生-1,1之間均勻分布的隨機數(shù)；當(dāng)ch=4 時，產(chǎn)生均值為0，方差為1 的正態(tài)分布隨機數(shù)。要求使用switch 函數(shù)。解：s=menu('ch','1','2','3','4');n=

17、;switch scase 1,n=20*rand(3)-10case 2,n=10*rand(3)-5case 3,n=2*rand(3)-1case 4,n=randn(3)otherwise disp('error')end8、企業(yè)發(fā)放的獎金按個人完成的利潤（I）提成。分段提成比例I K 為I K = 1 0 % , 1 05 % 2 02 % 0 4 01 % 4 0IIII < < > 萬元， 1 0 萬元， 2 萬元， 萬元即如王某完成25 萬元利潤時，個人可得y=10 x 10%+ 10 x 5% + 5 x 2%(萬元)據(jù)此編寫程序，求企業(yè)職工

18、的獎金。解：function bonus=bon(I)n=fix(I/100000)if(n>4)n=4;endbon1=100000*0.1;bon2=0.05*(200000-100000);bon3=0.02*(400000-200000);switch ncase 0,bonus=I*100000;case 1bonus=bon1+0.05*(I-100000);case 2,3bonus=bon1+bon2+0.02*(I-200000);case 4,bonus=bon1+bon2+bon3+0.01*(I-400000);end運行后：>> I=1700000;>> bonus=bon(I)n =179、有一分?jǐn)?shù)序列2/1，3/2，5/3/，8/5求前15 項和。解：s=1;t=2;sum=0;x=t/s;sum=sum+x;for i=1:15z=t;t=s+t;s=z;x=t/s;sum=sum+x;