Matlab7.0符號計算上機實驗1.doc

Matlab7.0符號計算上機實驗1.doc

ID:62033742

大?。?21.50 KB

頁數(shù):25頁

時間:2021-04-15

Matlab7.0符號計算上機實驗1.doc_第1頁
Matlab7.0符號計算上機實驗1.doc_第2頁
Matlab7.0符號計算上機實驗1.doc_第3頁
Matlab7.0符號計算上機實驗1.doc_第4頁
Matlab7.0符號計算上機實驗1.doc_第5頁
資源描述:

《Matlab7.0符號計算上機實驗1.doc》由會員上傳分享,免費在線閱讀,更多相關(guān)內(nèi)容在教育資源-天天文庫。

1、Matlab7。0符號計算實驗報告⒈目的本實驗旨在向?qū)W生介紹一種解決專業(yè)問題的快速有效且具有強大功能的科學(xué)與工程計算軟件。通過本實驗,應(yīng)使學(xué)生掌握的內(nèi)容是:MATLAB7。0符號對象的創(chuàng)建和使用,MATLAB7。0任意精度的計算, MATLAB7.0符號表達式的化簡和替換,MATLAB7.0符號矩陣計算,MATLAB7。0符號微積分,MATLAB7.0積分變換,MATLAB7。0符號代數(shù)方程和符號微分方程求解.該實驗主要為上機實驗,要求學(xué)生按要求上機實現(xiàn)相關(guān)的程序的設(shè)計,自己動手編寫程序并驗證程序的正確性.⒉實驗任務(wù)

2、分解通過一些實例初步掌握MATLAB7.0的基本符號計算功能和應(yīng)用。實驗任務(wù)可分解為:MATLAB7。0符號計算的基礎(chǔ),MATLAB7.0符號計算在高等數(shù)學(xué)中的應(yīng)用,MATLAB7.0符號方程(組)求解和符號矩陣計算。⒊實驗環(huán)境介紹長清校區(qū)數(shù)學(xué)實驗室實驗一MATLAB7。0符號計算的基礎(chǔ)實驗?zāi)康?1.掌握MATLAB7。0符號對象的創(chuàng)建和使用;2.掌握MATLAB7.0任意精度的計算;3.掌握MATLAB7.0符號表達式的化簡和替換。實驗要求:給出程序和實驗結(jié)果。實驗題目與結(jié)果:一、計算符號表達式在處的值,并將結(jié)果設(shè)

3、置為以下五種精度,即分別為小數(shù)點之后1位、2位、10位、20位、50位有效數(shù)字;>>x=pi—1;>〉y=x+cos(x)-sin(x)y = 0.7598>>y1=vpa(y,1)y1=.8>〉y2=vpa(y,2)y2=.76>〉y3=vpa(y,10)y3=.7598193629〉>y4=vpa(y,20)y4 =。75981936291375673509>〉y5=vpa(y,50)y5=.75981936291375673508952104384661652147769927978516二、設(shè)為符號變量,,,

4、試進行如下運算。(1),(2),(3)對進行因式分解,(4)求的反函數(shù)。>>symsxfg;>〉f=x^4+2*x^2+1;>>g=x^3+6*x^2+3*x+5;(1)>>a=f+ga?。健^4+8*x^2+6+x^3+3*x (2)〉>b=f*gb=(x^4+2*x^2+1)*(x^3+6*x^2+3*x+5)>> expand(b)ans =x^7+6*x^6+5*x^5+17*x^4+7*x^3+16*x^2+3*x+5(3)〉>factor(f)ans =(x^2+1)^2(4)〉>finverse(g)

5、Warning:finverse(x^3+6*x^2+3*x+5)isnot unique.>In sym。finverseat43ans= 1/2*(-60+4*x+4*(117—30*x+x^2)^(1/2))^(1/3)+6/(—60+4*x+4*(117-30*x+x^2)^(1/2))^(1/3)-2三、指出下面各條語句的輸出結(jié)果>〉f=sym(’2*u’);〉>subs(f,'u',2)>〉f2=subs(f,’u','u+2’)〉>a=3;〉> subs(f2,'u',a+2)>>subs(f2,’u’

6、,'a+2')>〉syms xy>〉f3=subs(f,’u’,x+y)>>subs(f3,[x,y],[1,2])〉〉subs(f3,[x,y],[x+y,x+y])>>f=sym('2*u');>>subs(f,'u’,2)ans=4〉〉 f2=subs(f,'u','u+2')f2=2*(u+2)?。荆尽=3;>> subs(f2,’u',a+2)ans = 14〉> subs(f2,’u','a+2')ans=2*((a+2)+2) 〉>symsxy>〉f3=subs(f,'u’,x+y)f3 =2*x+2

7、*y>>subs(f3,[x,y],[1,2])ans=  6〉〉subs(f3,[x,y],[x+y,x+y])ans=6*x+4*y四、簡化>>symsx;>〉f=((1/x^3)+(6/x^2)+12/x+8)^(1/3);〉>simplify(f)ans = ((2*x+1)^3/x^3)^(1/3)五、設(shè)a、b定義如下,試上機輸出factor(a)和factor(b)的結(jié)果,并指出哪個結(jié)果才是12345678901234567890的因式分解,為什么? a=sym('12345678901234567890

8、’) b=sym(12345678901234567890)y1=factor(a)?。?=factor(b)y1=?。?)*(3)^2*(5)*(101)*(3803)*(3607)*(27961)*(3541)y2=(2)^11*(7)*(7324703)*(117570121)?。尽怠=sym('12345678901234567890

當前文檔最多預(yù)覽五頁,下載文檔查看全文

此文檔下載收益歸作者所有

當前文檔最多預(yù)覽五頁,下載文檔查看全文
溫馨提示:
1. 部分包含數(shù)學(xué)公式或PPT動畫的文件,查看預(yù)覽時可能會顯示錯亂或異常,文件下載后無此問題,請放心下載。
2. 本文檔由用戶上傳,版權(quán)歸屬用戶,天天文庫負責(zé)整理代發(fā)布。如果您對本文檔版權(quán)有爭議請及時聯(lián)系客服。
3. 下載前請仔細閱讀文檔內(nèi)容,確認文檔內(nèi)容符合您的需求后進行下載,若出現(xiàn)內(nèi)容與標題不符可向本站投訴處理。
4. 下載文檔時可能由于網(wǎng)絡(luò)波動等原因無法下載或下載錯誤,付費完成后未能成功下載的用戶請聯(lián)系客服處理。