【正文】 Step”)*2)_ 1OS ?I ,5 2 ■ o x = 0::10。y = 2*x+3。subplot(411)。plot(x,y)。grid。title(“y 的原函數(shù)”)。subplot(412)。semilogy(x,y)。grid。title(“對 y 取對數(shù)”)。丫的原畫數(shù) 401111111l| I p | il ■ | i| I 九 _____ 1___ — ___ I ____ L ___ :…… ：___ J ________ L ___ u i| I |l I , il _■」 “ j I ■I __ h_____________ I I ■ Q 」 【 I 1F I I I I II I I IT 10 1□ 1 2 3 4 5 6 r 6 9 10 , 對 y 取對數(shù) 對弋觀對數(shù) subplot(413)。semilogx(x,y)。10 10 id 10 1lZ 10 w 10 10 40 20 0 對好對數(shù) ?掌握循環(huán)、分支語句的編寫，學會使用 look for、help 命令grid。title(”對 x 取對數(shù)“)。subplot(414)。loglog(x,y)。grid。title(”對 xy 均取對數(shù)“)。x =3::3。bar(x,exp(x.*x),”g“)實驗三 MATLAB 程序設(shè)計 實驗任務(wù) 1 ?熟悉 MATLAB 程序設(shè)計的方法和思路；程序舉例 f=[1,1]。i=1。while f(i)+f(i+1) m=3。n=4。for i=1:m for j=1:n a(i,j)=1/(i+j1)。end end format rat1/2 1/3 1/2 1/3 1/4 1/3 1/4 1/5(分數(shù)格式形式。用有理數(shù)逼近顯示數(shù)據(jù))m=5。n=4。for i=1:m for j=1:n a(i,j)=1/(i+j1)。end end format rat a a = 1 1/2 1/3 1/4 1/2 1/3 1/4 1/5 1/3 1/4 1/5 1/6 1/4 1/5 1/6 1/7 1/5 1/6 1/7 1/8 程序中沒有 format rat 命令時，如果上次運行結(jié)果沒有清除，輸出的結(jié)果就是上次運行的 結(jié)果！但是運用 clear 命令清楚之前的運行結(jié)果之后就會正常運行。 x=input（”請輸入 x 的值:’）。if x==10 y=cos(x+1)+sqrt(x*x+1)。else y=x*sqrt(x+sqrt(x))。1/4 1/5 1/6end y請輸入 x 的值 :2 y = 2391/647 x=input(“請輸入 x 的值:’)。if x==10 y= fprintf(” 不在定義域內(nèi)，請重新輸入：“)。return else y=1/(x10)。end y 請輸入 x 的值 :21/8 p=[0 0 0 1 3 0 2 0 0 9]。for i=1:length(p), if p(1)==0,p=p(2:length(p))。end。end。pp = Columns 1 through 5 300 Columns 6 through 70 9p=[0 0 0 1 3 0 2 0 0 9]。p(p==0)=[]。pp =3 e2(500)ans =1 2 3 5 8 13 21 34 55 89 144 233377 lookfor ffibno e2ffibno 計算斐波那契亞數(shù)列的函數(shù)文件 help e2 ffibno 計算斐波那契亞數(shù)列的函數(shù)文件 n 可取任意自然數(shù) 程序如下（用法：lookfor 關(guān)鍵詞 在所有 M 文件中找 “關(guān)鍵詞 ”，比如：lookfor max（即尋找關(guān)鍵詞 “ max”）其實就和我們平時用 CTRL+F 來查找關(guān)鍵詞”是一樣的 而 help 是顯示 matlab 內(nèi)置的幫助信息 用法：help 命令，比如 help inv，作用就是調(diào)用 inv 這個命令的幫助）程序設(shè)計題 用一個 MATLAB 語言編寫一個程序：輸入一個自然數(shù)，判斷它是否是素數(shù)，如果是，輸出“ It is one prime ”，如果不是，輸出“ It is n ot o ne prime.”。要求通 過調(diào)用子函數(shù)實現(xiàn)。最好能具有如下功能：①設(shè)計較好的人機對話界面，程序中 含有提示性的輸入輸出語句。②能實現(xiàn)循環(huán)操作，由操作者輸入相關(guān)命令來控制 是否繼續(xù)進行素數(shù)的判斷。如果操作者希望停止這種判斷，則可以退出程序。③ 如果所輸入的自然數(shù)是一個合數(shù)，除了給出其不是素數(shù)的結(jié)論外，還應(yīng)給出至少 一種其因數(shù)分解形式。例：輸入 6，因為 6 不是素數(shù)。則程序中除了有“ It is not one prime ”的結(jié)論外，還應(yīng)有：“ 6=2*3 ”的說明。function sushu while 1 x=input（” 請輸入一個自然數(shù) “）。if x“)for n=2:sqrt(x)if rem(x,n)==0 num3=x。num1=n。num2=x/n。disp([num2str(num3),”=“,num2str(num1),”x“,num2str(num2)])End End Endend y=input(” 是否繼續(xù)判斷？繼續(xù)請按 1 ,按任意鍵退出 ：“)； if y~=1 break end end 實驗四 MATLAB 的符號計算與 SIMULINK 的使用 實驗任務(wù) MATLAB 符號計算的特點和常用基本命令。 SIMULINK 的使用。程序舉例 a=sym(”[a11 a12。a21 a22]“)。da=det(a)ea=eig(a)da = a11*a22a12*a21 ea = 1/2*a11+1/2*a22+1/2*(a11 A 22*a11*a22+a22 A 2+4*a12*a21)A(1 /2)1/2*a11+1/2*a221/2*(a11 A 22*a11*a22+a22 A 2+4*a12*a21)A(1 /2)a=sym(”[2 3。1 5]“)。da=det(a)ea=eig(a)da = 7 ea = 7/2+1/2*21A(1/2)7/21/2*21A(1/2)(包括精確解和一定精度的解) r1=solve(”xA2+x1“)rv=vpa(r1)rv4=vpa(r1,4)rv20=vpa(r1,20)r1 = 1/2*5 八(1/2)1/21/2*5 八(1/2)1/2 rv =. rv4 =. rv20 =. 3 ． a=sym(”a“)。b=sym(”b“)。c=sym(w=10。x=5。y=8。z=11。A=[a,b。c,d] B=[w,x。y,z] det(A)det(B)A = [ a, b] ”c“)。d=sym(”d“)。%定義 4 個符號變量 %定義 4 個數(shù)值變量 %建立符號矩陣 A %建立數(shù)值矩陣 B %計算符號矩陣 A 的行列式 %計算數(shù)值矩陣 B 的行列式[ c, d]X108 11 ans = a*db*c ans = 150 4.187。 syms x y。s=(7*x A 28*y A 2)*(x A 2+3*y A 2)。expand(s)% 對 s 展開 collect(s,x)% 對 s 按變量 x 合并同類項（無同類項）factor(ans)% 對 ans 分解因式ans = 7*xA413*xA2*yA224*yA4 ans = 7*xA413*xA2*yA224*yA4 ans =(8*yA2+7*xA2)*(xA23*yA2) AX=b 求解 A=[3