matlab平時作業(yè)

1、MATLAB程序設計基礎平時作業(yè)說明：1-17題為必做題，須分別提交紙質文檔及與紙質文檔對應的電子文檔.doc版。1. 對于每輸入的一個百分制成績，要求輸出對應的成績等級A、B、C、D、E。其中A的分值在區(qū)間90,100)內；B的分值在區(qū)間80,90)內；C的分值在區(qū)間70,80)內；D的分值在區(qū)間60,70)內；E的分值在區(qū)間0,60)內；要求：1）要有對輸入成績是否合理的判斷：輸入成績如果含有小數(shù)，只能是“x.5”的形式，其實x為合理的整數(shù)。成績必須在百分制范圍內。如果輸入的成績不合理，應向命令窗口輸出報錯信息“輸入的成績不合理”2）用if語句實現(xiàn)題目要求x=input('輸入一個

2、數(shù)值:'); %輸入if mod(x,0.5)=0&&x>=0&&x<100; %mod(x,0.5)表示x除以0.5的余數(shù),用以判斷成績含小數(shù)是否為x.5 if x<0&x>=100 %不合理成績disp('輸入的成績不合理') elseif x<100&x>=90%成績90 100)disp('成績等級為：A')elseif x<90&x>=80 %成績80 90)disp('成績等級為：B') elseif x<80&

3、;x>=70%成績70 80) disp('成績等級為：C') elseif x<70&x>=60%成績60 70)disp('成績等級為：D')elsedisp('成績等級為：E')%成績0 60) endelse disp('輸入的成績不合理') end運行結果：輸入一個數(shù)值:57成績等級為：E輸入一個數(shù)值:93.3輸入的成績不合理3）用switch語句實現(xiàn)題目要求x=input('請輸入一個成績:'); %輸入if mod(x,0.5)=0&&x>=0&

4、&x<100 %mod(x,0.5)表示x除以0.5的余數(shù),用以判斷成績含小數(shù)是否為x.5 switch floor(x/10) %x除以10向下取整 case 9 disp('成績等級為：A'); case 8 disp('成績等級為：B'); case 7 disp('成績等級為：C'); case 6 disp('成績等級為：D'); case 5,4,3,2,1,0 disp('成績等級為：E'); otherwise disp('輸入的成績不合理'); endelse dis

5、p('輸入的成績不合理');end運行結果：請輸入一個成績:0成績等級為：E請輸入一個成績:120輸入的成績不合理請輸入一個成績:95.5成績等級為：A2. 設計程序，完成10,99之間兩位整數(shù)的加、減、乘、除四則運算。要求：1）隨機產(chǎn)生區(qū)間10,99之內的兩位隨機整數(shù)2）在產(chǎn)生符合要求的隨機數(shù)后，輸入一個運算符號，做相應的運算3）運算結果顯示在命令窗口，其顯示類似于“a=x+y=68”的形式。x=round(99-10)*rand(1)+10); %在區(qū)間10 99隨機選一個數(shù)x y=round(99-10)*rand(1)+10); %在區(qū)間10 99隨機選一個數(shù)ysym=

