matlab簡單編程21個題目及問題詳解

1、實用文檔cl 3sin xy = cos x 0.5 -1、 設(shè)-(1+X )，把x=02冗間分為101點，畫出以x為橫坐標(biāo)，y為縱坐標(biāo)的曲線。第一題的matlab源程序： 考慮cos (x)為一個整體，然后乘以中括號里面的全部x=0:2*pi/100:2*pi;%x的步長以及范圍從 0到2*piy=cos(x).*(0.5+3*sin(x)./(1+xA2);%y 的表達(dá)式plot(x,y)%畫出圖形圖如下：考慮對整體求解cos,先求x乘以括號中的部分x=0:2*pi/100:2*pi;%x的步長以及范圍從 0至U2*piy=cos(x.*(0.5+3*sin(x)./(1+x.A2);%y

2、 的表達(dá)式plot(x,y)%畫出圖形文案大全圖如下:R1,求其各列的平均值和均方差。并求%產(chǎn)生正態(tài)分布隨機(jī)矩陣2、產(chǎn)生8X6階的正態(tài)分布隨機(jī)數(shù)矩陣 該矩陣全體數(shù)的平均值和均方差。第二題的matlab源程序如下：R1=randn(8,6)R1 =aver=(sum(R1(1:end,1:end)./8 aver =0.07680.13630.1022a=std(R1(1:end,1:end)a =1.08190.80931.3456aver1=(sum(R1(:)./48aver1 =%產(chǎn)生各行的平均值-0.34730.49750.1044%產(chǎn)生各行的均方差也就是標(biāo)準(zhǔn)差0.82330.8079

3、1.2150%全體數(shù)的平均值1.0933-0.76971.5442-0.19241.41930.21571.10930.37140.08590.88860.2916-1.1658-0.8637-0.2256-1.4916-0.76480.1978-1.14800.07741.1174-0.7423-1.40231.58770.1049-1.2141-1.0891-1.0616-1.4224-0.80450.7223-1.11350.03262.35050.48820.69662.5855-0.00680.5525-0.6156-0.17740.8351-0.66691.53261.10060.

4、7481-0.1961-0.24370.18730.0950 b=std(R1(:) b =%全體數(shù)的均方差即標(biāo)準(zhǔn)差1.0103 3、設(shè) x=rcost+3t,y=rsint+3, 分別令 r=2,3,4，畫出參數(shù) t=010 區(qū)間生成的 xy 曲線。第三題的matlab源程序如下: t=0：0.1：10； r=2;x=r*cos(t)+3*t;y=r*sin+3; plot(x,y,'r-') hold on t=0:0.1:10;r=3;x=r*cos(t)+3*t;y=r*sin+3; plot(x,y,'b-.') hold ont=0:0.1:10;r

5、=4;x=r*cos(t)+3*t;y=r*sin+3;%t的取值以0.1的步長從0到10%r=2%x的表達(dá)式%y的表達(dá)式%用紅色線繪制曲線%保持窗口巾gure窗口不關(guān)閉%t的取值以0.1的步長從0到10%r=3%x的表達(dá)式%y的表達(dá)式%用藍(lán)色點劃線繪制曲線%保持窗口巾gure窗口不關(guān)閉%t的取值以0.1的步長從0到10%r=4%x的表達(dá)式%y的表達(dá)式用紅紫色線繪制曲線plot(x,y,'m-')%圖形如下:76532100510162025304、設(shè) f(x)=x 5- 4x 4 +3x 2- 2x+ 6(1)在x=-2,8之間取100個點，畫出曲線，看它有幾個過零點。(提示

6、:用polyval 函數(shù))(2)用roots函數(shù)求此多項式的根。(1)matalab的源程序如下：(利用的是在同一圖中展現(xiàn)兩條線，通過放大看交點)x=-2:0.1:8;%x的取值范圍y=0;%y=0plot(x,y,'m')%用紅紫色表示繪制出的y=0的圖像hold on%保持figure窗口不關(guān)閉x=-2:0.1:8;%x的取值范圍p=1,-4,0,3,-2,6;%p表示函數(shù)表達(dá)式的系數(shù)矩陣y=polyval(p,x),plot(x,y) %y=polyval表示求多項式x處的y值，繪制圖像圖如下:將圖形進(jìn)行放大（如下），發(fā)現(xiàn)在實數(shù)域上有3個零點(2) matlab源程序如下

7、：roots(p)%求函數(shù)的根ans =3.7999-1.26071.34790.0564 + 0.9623i0.0564 - 0.9623i可見實數(shù)解有三個，復(fù)數(shù)解有兩個。5、建立一個字符串向量，然后對該向量做如下處理: (1)取第15個字符組成的子字符串。(2)將字符串倒過來重新排列。(3)統(tǒng)計字符串中小寫字母的個數(shù) 整個題的源程序如下所示：ch='abcj245DFGK'subch=ch(1:5)%建立一個字符串%取15個字符subch = abcj2revch=ch(end:-1:1)revch =KGFD542jcbaa=find(ch>='a'

8、&ch<='z');length(a)ans =46、設(shè) x=sint, y=sin(nt+a) , (1)若 a=1,令 n =1,2,3,4, (2)若 n=2,取 a=0,兀 /3 , (1) matlab源程序如下;t=-pi:0.1:pi;a=1;n=1;x=sin(t);y=sin(n*t+a);subplot(221),plot(x,y)hold ont=-pi:0.1:pi;a=1;n=2;x=sin(t);y=sin(n*t+a);subplot(222),plot(x,y)hold ont=-pi:0.1:pi;a=1;n=3;x=sin(t)

9、;y=sin(n*t+a);subplot(223),plot(x,y)hold ont=-pi:0.1:pi;a=1;n=4;%倒序排列%找到小寫字母的位置 %輸出小寫字母的個數(shù)在四個子圖中分別畫出其曲線。冗/2 ,及冗，在四個子圖中分別畫出其曲線%t的取值范圍%a=1, n=1%x的表達(dá)式%y的表達(dá)式%在子圖第一張，畫出圖像%保持figure不關(guān)閉%t的取值范圍%a=1, n=2%x的表達(dá)式%y的表達(dá)式%在子圖的第二張畫出圖像%保持figure不關(guān)閉%t的取值范圍%a=1, n=3%x的表達(dá)式%y的表達(dá)式%在子圖的第三張，畫出圖像%保持figure不關(guān)閉%t的取值范圍%a=1, n=4%x

10、的表達(dá)式x=sin(t);y=sin(n*t+a);subplot(224),plot(x,y)%y的表達(dá)式%在子圖的第四張上，繪制圖像圖形如下所示：(2) matlab源程序如下所示: t=-pi:0.1:pi;a=0;n=2;x=sin(t);y=sin(n*t+a);subplot(221),plot(x,y)hold ont=-pi:0.1:pi;a=pi./3;n=2;x=sin(t);y=sin(n*t+a);subplot(222),plot(x,y)hold ont=-pi:0.1:pi;a=pi./2;n=2;x=sin(t);y=sin(n*t+a);subplot(223

11、),plot(x,y)hold ont=-pi:0.1:pi;a=pi;n=2;%t的取值%a與n的取值%x的表達(dá)式%y的表達(dá)式%在子圖的第一張繪制圖像%保持日gure窗口不關(guān)閉%t的取值%a與n的取值%x的表達(dá)式%y的表達(dá)式%在子圖的第二張繪制圖像%保持日gure窗口不關(guān)閉%t的取值%a與n的取值%x的表達(dá)式%y的表達(dá)式%在子圖的第三張繪制圖像%保持日gure窗口不關(guān)閉%t的取值范圍%a與n的取值x=sin； y=sin(n*t+a);%x的表達(dá)式%y的表達(dá)式%在四個字圖的第四張畫出圖像subplot(224),plot(x,y)22、7、繪制z =x2e<x*）在定義域x=-2,2,

12、y=-2,2 內(nèi)的曲面Matlab源程序如下： x,y=meshgrid(-2:0.1:2); z=x.A2.*exp(-(x.A2+y.A2); surf(x,y,z)%此處描述x, y的范圍%z的表達(dá)式%繪制曲面圖形如下所示:經(jīng)過旋轉(zhuǎn)之后的圖形:0.40.350.30.250 2D.150.1口050(_x2 4x -3)/2-3 < x < -18、y = «x2 +1-1 <x <1 ，畫出在匚3,3】上的曲線2一 -一(-x +4x3)/21 <x<3所建立的M文件如下所示: function y=A(x);if x>=-3&

13、;x<-1y=(-xA2-4*x-3)./2;elseif x>=-1&x<1y=-x.A2+1;else x>=1&x<=3y=(-x.A2+4*x-3)./2;%文件名為A%在-3x<-1的條件下%y所對應(yīng)上述條件下的表達(dá)式%-1 «x<1%y所對應(yīng)的上述條件下的表達(dá)式%1 «x «3end在 command windows 中輸入:%y所對應(yīng)的上述條件下的表達(dá)式%結(jié)束fplot('A',-3,3) % 繪制在-3到3上的曲線所得圖形如下所示：syms xy=x.*exp(2.*x);s

14、=diff(y,x,20)9、(1)設(shè)¥=£產(chǎn)，求d網(wǎng)%定義變量x%y的表達(dá)式%對y求x的20階導(dǎo)數(shù)s =10485760*exp(2*x) + 1048576*x*exp(2*x)(2)設(shè)貫=兌?工十ex%,求一syms x%定義變量xy=sin(x).A4+cos(x).A4;%y 的表達(dá)式s=diff(y,x,10)%對 y 求 x 的 10 階導(dǎo)數(shù)s =-262144*cos(x)A4 + 1572864*cos(x)A2*sin(x)A2 - 262144*sin(x4現(xiàn)(3)設(shè)工=/與粉,十工外。節(jié)，求去syms x y%定義自變量x, yz=y.A4.*(si

15、n(x).A4)+x.A4.*(cos(y).A4);%z 的表達(dá)式s=diff(z,x)%對z求x的偏導(dǎo)數(shù)s =4*xA3*cos(y)A4 + 4*yA4*cos(x)*sin(x)A3(4)(x2 1)(x2 x)dxsyms xi=int(1/(x.A2+1).*(x.A2+x);pretty(i)%定義自變量x%不定積分表達(dá)式%求解不定積分/ i log(x + 1)log(x) + log(x + i) | 1/4 | -+ log(x - i) | - - 1/4 |(5)(9-x2)3dxsyms xi=int(sqrt(9-xA2)A3)./(xA6); pretty(i)/

16、 923 1/2I - - 1 I (- (x - 9)I 2 | x /%定義自變量x%不定積分表達(dá)式%求解不定積分(6)2 dx1 x x3345 xsyms xi=int(1./(x+x.A3),2,1) i =log(10A(1/2)/4)%定義自變量x%定積分表達(dá)式%定積分的解(7)求級數(shù)的和：n 2 nm 2nsyms n%定義自變量nf=(n+2)./(2.An);%級數(shù)的表達(dá)式symsum(f,n,n-1,inf)%對級數(shù)求和ans = 2*2A(1 - n)*(n + 2)(8)利用函數(shù)int計算二重不定積分看即syms x yi=int(int(x+y)*exp(-x.*y

17、);pretty(i)%定義自變量x, y%二重不定積分表達(dá)式%求解二重不定積分y exp(x y),顯示所有小于n的質(zhì)數(shù)。%定義文件名為xin10%清屏，清內(nèi)存%輸入數(shù)據(jù)%i的變化范圍%判斷是否為質(zhì)數(shù)%輸出質(zhì)數(shù)%結(jié)束if%結(jié)束for% 結(jié)束 function0時結(jié)束輸入，求這些數(shù)的平均值以及平方10、編制m文件，輸入n （正整數(shù)） 編制的m文件如下所示：function xin10;clc,clearn=input（'輸入一個正整數(shù)：'）for i=1:1:nif isprime（i）disp（i）endendend11、從鍵盤輸入若干個數(shù)，當(dāng)輸入和。編制的m文件如下所示：假

18、設(shè)輸入數(shù)據(jù)結(jié)束時輸入0作為結(jié)束，不考慮輸入的數(shù)據(jù)中出現(xiàn)會結(jié)束輸入的現(xiàn)象，則程序如下:%定義文件名為lin11%清屏，清內(nèi)存function lin11;clc,cleark=size(m) for i=1:k(2)if m(i)=0 k(2)=i-1;end end M=m(1:k(2) B=sum(M)./k(2) C=sum(MA2) endm=input（'以1*n矩陣的形式輸入數(shù)據(jù)并以0結(jié)尾：）輸入數(shù)據(jù)%m的大小%i的變化范圍%m=0時的每一個k值%另卜（2） =i-1%結(jié)束if循環(huán)%結(jié)束for循環(huán)%將若干個數(shù)輸出%對若干個數(shù)求和再求平均值%輸出m的值，進(jìn)行求值% 結(jié)束 fun

19、ction編制的m儀件如下所示:function lin12;clc,cleara=rand(10,1)for i=1:length(a)12、利用rand產(chǎn)生10個隨機(jī)數(shù)，利用for循環(huán)對其進(jìn)行排序（從大到?。?定義文件名為lin12%清屏，清內(nèi)存%產(chǎn)生10個隨機(jī)數(shù)，這個地方可加可不加；%i的取值范圍for j=1:length(a)-1if a(j)<a(j+1) t=a(j)； a(j)=a(j+1); a(j+1)=t；endendendaendfunction lin13;clc,cleara=input('請輸入數(shù)據(jù): n=length(a);for i=1:n/2b

20、=a(i);a(i)=a(n+1-i);a(n+1-i)=b;endfprintf(' %d',a) endfunction lin14;clc,cleara=input('請輸入數(shù)據(jù)：','s) n=length(a);fprintf('位數(shù) %d',n);end%定義文件名為lin15%清屏，清內(nèi)存%輸入a的值%輸入b的值%定義i, sum, he的初始值%i的變化范圍%通過i值得變化得到不同的數(shù)值%令sum為總和%i的變化%j的取值范圍%冒泡法進(jìn)行排序，判斷兩個數(shù)的大小%將小的數(shù)賦值給t%將大數(shù)賦值給a (j)%將1賦值給a (j+

21、1)%結(jié)束if循環(huán)%結(jié)束for循環(huán)%結(jié)束for循環(huán)%輸出數(shù)據(jù)命匕處一定要寫上a,不然沒有數(shù)據(jù)輸出% 結(jié)束 function13、將一個整型數(shù)組的元素按逆序重新存放(如原序為：8, 6, 5, 4改為：4, 5, 6, 8)。編制的m文件如下所示:%定義文件名為lin13%清屏，清內(nèi)存')%輸入數(shù)據(jù)%判斷數(shù)據(jù)的長度%i的變化范圍%將2賦值給b%將對應(yīng)位置的值復(fù)制給a (i)%將b值賦值給a (n+1-i)%結(jié)束for循環(huán)%這個地方d之前有一個空格，保證輸出數(shù)據(jù)時分開的% 結(jié)束 function14、輸入一個整數(shù)，寫一程序輸出它是幾位數(shù)編制的m文件如下所示：%定義文彳名為lin14%清屏，

22、清內(nèi)存%此處由于輸入的整數(shù)是以整體的形式輸入，所以s'%計算a的長度%輸出a的位數(shù)% 結(jié)束 function 15、從鍵盤上輸入 a與n的值，計算sum=a+aa+aaa+aaaa+(共n項)的和。 例 a=2, n=4,則 sum=2+22+222+2222.編制的m文件如下所示：function lin15;clc,cleara=input('請輸入 a=) n=input('請輸入 n=') i=1;sum=0;he=0;while i:nhe=he+a*10A(i-1); sum=sum+he;i=i+1;end% 結(jié)束 whilefprintf(

23、9;sum is %ld',sum);%輸出總和end% 結(jié)束 function16、設(shè)計一個用于計算個人所有稅的程序。假設(shè)個人所得稅的繳納標(biāo)準(zhǔn)為：月 收入少于等于800元者不納稅；超出800元的部分，納稅5%超出2000元的部 分，納稅10%超出5000元的部分，納稅20%超出10000元的部分，納稅30% 超出100000元的部分，納稅40%編制的m文件如下所示：function xin16;%定義文件名為 xin16clc,clear%清屏，清內(nèi)存a=input('請輸入月收入：')%輸入a的值if a>100000%工資數(shù)目大于100000時j=(a-10

24、0000)*0.4+90000*0.3+5000*0.2+3000*0.1+1200*0.5; % 交稅計算公式elseif a>10000%工資數(shù)目大于10000且小于100000的條件下時j=(a-10000)*0.3+5000*0.2+3000*0.1+1200*0.5;% 交稅計算公式%交稅計算公式elseif a>5000%工資數(shù)目大于5000且小于10000時j=(a-5000)*0.2+3000*0.1+1200*0.5;elseif a>2000%工資數(shù)目大于2000且小于5000時j=(a-2000)*0.1+1200*0.05；%交稅計算公式elseif

25、a>800%工資數(shù)目大于800且小于2000時j=(a-800)*0.05; elsej=0 end fprintf(' %d',j) end%工資低于800時%不交稅%結(jié)束if循環(huán)%輸出交稅數(shù)目% 結(jié)束 function%交稅計算公式17、求愛因斯坦數(shù)學(xué)題。有一條長階，若每步跨 2階，則最后剩余1階；若每 步跨3階，則最后剩2階；若每步跨5階，則最后剩4階；若每步跨6階，則 最后剩5階；若每步跨7步，最后正好一階不剩。編制的m文件如下所示:%建立文件xin17%清屏，清內(nèi)存%2的取值變化范圍function xin17; clc,clearfor a=7:1000if

26、rem(a,2)=1&rem(a,3)=2&rem(a,5)=4&rem(a,6)=5&rem(a,7)=0%判斷條件就是題目的要求break%跳出if,再進(jìn)行循環(huán)end%結(jié)束if循環(huán)end%結(jié)束for循環(huán)fprintf('臺階is %d',a);%輸出臺階數(shù)目end% 結(jié)束 function通過計算結(jié)果顯示為119階臺階。18、將一個正整數(shù)n輸出成“千分位”形式，即從個位數(shù)起，每三位之間加一 個逗號，例如，將 7654321輸出成7,654,321.編制的m文件如下所示:function xin18;clc,cleara=input('

27、請輸入數(shù)據(jù)','s) n=length(a);b=；e二口;d=','；%建立文件名為xin18%清屏，清內(nèi)存%輸入數(shù)據(jù)%計算a的長度%令b為空集%令6為空集%令為，z=fix(n/3); %z為對三求商然后對零取整，z是添加幾個，的標(biāo)志if rem(n,3)=0;c=rem(n,3);for i=1:c%如果n對三取余不等于0%c為對三取余數(shù)之后的余數(shù)，為1或2%i的取值范圍b=b,a(i);%將2的最高位賦值給b,如果c=1, b=a最高位上的一位數(shù)，如果c=2, b=a最高位上的兩位數(shù) endif c=1%討論c=1的情況for i=1:z%i的變化范圍e

28、=e,d;%將逗號放到e中去for i=(3*i-1):(3*i+1)%i 的變化范圍e=e,a(i); %根據(jù)i隨z的變化，將2,3,4。5,6,7。等位置上的 數(shù)據(jù)放到e中的逗號后面end%結(jié)束for循環(huán)end%結(jié)束for循環(huán)。此處的for循環(huán)，就是輸入一個逗號,輸入三個數(shù)據(jù)，以此類推else c=2%討論c=2的情況for i=1:z%i的變化范圍e=e,d;%將逗號放到e中去for i=(3*i):(3*i+2) %i 的取值范圍e=e,a(i);%將3,4,5。6,7,8。等位置上的數(shù)據(jù)放到逗號后面end%結(jié)束for循環(huán)end%結(jié)束for循環(huán)end%結(jié)束if循環(huán)else%此處討論恰好

29、整除的情況for i=1:3%i的取值范圍b=b,a(i);%整除的情況下把前三位直接賦值到b中去end%結(jié)束for循環(huán)for i=1:(z-1)%i的取值范圍。由于恰好整除，所以前三位不用加逗號，所以逗號數(shù)為z-1e=e,d;%將逗號放到e中for i=3*i+1:3*i+3%i 的取值范圍e=e,a(i);%根據(jù)i隨z的變化,將4,5,6。7,8,9。等位置上的數(shù)據(jù)放到e中的逗號后面endend%結(jié)束for循環(huán)%結(jié)束for循環(huán)end f=b,e end%結(jié)束if循環(huán)%將b得到的前幾位與eW到的逗號數(shù)據(jù)放入f中去% 結(jié)束 function19、求出矩陣a中每行元素的平均值和平均最大的行號。編

30、制的m文件如下所示:function lin19;clc,cleara=input('輸入矩陣:');b=；i,j=size(a);for i=1:is=0;for j=1:j s=s+a(i,j);end%定義文件名為lin19%清屏，清內(nèi)存%輸入矩陣%定義b為空集%對2求行列數(shù)，放到i, j中%i的取值范圍%令$二0%j的取值范圍%對i 一定時對行累加，改變i值，不同的行重新進(jìn)行求和%結(jié)束for循環(huán)%輸出每行平均值fprintf('平均值(i=%d)=%dn',i,s/j)b=b,s/j;endw=0;for p=1:iif b(p)>w w=b(p)

31、;endendfprintf('平均最大行1end20、%平均值放入b中去%結(jié)束for循環(huán)%令0=0%p的取值范圍%每一行的平均值與w比較;%如果平均值大于w,將平均值賦值給w%結(jié)束if循環(huán)%結(jié)束for循環(huán)%dn',p)%輸出平均最大的行號% 結(jié)束 function考慮甬?dāng)?shù)山三一二.就在r = 3:0L8處生戊-綱梏£果用馥三次樣 條和B樣條分別棧據(jù)進(jìn)體合，并由數(shù)果求取二階導(dǎo)粼將鼬的給粘理論出自 (使用函數(shù)csapi和spapi及其相關(guān)函數(shù))x=3:0.4:8;%x 的取值y=(sqrt(1+x)-sqrt(x-1)./(sqrt(2+x)+sqrt(x-1); %y

32、 的表達(dá)式Sp=csapi(x,y)%得到分段三次樣條的有關(guān)數(shù)據(jù)Sp = form: 'pp' breaks: 3 3.4000 3.8000 4.2000 4.6000 5 5.4000 5.8000 6.2000 6.6000 7.0000 7.4000 7.8000 coefs: 12x4 double pieces: 12order: 4dim: 1Sp.coefs%顯示每個區(qū)間上三次多項式的系數(shù)ans =-0.00340.0169-0.05340.1605-0.00340.0127-0.04150.1416-0.00180.0086-0.03300.1268-0.00

33、130.0064-0.02700.1149-0.00090.0048-0.02250.1050-0.00070.0038-0.01900.0968-0.0005-0.0004-0.0003-0.0002-0.00020.00300.00240.00200.00160.0014-0.0163-0.0142-0.0124-0.0110-0.00980.08970.08360.07830.07370.0695-0.00020.0012-0.00870.0658ezplot('(sqrt(1+x)-sqrt(x-1)./(sqrt(2+x)+sqrt(x-1)',3,8)%在區(qū)間繪制 yhold on%保持figure不關(guān)閉fnplt(Sp)%繪制出插值效果(sqrt(1 -+x)-£qrt(K-1 )/(sqrl(2*x)+sqrt(M-1)分段三次樣條插值擬合效果