精通matlab科學(xué)計(jì)算

第2 的擴(kuò)展運(yùn)算，它還支持復(fù)數(shù)這一數(shù)值元素，這也是區(qū)別于其他高級(jí)語(yǔ)言的最大數(shù)值類數(shù)值類4變量是數(shù)值計(jì)算的基本單元。與C語(yǔ)言等其他高級(jí)語(yǔ)言不同，語(yǔ)言中的變量中不允許有未定義的變量，也不需要預(yù)先定義變量的類型，會(huì)自動(dòng)生成變量，1）變量命名規(guī)則中的變量命名規(guī)則如下：變量名區(qū)分大小寫，因此A與a某些常量也可以作為變量使用，如i在中表示虛數(shù)單位，但也可以作常量是指那些在中已預(yù)先定義其數(shù)值的變量，默認(rèn)的常量如表2.1所示表2.1默認(rèn)常 圓周浮點(diǎn)數(shù)的相對(duì)誤INF（或無(wú)窮i（或虛數(shù)單位，定義為NaN（或代表不定值（即函數(shù)實(shí)際輸入?yún)?shù)個(gè)最大的正實(shí)函數(shù)實(shí)際輸出參數(shù)個(gè)最小的正實(shí)ANS（或默認(rèn)變量名，以應(yīng)答最近一次操作運(yùn)算 指定變量時(shí)，賦值給一個(gè)特殊的變量ans，數(shù)據(jù)的顯示格式由format命令控制。format只表2.2的數(shù)據(jù)顯示格 format短格式（5位定點(diǎn)數(shù)formatlong長(zhǎng)格式e方format長(zhǎng)格式（15位定點(diǎn)數(shù)format2位十進(jìn)制格formatshort短格式e方format十六進(jìn)制格工作空間中的變量可以用save命令到磁盤文件中。鍵入命令“save<文件名 用load命令可將變量從磁盤文件讀入 字符是 中符號(hào)運(yùn)算的基本元素，也是文字等表達(dá)方式的基本元素，在中字符串作為字?jǐn)?shù)組用單引號(hào)（）到程序中，還可以過字符串運(yùn)算元胞是元胞數(shù)組（CellArray）的基本組成部分。元胞數(shù)組與數(shù)字?jǐn)?shù)組相似，以下標(biāo)來(lái)區(qū)分，單元胞數(shù)組由元胞和元胞內(nèi)容兩部分組成。用花括號(hào){}表示元胞數(shù)組的內(nèi)容，用圓2-1】解 A(1,1){'AnexampleofcellA(1,2){[12;34]};A{2,1}tf(1,[1,8]);A{2,2}{A(1,2);'Thisiscelldisp(A)%顯示該元胞數(shù)元胞數(shù)組A12行用元胞內(nèi)容串，最后，用celldisp函數(shù)顯示該元胞數(shù)組A。與元胞數(shù)組相似，結(jié)構(gòu)數(shù)組（StructureArray）也能存放各類數(shù)據(jù)，使用指針方式傳遞例如，可用parameter.temperature表示某一對(duì)象的溫度參數(shù)，用parameter.humidity表示某一對(duì)象的濕度參數(shù)等，因此，該結(jié)構(gòu)數(shù)組parameter由兩個(gè)屬性組成。面象的語(yǔ)言采用了多種對(duì)象，如自動(dòng)控制中常用的傳遞函數(shù)模型對(duì)objectobject2.2關(guān)系運(yùn)算和邏輯運(yùn)，除了傳統(tǒng)的數(shù)算外還支持關(guān)系運(yùn)算和邏輯運(yùn)算如果你已經(jīng)有了一些關(guān)系運(yùn)算和邏輯運(yùn)算主要用于控制基于真/假命題的各命令（通常在M文件中），2.32.3關(guān)系運(yùn)算 <小大于等小于等等大不等2.4邏輯運(yùn)算 &邏輯~邏輯|邏輯此外，還提供了若干關(guān)系運(yùn)算函數(shù)和邏輯運(yùn)算函數(shù)，分別如表2.5和表2.5關(guān)系運(yùn)算函 所有向量為非零元素時(shí)為邏輯異或運(yùn)任一向量為非零元素時(shí)為2.6邏輯運(yùn)算函 位方式的邏輯與運(yùn)位比較運(yùn)位方式的邏輯或運(yùn)最大無(wú)符號(hào)浮點(diǎn)整位方式的邏輯異或運(yùn)將二進(jìn)制移位運(yùn)2.3矩陣及其運(yùn) 都是以矩陣為單元進(jìn)行的，可見矩陣是的。表2.7列出提供的每mn列構(gòu)成的數(shù)組amn階矩陣，它總共由mn個(gè)元素組成，矩陣元素記為aij，其中i表示行，j表示列。mn時(shí)，矩陣a稱為方陣。當(dāng)ij時(shí)，所有的aij0，且mn，得到的矩陣稱為對(duì)于(mnwwijaji時(shí)，稱wa的轉(zhuǎn)置矩陣，記為wa對(duì)于a為(m1的形式時(shí)，稱am個(gè)元素的列向量，對(duì)于a為(1n的形式時(shí)，稱n“;）2-2】解 a[123;452×3a，a11、2、33個(gè)元素組成，第4、5、63a23456b[a111213]%添加一行元素[11123×3b，b矩陣是在a11、12、13，組成一個(gè)3×3矩陣，輸出結(jié)果如下：b 單個(gè)元素的：b(3,2)→12，了第3行和第2列交叉的元素整列元素的：b(:,3)→[3,6,13]，了第3列中的所有元素整行元素的：b(1,:)→[1,2,3]，了第1行中的所有元素整塊元素的：b(2:3,2:3)→[5,6;12,13]，了一個(gè)（2×2）的子塊矩陣表2.7常用特殊矩陣生成函 功能說 功能說zeros(生成元素全為0的矩tril(生成下三角ones(生成元素全為1的矩eye(生成單位矩rand(生成均勻分布隨機(jī)矩company(生成伴隨矩randn(生成正態(tài)分布隨機(jī)矩hilb(生成Hilbert矩magic(生成魔方矩vander(生成vander矩diag(生成對(duì)角矩hankel(生成hankel矩triu(生成上三角hadamard(生成hadamard矩2-3】解 a[12,34,56;78,9];btril(a)%生成下三角矩btril生成的a矩陣的下三角b004507892.8矩陣基本運(yùn)操作符功能說操作符功能說矩陣加/矩陣的左矩陣減‘矩陣轉(zhuǎn)*矩陣乘ogm(矩陣對(duì)數(shù)運(yùn)^矩陣的expm(矩陣指數(shù)運(yùn)\矩陣的右inv(矩陣求2-4】解 a[12;34];b[35;29];div1a/b;%div2b\a%abdiv12.9常用矩陣函數(shù)功能說功能說rot90(矩陣逆時(shí)針旋轉(zhuǎn)eig(計(jì)算矩陣的特征值和特征向flipud(矩陣上下翻rank(計(jì)算矩陣的fliplr(矩陣左右翻trace(計(jì)算矩陣的flipdim(矩陣的某維元素翻norm(計(jì)算矩陣的shiftdim(矩陣的元素poly(計(jì)算矩陣的特征方程2-5】解 a[1352467913];[bc]eig(a)%求取矩陣的特征值和特征通過函數(shù)eig(ab和特征值cb0.3008 c0000002.10常用矩陣分解運(yùn)算功能說功能說eig(矩陣的特征值分svd(矩陣的奇異值分qr(矩陣的QR分chol(矩陣的Cholesky分schur(矩陣的Schur分lu(矩陣的LU分2-6】解：程序代碼如下a[62,12,31;11,1];[L,UP]lu(a)%對(duì)矩陣進(jìn)行LU分lu(aLUULP，L00U0000P 0 0 1中的數(shù)據(jù)中的數(shù)據(jù)精 的數(shù)據(jù)類型表2.11整數(shù)的表示范數(shù)據(jù)類表示范數(shù)據(jù)類表示范-2727-028--215215-0216--231231-0232--263263-0264-當(dāng)數(shù)據(jù)超過類型的表示范圍時(shí)，將數(shù)據(jù)表示成該類型的最大值或最小值，2-7】234與-234int8解：在命令窗口輸入下列命令int8(234)ans int8(-234)ans-128由于int8表示的整數(shù)范圍為- 27-1，而234127，所以輸出結(jié)果為int8能夠127，同理由于-234-128int8能夠表示的最大整數(shù)-128。范圍如表2.12所示。表2.12浮點(diǎn)數(shù)數(shù)的表示范浮點(diǎn)類表示范浮點(diǎn)類表示范單精-2128-2-126,2-126雙精-21024-2-1022,2-1022表2.12的數(shù)據(jù)來(lái)自于的幫助文件，這些數(shù)據(jù)都是遵循IEEE標(biāo)準(zhǔn)得出的數(shù)據(jù)。從表2.12可以看出，雙精度數(shù)據(jù)能夠表示的最小實(shí)數(shù)為2-1022，即精度表示數(shù)據(jù)的精度為2-1022。當(dāng)然， 2-8】解：在命令窗口輸入下列命令str'雙精度數(shù)據(jù)的表示范圍為:\n\t%gtogn\t%gto%sprintf(str,-realmax,-realmin,realmin,realmax)ans-1.79769e308to-2.22507e-3082.22507e-308to1.79769e308 的數(shù)值精度的數(shù)值精度也就是能夠表示的最小實(shí)數(shù)，任何一個(gè)絕對(duì)值小于的數(shù)值精度的實(shí)數(shù)都會(huì)被當(dāng)成0處理。經(jīng)過實(shí)驗(yàn)，可以知道，在7版本中 能夠表示的最小實(shí)數(shù)為2-1074，任何絕對(duì)值小于2-1074的實(shí)數(shù)0【例2-9】數(shù)值精度實(shí)例。7的數(shù)值精度解：在命令窗口輸入下列命x2^(1074)x4.9407e324x2^(1075)x x0ans 從上面的例子可以看出，由于2-10752-1074，7視其為0，而且通過和0的因此，如果在程序中要判斷某個(gè)實(shí)數(shù)是否等于0，最可靠的辦法是將其2-1074相比較，看它的絕對(duì)值是否小于等于2-1074當(dāng)然對(duì)于具體的問題可以根據(jù)問題的需要和精度，采取不同的比較對(duì)象，例如中有一個(gè)內(nèi)置常量eps，其值為2.220410-16，它可以作為實(shí)數(shù)是否等于0的一2.4.3的顯示精的顯示精度是指顯示的有效位數(shù)。中的顯示精度是可以修改的，顯示精度修改了，原來(lái)的數(shù)據(jù)并沒有變，只是數(shù)據(jù)在命令窗口中顯體用法可參考幫助文檔?！纠?-10】顯示精度實(shí)例。的顯示精度解：在命令窗口輸入下列命令x %默認(rèn)的顯示精x format %設(shè)置為long顯示精xx format x1/3x %顯示10位有效數(shù)ans %顯示20位有效數(shù)ans2.5符號(hào)oolbox符號(hào)表達(dá)式是代表數(shù)字、函數(shù)、算子和變量的字符串，或字符串?dāng)?shù)組。不符號(hào)表達(dá)式幾乎完全像基本的命令。和symsyms函數(shù)sym一次只能定義一個(gè)符號(hào)變量使用不方便提供了另一個(gè)函數(shù)syms，syms符號(hào)變量名1符號(hào)變量名2符號(hào)變量名sym2-11】解：在窗口，輸入下列命令 symsxy; %建立符號(hào)變量x、y %定義符號(hào)表達(dá)式 例如，在命令窗口中輸入Asym('[a,2*b;3*a,0]')A

2b0 A[a,[3*a,需要注意的是：符號(hào)矩陣的每一行的兩端都有方括號(hào)，這是與數(shù)值矩陣的。 還有。transpose(s)：返回sdeterm(s)sfactor(S)：對(duì)S分解因式，Sexpand(S)S進(jìn)行展開，Scollect(S,v)Sv合并同類項(xiàng)，S diff是求微分最常用的函數(shù)，其輸入?yún)?shù)既可以是函數(shù)表達(dá)式，也可以是符號(hào)矩陣。常用的格式是：diff(f,x,n)，表示f關(guān)于x求n階導(dǎo)數(shù)。常用的格式是：int(frx0,x1)。其中，f為所要積分的表達(dá)式，r為積分變量，若為定積分，則x0與x1為積分上下限。2-12】1fsinax)xa symsa %定義符號(hào)變量a和 %創(chuàng)建函數(shù)dfxdiff(f %對(duì)x求dfadiff(f %對(duì)a求fdfx %f對(duì)x求導(dǎo)的結(jié)dfa %f對(duì)a求導(dǎo)的結(jié)2-13】2。已知表達(dá)式fxlog(1xxx在(0,1)解：輸入如下程序代碼syms %定義符號(hào)變量 %創(chuàng)建函數(shù) %對(duì)x積 int11/2*(1x)^2*log(1x)3/41/2*x-1/4*x^2-(1x)*log(1x)%積分表達(dá)式int21/4 2.6復(fù)數(shù)及其運(yùn)ij字元來(lái)代表虛部復(fù)數(shù)運(yùn)算的。xabia稱為實(shí)部，bxrei。其中r稱為復(fù)數(shù)的模，又記為|x|；稱為復(fù)數(shù)的幅角，又記為Arg(x)，且滿足如下關(guān)系：a2r ，tana2直接法就是利用符號(hào)i或j來(lái)表示復(fù)數(shù)單位，將復(fù)數(shù)看作完整的一個(gè)表達(dá)式輸入。其所謂符號(hào)函數(shù)法就是將復(fù)數(shù)看成是函數(shù)形式，其實(shí)部和虛部看作自變量，用syms來(lái)構(gòu)造，用subs對(duì)符號(hào)函數(shù)中自變量進(jìn)行賦值?！纠?-14】復(fù)數(shù)構(gòu)造實(shí)例。試分別使用上述兩種方法，在中構(gòu)造復(fù)x1i解：在命令窗口中輸入 symsab %a，b為實(shí)數(shù) symsrct %r，ct為實(shí)數(shù) x11.00001.0000ix21.00001.0000ix31.00001.0000ix41.0000【例 】復(fù)數(shù)矩陣構(gòu)造實(shí)例。構(gòu)造復(fù)數(shù)矩

13i11 在命令窗口中輸入：A1[sqrt(2)*exp((pi/4)*i)12i13i;sqrt(2)*exp((-pi/4)*i)1-2i1-3i]A2re[111;111];A2im[123;-1-2-3];輸出結(jié)果A1和A21.00001.0000i1.00002.0000i1.00001.00001.0000i1.00002.0000i1.0000提供了繪制極坐標(biāo)圖的函數(shù)polar，它還可以繪制出極坐標(biāo)柵格線，其常用其中，theta為極坐標(biāo)極角，rho2-16】ytitsin(t 命令窗口中輸入以下命令，最終顯示結(jié)果如圖2-1所示。圖-1極坐標(biāo)下復(fù)數(shù)的t0:0.01:2*pi;yti*t.*sin(t); title('直角坐標(biāo)圖 對(duì)于一個(gè)形如xabi的復(fù)數(shù)通常希望能夠了解它自身的結(jié)構(gòu)性質(zhì)包括實(shí)部、 提供了方便的操作函數(shù)，如表2.13所示。下面利用實(shí)例來(lái)2.13常用矩陣分解運(yùn)算 求復(fù)數(shù)或復(fù)數(shù)矩陣A的實(shí)求復(fù)數(shù)或復(fù)數(shù)矩陣A的求復(fù)數(shù)或復(fù)數(shù)矩陣A的虛求復(fù)數(shù)或復(fù)數(shù)矩陣A的相角，單位為弧求復(fù)數(shù)或復(fù)數(shù)矩陣A的共 a是f(z的孤立奇點(diǎn)，Caa的閉路，積分

Cf(z)dz稱為f(z)a點(diǎn)的留數(shù)或殘數(shù)，記作Resf(za留數(shù)定理如果函數(shù)f(z)在閉路C上解析在C的除去n個(gè)孤立奇點(diǎn)a1,a2nnCf(z)dz2iRes[f(z),akk如果能夠求得函數(shù)f(z)在各極點(diǎn)的留數(shù)的顯式表達(dá)，則該閉路積分很容易獲得。由羅朗展開：若a是f(z)的m重極點(diǎn)，則函數(shù)在該點(diǎn)的留數(shù)可以表示為： Res[f(z),a] (m1)!zadz

f由于在工程中遇到的f(zanznan1zn1...a1za0bmzmbm1zm1...bmzb0residue可以求得該有理式的留數(shù)，residue的返回參數(shù)有3個(gè)，分別對(duì)應(yīng)留數(shù)向調(diào)用格式：[rp,k]residue

a0]

z2-17】留數(shù)計(jì)算實(shí)例。求如下函數(shù)的奇點(diǎn)處的留數(shù)z22z解：在命令窗口中輸入r=p20kz因此，可得Resf(z21.5;Resf(z00.5z【例2-18】留數(shù)計(jì)算實(shí)例2|z|2

Cz4

dz，其中C解：在命令窗口中輸入r=-0.2500--0.2500+p0.0000+0.0000-k因此，可知在圓周

|z|

內(nèi)有四個(gè)極點(diǎn)，因此積分值等于2pi0.250.250.250.2502.7復(fù)變函數(shù)及其運(yùn)Taylor級(jí)數(shù)展開在復(fù)變函數(shù)中有很重要的地位，如分析復(fù)變函數(shù)的解析性等。變換、拉斯變換與z變換，是理工科中的三大積分變換。Z變換可以說是針對(duì)離散信號(hào)和系統(tǒng)的拉斯變換，由此就很容易理解Z變換的重要性，也很容易理解Z變換和變換之間的關(guān)系。在數(shù)字信號(hào)處理中，Z變換是Taylor展開、積分變換等復(fù)變函數(shù)運(yùn)算提供了支持，f(x)x=x01(階)直到(n+1)階的導(dǎo)數(shù)，則在該鄰域內(nèi)，函數(shù)f(x)在點(diǎn)x=x0時(shí)，項(xiàng)數(shù)趨向無(wú)窮的冪級(jí)數(shù)如下： f f(x)

f(x)f‘() 0 0 這個(gè)冪級(jí)數(shù)叫做函數(shù)f(x)的Taylor（泰勒）級(jí)數(shù)，在中可由函數(shù)taylor來(lái)taylor(f)：返回f函數(shù)的五次冪多項(xiàng)式近似，函數(shù)ffindsym確定。taylor(f,n)：返回n-1次冪多項(xiàng)式，函數(shù)f的自變量由函數(shù)findsym()確定。taylor(f,a)：返回a點(diǎn)附近的五次冪多項(xiàng)式近似。taylor(f,x)ffx，不是由函數(shù)findsym()確定。2-19】泰勒級(jí)數(shù)展開實(shí)例。求下列函數(shù)的泰勒展式：f1(x)=ln(1+x)、f2(x)=(1+x)mf3(x)=sin(x)在處的泰4解：在命令窗口中輸入symsxm;f2=taylor((1+x)^m);f3=taylor(sin(x),pi/4,f1=x^5/5-x^4/4+x^3/3-x^2/2+f2=(m/5-(m*(m/4-m^2/6))/3+(m*(m*(m/12-m^2/24)-m/6+m^2/12))/2-m*(m/8-(m*(m/12-m^2/24))/2+m*(m*(m/48-m^2/120)-m/18+m^2/48)-m^2/18)-m^2/8)*x^5+((m*(m/4-m^2/6))/2-m/4-m*(m*(m/12-m^2/24)-m/6+m^2/12)+m^2/6)*x^4+(m/3-m*(m/4-m^2/6)-m^2/4)*x^3+(m^2/2-m/2)*x^2+m*x+1f3=(2^(1/2)*(pi/4-x)^3)/12-(2^(1/2)*(pi/4-x)^2)/4+2^(1/2)/2-(2^(1/2)*(pi/4-55 x

2x2

1 32sin(x) 1 2

4

4

4Fourierifourier1、fourierF=fourier(f)x為數(shù)量符號(hào)fFourierw的函f=f(w)fourier函數(shù)返回t的函數(shù)F=F(t)F(w)=int(f(x)*exp(-i*w*x),x,-inf,inf)為對(duì)x的積分。F=fourier(f,v)v代替默認(rèn)變量wFourierfourier(f,v)F(v)=int(f(x)*fourier(f,u,v)vxu積分，且有fourier(f,u,v)等價(jià)于F(v)=int(f(u)*exp(-i*v*u),u,-inf,inf)。2、ifourierf=ifourier(F)：返回默認(rèn)獨(dú)立變量wFFourierx的函數(shù)。如果F=F(x)，則Ifourier函數(shù)返回t的函數(shù)f=f(t)。一般來(lái)說f(x)=1/(2*pi)*int(F(w)*exp(i*w*x),w,-inf,inf)，對(duì)w積分。f=ifourier(F,u)：以變量u代替x，且ifourier(F,u)等價(jià)于f(u)=1/(2*pi)*int(F(w)*exp(i*w*u,w,-inf,inf)。對(duì)w積分。=pi)*int(F(v)*exp(i*v*u,v,-inf,inf)v。解：在命令窗口中輸入clearall;symstbwF2=fourier(diff(sym('f(t)')),t,w)F1=pi*(2*heaviside(-w)-1)*iF2=w*transform::fourier(f(t),t,-w)*iF3=heaviside(b+w)-heaviside(w-b)

t

0,其

1, t

0,其【例2-21】逆變換實(shí)例。求下列函數(shù)的逆變換：F1()

b解：在命令窗口中輸入clearall;symstbwF2=1/(b+j*w);f1=heaviside(t+1)/2-heaviside(t-f2=((2*pi*heaviside(t))/exp(b*t)-pi/exp(b*t)+ 2

t0 tLaplace在中可由函數(shù)laplace來(lái)實(shí)現(xiàn)拉斯變換，ilaplace函數(shù)來(lái)實(shí)現(xiàn)拉斯1、laplace斯變換。函數(shù)返回斯變換。函數(shù)返回默認(rèn)為S的數(shù)。如果F=F(s)，則 斯函數(shù)返回t的函數(shù)L=L(t)。其中定義L為對(duì)t的積分L=laplace(F,t)：以t代替s為變量的拉int(F(x)*exp(-t*x)0,inf)斯變換。laplace(F,t等價(jià)于 于L(z)int(F(w)*exp(-z*w),0,inf)其中，F(xiàn)是時(shí)域函數(shù)表達(dá)式，約定的自變量是t,得到的 斯變換函數(shù)是L(s)2、ilaplace斯變換，默認(rèn)返回斯變換，默認(rèn)返回tF=ilaplace(L)：返回以默認(rèn)獨(dú)立變量s的符號(hào)表達(dá)式L如果L=L(t)，則ilaplace返回x的函數(shù)F=F(xF(xs都在直線scF=ilaplace(L,y)y代替默都在直線scF=ilaplace(L,y)y代替默認(rèn)的t的函數(shù)且有ilaplace(L,y)等價(jià)于F(y)=int(L(y)*exp(s*y),s,c-i*inf,c+i*inf)。這里y是個(gè)數(shù)量符號(hào)。F=ilaplace(L,y,x)xt的函數(shù)，有ilaplace(L,y,x)等價(jià)于F(y)=int(L(y)*exp(x*y),y,c-i*inf,c+i*inf)，對(duì)y取積分。它將L 斯函數(shù)變換為時(shí)域函數(shù)F【例2-22】拉斯變換實(shí)例。求函數(shù)f(t)eat（a為實(shí)數(shù)、f2(t)tsint的拉斯變換1解：在命令窗口中輸入1symssymsts f1=exp(a*tf2=t- L1=laplace(f1); L1=laplace(f1); 程序運(yùn)行后，輸出的結(jié)果為L(zhǎng)1L1L2由運(yùn)行結(jié)果可知Lf1(t

s

(t)]

1 s2【例2-23】拉斯逆變換實(shí)例。求函數(shù)F1(s)

s(1s2

F2(s

s(s1)(s

的拉斯變換解：在命令窗口中輸入symst f1=ilaplace(F1); 斯逆變 f1=ilaplace(F1); 斯逆變 程序運(yùn)行后，輸出的結(jié)果為f2由運(yùn)行結(jié)果可知L1[F(S1cos(tL1[F(S2ete2t

t…0 Z提供了符號(hào)運(yùn)算工Z變換和Z逆變換Z變換的函數(shù)是ztrans，進(jìn)行Z逆變換的函數(shù)是iztrans。ztrans函數(shù)常用的調(diào)用格式如下F=ztrans(f)：函數(shù)返回獨(dú)立變量n關(guān)于符號(hào)向量f的Zztrans(fsymsum(f(n)/znn,0,inf)F=ztrans(fw)nfZ變換函數(shù)，只是用w代替了默認(rèn)的z：ztrasum(f(n)/wn,n,0,inf)。)=iztrans函數(shù)常用的調(diào)用格式如下f=iztrans(F)zF的Z逆變換函數(shù)，這是默認(rèn)的調(diào)用f=iztrans(Fk)：函數(shù)返回獨(dú)立變量kFZk代替