6、input('請輸入運算符號(+-*/):','s'); switch sym %switch語句 case'+' a=x+y; disp('輸入的兩數(shù)之和b=',num2str(x),'+',num2str(y),'=',num2str(a); %輸出顯示類似于“a=x+y=68”的形式 case'-' a=x-y; disp('輸入的兩數(shù)之差c=',num2str(x),'-',num2str(y),'=',num2str(a);

7、case'*' a=x*y; disp('輸入的兩數(shù)之積d=',num2str(x),'*',num2str(y),'=',num2str(a); case'/' a=x/y; disp('輸入的兩數(shù)之商e=',num2str(x),'/',num2str(y),'=',num2str(a); otherwise disp('輸入符號錯誤') end運行結果：請輸入運算符號(+-*/):+輸入的兩數(shù)之和b=95+31=126請輸入運算符號(+-*/):

8、-輸入的兩數(shù)之差c=64-53=11請輸入運算符號(+-*/):*輸入的兩數(shù)之積d=89*78=6942請輸入運算符號(+-*/):/輸入的兩數(shù)之商e=51/12=4.253. 求下列分段函數(shù)的值并繪制函數(shù)曲線。 要求：1）x的值以向量的形式從鍵盤輸入。2）用if語句實現(xiàn)，分別輸出x =5.0，3.0，1.0，2.0，2.5，3.0，5.0時的y值。x=input('please input x=');if x<0 & x=3 y=x.2+x-6;elseif x>=0 & x<10 & x=2 & x=3 y=x.2-5.*x

9、+6;else y=x.2-x-1;enddisp('y=',num2str(y)運行結果：please input x=-5.0 -3.0 1.0 2.0 2.5 3.0 5.0 y=14 11 -1 1 2.75 5 193）自行擬定自變量的取值，繪制該分段函數(shù)的二維曲線。x=-10:1:15; %自擬自變量的取值if x<0 & x=3 y=x.2+x+6;elseif x>=0 & x<10 & x=2 & x=3 y=x.2-5.*x+6;else y=x.2-x-1;endplot(x,y,'r-')

10、 %繪制分段函數(shù)title('分段函數(shù)')xlabel('x')ylabel('y')hold on 圖3（3）分段函數(shù)的二維曲線4. 分別使用for循環(huán)語句和while循環(huán)語句找出最小的n值，使得n!>10100,并求出n!。for循環(huán)語句如下：s=1; %對s賦初值for n=1:100 %默認步長為1 s=s*n; %s=n! if(s>10100) break; %如果s>10100,結束運行 endenddisp('n=',num2str(n) %輸出n disp('s=',num2st

11、r(s) %輸出s運行結果：n=8s=40320while循環(huán)語句如下：s=1; %對s賦初值n=0; %循環(huán)次數(shù)賦初值while(s<=10100) n=n+1; %不能使用+ s=s*n; enddisp('n=',num2str(n);disp('s=',num2str(s); %組合字符串 運行結果：n=8s=403205. 已知f1=1，n=1f2=0，n=2f3=1，n=3fn=fn12fn2fn3，n>3求f1 f50中：1） 最大值和最小值及它們的位置、各數(shù)之和（最大值給變量MAX，其位置給變量x，最小值給變量MIN，其位置給變量y，

12、各數(shù)之和給變量SUM）。2） 正數(shù)、零、負數(shù)的個數(shù)（依次賦值給變量N_positive，N_zero，N_negative）。f(1)=1; f(2)=0;f(3)=1; %輸入n=3; %賦值while (n<50) n=n+1; f(n)=f(n-1)-2*f(n-2)+f(n-3);endMAX,x=max(f) %MAX為向量f中最大值元素，x為最大值元素的索引號MIN,y=min(f) %MAX為向量f中最小值元素，y為最小值元素的索引號SUM=sum(f)N_positive=numel(find(f>0) %numel表示返還數(shù)組中的元素個數(shù)，N_zero=numel

13、(find(f=0) %返還f中等于0的位置N_negative=numel(find(f<0)運行結果：MAX =406631x = 49MIN =-403795y =50SUM =76861N_positive =25N_zero =2N_negative =236. Fibonacci數(shù)組的元素滿足Fibonacci規(guī)則：ak+2=ak+ak+1,（k=1，2，3，）；且a1= a2=1。請設計一段程序，求出該數(shù)組中第一個大于10000的元素對應的k值以及該元素值。a(1)=1; a(2)=1;k=1;while a(k)<=10000 a(k+2)=a(k)+a(k+1);

14、 k=k+1;enddisp('k=',num2str(k)disp('a(k)=',num2str(a(k)運行結果：k=21a(k)=109467. 定義自變量x的取值區(qū)間為5,5，編寫一個函數(shù)文件，實現(xiàn)繪制函數(shù)y=ax3+bx2+cx+d的曲線，其中給定各系數(shù)值分別為a=1，b=2， c=5，d=3;a=3， b=7，c=1， d=4;a=0， b=4， c=0， d=5 要另外編寫調用該函數(shù)文件的m文件，實現(xiàn)調用。function my1(a,b,c,d) %定義主函數(shù)，函數(shù)文件名為my1.mx=-5:1:5; %輸入xy=a*x.3+b*x.2+c*x

15、+d;plot(x,y); %繪制曲線xlabel('x軸'); %標志x軸ylabel('y軸'); %標志y軸title('函數(shù)圖像'); %設置圖名稱a=input('please input a='); %子函數(shù)b=input('please input b=');c=input('please input c=');d=input('please input d='); %利用輸出參數(shù)a,b,c,d向my1傳遞數(shù)據(jù)my1(a,b,c,d); %調用的是啊,a,b,c,d四個

16、參數(shù)圖7（1） 圖7（2）圖7（3） 8. 設計程序，在一個圖形窗體中畫出函數(shù)y=ex和y=logx的特性曲線。其中y=ex中x取值區(qū)間為(2,2)，y=logx中x取值區(qū)間為(0.1,5)。要求：1）要有圖例、標題（'二維圖'）、坐標軸標簽（'x軸數(shù)據(jù)'和'y軸數(shù)據(jù)'）；2）每條曲線的屬性自己設置，但不能全用默認值。x1=-2:0.1:2;y1=exp(x1);x2=0.1:0.1:5;y2=log(x2);plot(x1,y1,'r*-',x2,y2,'ko-')legend('函數(shù)y=exp(x)的特

17、性曲線','函數(shù)y=logx的特性曲線')title('二維圖')xlabel('x軸數(shù)據(jù)')ylabel('y軸數(shù)據(jù)')圖89. 采用2行2列子圖的方式進行二維圖形繪制，要求：1）自變量x的取值區(qū)間為0,22）三個函數(shù)分別為：y1=sin(x)，y2=cos(x)，y3=sin(x)+cos(x)3）在第1幅子圖繪制y1的曲線；在第2幅子圖繪制y2的曲線；在第3幅子圖繪制y3的曲線；在第4幅子圖中繪制所有3條曲線4）對每一幅子圖，都要有圖例、坐標軸標簽；曲線的屬性自己設置，不能全用默認值。x=0:0.1:2*pi;y1=

18、sin(x);y2=cos(x);y3=sin(x)+cos(x);subplot(2,2,1);plot(x,y1,'ro-');legend('sin曲線');title('y=sin(x)');xlabel('橫坐標1');ylabel('縱坐標1');subplot(222);plot(x,y2,'g*-');legend('cos曲線');title('y=cos(x)');xlabel('橫坐標2');ylabel('縱坐標2

19、9;);subplot(223);plot(x,y3,'b-p');legend('sin(x)+cos(x)曲線');title('y=sin(x)+cos(x)');xlabel('橫坐標3');ylabel('縱坐標3');subplot(224);plot(x,y1,'ro-',x,y2,'g*-',x,y3,'b-p'); %同時繪制三個函數(shù)的圖像legend('sin曲線','cos曲線','sin(x)+cos(x

20、)曲線');title('三個函數(shù)的曲線');xlabel('橫坐標');ylabel('縱坐標');圖（9）10. 采用1行2列子圖的方式進行二維圖形繪制，要求： 1）在第1幅子圖繪制三維曲線：x0,20y=sin(x)z=cos(x)2）在第2幅子圖繪制三維曲面：x20,20y10,103） 對每一幅子圖，都要有標題、坐標軸標簽；曲線的屬性自己設置，不能全用默認值。subplot(121); %在第一個子圖中繪圖x1=0:0.1:20;y1=sin(x1);z1=cos(x1);plot(x1,y1,'ro-',x1,

21、z1,'b*-');legend('sin曲線','cos曲線'); %設置圖例title('圖1');xlabel('x1');ylabel('y1');subplot(122); %在第二個子圖中繪圖x2=-20:1:20;y2=0.5*x2;X,Y=meshgrid(x2,y2); %生成X和Y網(wǎng)絡坐標矩陣z=(X.2)/9-(Y.2)/4; %網(wǎng)格化后直接輸入zmesh(X,Y,z); %曲面繪圖legend('雙曲面');title('圖2');xlabel

22、('x2');ylabel('y2');zlabel('z')grid on;圖（10）11. 已知函數(shù)，其中x的取值區(qū)間為2,2，試編程繪制該二維曲線，要求：要有標題、坐標軸標簽，曲線顏色用紅色。 p=ezplot('(x.2+y.2-1).3-x.2.*y.3',-2 2); %在區(qū)間-2 2繪制(x2+y2-1)3-x2y2=0的圖像set(p,'color','r'); %指定曲線顏色為紅色title('曲線圖');xlabel('x軸');ylabel(&#

23、39;y軸');圖（11）12. 已知原始數(shù)據(jù)，插值點為，試用4種不同的插值方法進行插值，并在一幅圖的4個子圖中分別繪出4種插值曲線，并各自注明所使用的插值方法。x=0:1:4*pi;y=sin(x).*exp(-x/5);xi=0:0.3:4*pi; %設置新的插點subplot(221)yi=interp1(x,y,xi,'nearset'); %使用最近點插值plot(xi,yi)title('最近點插值')subplot(222)yi=interp1(x,y,xi); %默認為線性插值plot(xi,yi)title('線性插值'

24、)subplot(223)yi=interp1(x,y,xi,'cubic'); %使用三次立方插值plot(xi,yi)title('三次立方插值')subplot(224)yi=interp1(x,y,xi,'spline'); %使用三次樣條插值plot(xi,yi)title('三次樣條插值')圖（12）13考慮表1中的數(shù)據(jù)：表1 一周氣溫溫度數(shù)據(jù)周一周二周三周四周五周六周日早晨15151813131516中午20212519202225夜晚16171214161819完成以下功能：1）將一周的氣溫溫度數(shù)據(jù)，使用一個矩陣A

25、進行保存。2）取出周二中午、周四上午、周日中午和夜晚的溫度值，組成新的溫度值向量。3）取出矩陣A中的前六個數(shù)值，組成新的溫度值向量。4）取出周二、周三、周四及周六的數(shù)據(jù)，組成新的溫度值矩陣。5）求每天溫度的平均值6）求一周的早晨、中午和夜晚溫度的平均值7）求一周的溫度平均值及最高溫度及最低溫度A=15 15 18 13 13 15 16;20 21 25 19 20 22 25;16 17 12 14 16 18 19 %將一周的氣溫溫度數(shù)據(jù)使用A矩陣保存B=A(5 10 20 21) %取出周二中午周四上午周日中午和夜晚的溫度組成新矩陣C=A(1:1:6) %取出A中的前六個組成新的矩陣D=

26、A(:,2 3 4 6) %取出周二周三周四和周六的數(shù)據(jù)組成新的矩陣everyday=mean(A) %表示每一列的平均值，即每天氣溫的平均值mae=mean(A') %一周早上，中午，晚上的平均值all=mean(A(:) %一周的溫度平均值MAX=max(A(:) %一周的最高溫度MIN=min(A(:) %一周的最低溫度14考慮表2中的數(shù)據(jù)，這是一位舉重教練收集的不同性別的舉重動員每個年齡的最大舉重磅數(shù)，據(jù)信這三者之間存在著某種函數(shù)關系，請估算當前隊伍中隊員（要求1320歲每一歲都有）的最大舉重磅數(shù)。表2 不同性別的舉重動員每個年齡的最大舉重磅數(shù)性別年齡重舉大最131517192

27、0男270290330350380女250270310330360t=13 15 17 19 20'y=270 290 330 350 380;250 270 310 330 360'X=ones(size(t) t t.2;%ones（size（t）以t的長度為矩陣的大小設置全1矩陣Z=Xy; %得到系數(shù)T=(13:20)'Y=ones(size(T) T T.2*Z;MAXT=max(Y)運行結果：MAXT =14305/38 13545/3815.某鄉(xiāng)鎮(zhèn)企業(yè)2009-2015年的生產(chǎn)利潤如表3所示。表3 過去7年生產(chǎn)利潤年份2009201020112012201320142015利潤（萬元）70122144152174196202要求用多項式擬合法預測出該企業(yè)2016年和2017年的利潤，并繪出擬合圖。x=2009:2015; %x在2009 2015取值，間隔為1y=70 122 144 152 174 196 202; %得到數(shù)據(jù)（x,y)p=polyfit(x,y,1); %對（x,y)數(shù)據(jù)進行1階多項式擬合，得到多項式系數(shù)px1=2009:2018;y1=polyval(p,x1); %根據(jù)得到的多項式，計算x處的