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

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

2、t;-0.5);C=sum(sum(B)運(yùn)行結(jié)果（每次運(yùn)行結(jié)果是不同的，僅作參考）：C=488、有一矩陣A，找出矩陣中值等于1的元素，并將它們重新排列成列向量B。解：A=2*rand(4);k=find(A<=1);A(k)=;%刪除下標(biāo)為k的元素B=A'9、在一測(cè)量矩陣A（100*3）中，存在有奇異值（假設(shè)大于100的置認(rèn)為是奇異值），編程實(shí)現(xiàn)刪去奇異值所在的行。解：A=120*randn(10,3);i,j=find(A>100);A(i,:)= %刪去存在奇異值的行運(yùn)行結(jié)果（每次運(yùn)行結(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的行運(yùn)行結(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、將圖形窗口分成兩個(gè)，分別繪制正割和余割曲線，并加上標(biāo)注。程序?yàn)椋簒1=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，在半對(duì)數(shù)坐標(biāo)系中繪制曲線。程序?yàn)椋簒=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 %半對(duì)數(shù)坐標(biāo)軸legend('ity=exp(x+5)+x.3')title('半對(duì)數(shù)坐標(biāo)')xlabel('x'), ylabel('y')5、將圖形窗口分成兩個(gè)，分別繪制以下函數(shù)在-3,3區(qū)間上的曲線，并利用axis 調(diào)整軸刻度，使他們具有相同縮放尺度。y1=2x+5； y2=x2-3x+1。程序?yàn)椋簒=-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、制餅圖。程序?yàn)椋簒=190 33 45 42 45;explode=0 1 0 0 0;figure(1)subplot(2,1,1)colormap hsvpie(x,explode)gtext('生活費(fèi)')gtext('資料費(fèi)')gtext('電話費(fèi)')gtext('衣服')gtext('其它')title('二維餅圖')subplot(2,1,2)colormap hsvpie3(x,explode)title('三維餅圖')8、畫(huà)出函數(shù)z=（x-2）2+（y-1.2）2+si

10、n(xy)的三維曲線和網(wǎng)格曲線。程序?yàn)椋簒,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、畫(huà)出下列函數(shù)的曲面及等高線圖z=x2+y2+sin(xy)。程序?yàn)椋簒,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、編寫(xiě)M 函數(shù)實(shí)現(xiàn)：求一個(gè)數(shù)是否為素?cái)?shù)，在編寫(xiě)一主程序，要求通過(guò)鍵盤(pán)輸入一個(gè)整數(shù)，然后完成判斷其是否為素?cái)?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('原方程存在兩個(gè)不同的根&

13、#39;) disp('x1=',num2str(x1) disp('x2=',num2str(x2);elseif x1=x2; disp('原方程存在兩個(gè)相同的根') disp('x1=x2=',num2str(x1)else disp('原方程的根不存在')end%MATLAB上機(jī)作業(yè)，184頁(yè)習(xí)題4，求二次方程的實(shí)根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;原方程有兩個(gè)不同的根'num2str(s1)'和'nums2tr(s2)elseif d=0 s1=(-b-sqrt(d)/(2*a); s1=(-b-sqrt(d)/(2*a); %disp('原方程有兩個(gè)相同的根'num2str(s1)else disp('原方程無(wú)實(shí)根')endend5、編寫(xiě)程序計(jì)算x 在（-3,3）上，并畫(huà)出曲線。解：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 時(shí)，產(chǎn)生-10,10之間均勻分布的隨機(jī)數(shù)；當(dāng)ch=2 時(shí)，產(chǎn)生-5,5之間均勻分布的隨機(jī)數(shù)；當(dāng)ch=3 時(shí)，產(chǎn)生-1,1之間均勻分布的隨機(jī)數(shù)；當(dāng)ch=4 時(shí)，產(chǎn)生均值為0，方差為1 的正態(tài)分布隨機(jī)數(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ā)放的獎(jiǎng)金按個(gè)人完成的利潤(rùn)（I）提成。分段提成比例I K 為I K = 1 0 % , 1 05 % 2 02 % 0 4 01 % 4 0IIII < < > 萬(wàn)元， 1 0 萬(wàn)元， 2 萬(wàn)元， 萬(wàn)元即如王某完成25 萬(wàn)元利潤(rùn)時(shí)，個(gè)人可得y=10 x 10%+ 10 x 5% + 5 x 2%(萬(wàn)元)據(jù)此編寫(xiě)程序，求企業(yè)職工

18、的獎(jiǎng)金。解：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運(yùn)行后：>> I=1700000;>> bonus=bon(I)n =179、有一分?jǐn)?shù)序列2/1，3/2，5/3/，8/5求前15 項(xiàng)和。解：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;