MATLAB 試驗 符號運算

本文格式為Word版，下載可任意編輯——MATLAB試驗符號運算寧德師范學院物理與電氣工程系

試驗報告

（2023—2023學年其次學期）課程名稱MATLAB試驗名稱試驗一符號運算專業(yè)電氣工程及其自動化年級11電氣2班學號姓名指導教師試驗日期2023年5月13日

試驗一符號運算

一、試驗目的

1、了解符號、數(shù)值、字符等數(shù)據(jù)類型的區(qū)別2、了解符號運算的特點、優(yōu)缺點

3、把握符號變量的創(chuàng)立和運算，以及其運算的基本應用4、把握基本的符號繪圖指令二、試驗內容

1、指出下面的M1，M2，M3分別是什么，并上機驗證。

取a=1、b=2、c=3、d=4，M1=[a,b;c,d]，M2='[a,b;c,d]'M3=sym('[a,b;c,d]');

程序：

>>a=1,b=2,c=3,d=4a=1b=2c=3d=4

>>M1=[a,b;c,d]

，M1=1234>>M2='[a,b;c,d]'M2=[a,b;c,d]

>>M3=sym('[a,b;c,d]')M3=[a,b][c,d]

結論：M1是矩陣，M2是向量，M3是數(shù)組。

a?b?a2、下面2種取值狀況下，計算b并賦給相應狀況下的c1、c2，c1、c2相等嗎，為什么？上機驗證。

（1）a1=10^10;b1=10^(-10);

（2）將a1、a2作為符號變量賦給a2、b2;程序：

>>a1=1e10;b1=1e-10;c1=(a1+b1-a1)/b1a2=sym(a1);b2=sym(b1);c2=(a2+b2-a2)/b2c1=0

問

c2=1

由于c1=0，c2=1，所以c1不等于c2。a1、b1是具體的數(shù)值，a2、b2是符號變量。

3、符號表達式中自由變量的確定生成符號變量a、b、x、X、Y、k=3、

2f?azX?(bx?k)Y。cw?ysinaz=，表達式為

（1）找出f中的全部自由符號變量（2）在f中確定最優(yōu)先的自由符號變量

（3）在f中確定2個和3個自由變量時的執(zhí)行狀況

（4）試通過對各符號變量與x的ASCII值做絕對差值，分析自由變量優(yōu)

秀順序，能得出什么結論？程序：

>>symsabxXYk=sym('a');

z=sym('c*sqrt(w)+y*sin(a)');f=a*z*X+(b*x^2+k)*Y;

>>findsym(f)%找出f中的全部自由符號變量ans=

X,Y,a,b,c,w,x,y

>>findsym(f,1)%在f中確定最優(yōu)的自由符號變量ans=

x

>>findsym(f,2)%在f中確定2個自由變量時的執(zhí)行狀況ans=

x,y

>>findsym(f,3)%在f中確定3個自由變量時的執(zhí)行狀況ans=

x,y,w>>double('a')ans=97>>double('b')ans=98

>>findsym(f,8)ans=

x,y,w,c,b,a,X,Y>>A=abs(97-120)A=23