matlab實(shí)驗(yàn)內(nèi)容答案

1、實(shí)驗(yàn)報(bào)告說(shuō)明：matlab 課程實(shí)驗(yàn)需撰寫(xiě) 8 個(gè)實(shí)驗(yàn)報(bào)告， 每個(gè)實(shí)驗(yàn)報(bào)告內(nèi)容寫(xiě)每次 實(shí)驗(yàn)內(nèi)容中標(biāo)號(hào)呈黑體大號(hào)字顯示的題目。第一次實(shí)驗(yàn)內(nèi)容：實(shí)驗(yàn)一 MATLAB 運(yùn)算基礎(chǔ)一、實(shí)驗(yàn)?zāi)康?. 熟悉啟動(dòng)和退出 MATLAB 的方法。2. 熟悉 MATLAB 命令窗口的組成。3. 掌握建立矩陣的方法。4. 掌握 MATLAB 各種表達(dá)式的書(shū)寫(xiě)規(guī)則以及常用函數(shù)的使用。二、實(shí)驗(yàn)內(nèi)容1.先求下列表達(dá)式的值，然后顯示MATLAB工作空間的使用情況并保存全部變量。1) z1 2s1ine8252)z21 ln( x 1 x2 ) ，其中212i20.4550.3a0.3 a3)z3eegsin( a0.3),

2、a3.0,2.9,2.8,L ,2.8, 2.9,3.02提示：利用冒號(hào)表達(dá)式生成 a 向量，求各點(diǎn)的函數(shù)值時(shí)用點(diǎn)乘運(yùn)算。t0t14) z4 t 211t2 ，其中 t 0t22t 12t3提示：用邏輯表達(dá)式求分段函數(shù)值。2.已知12 34 41 3 1A 34 7 87 ， B2033 65 73 2 7求下列表達(dá)式的值：（1）A+6=B和 A-B+I（其中 I 為單位矩陣）。（2）A*B 和 A.*B 。（3）A3 和 A.3 。（4）A/B 和 BA 。（5）A， B和 A(1,3,;);B2 。3.設(shè)有矩陣 A 和 B1234530166789101769A 1112131415 ，

3、B023416171819209702122232425413111）求它們的乘積 C 。2）將矩陣 C的右下角 3×2 子矩陣賦給 D3）查看 MATLAB工作空間使用情況。4.完成下列操作：（1）求100，999之間能被 21 整除的數(shù)的個(gè)數(shù)。 提示：先利用冒號(hào)表達(dá)式，再利用 find 和 length 函數(shù)。（2）建立一個(gè)字符串向量，刪除其中的大寫(xiě)字母。提示：利用 find 函數(shù)和空矩陣。第二次實(shí)驗(yàn)內(nèi)容：實(shí)驗(yàn)三 選擇結(jié)構(gòu)程序設(shè)計(jì)一、實(shí)驗(yàn)?zāi)康?. 掌握建立和執(zhí)行 M 文件的方法。2. 掌握利用 if 語(yǔ)句實(shí)現(xiàn)選擇結(jié)構(gòu)的方法。3. 掌握利用 switch 語(yǔ)句實(shí)現(xiàn)多分支選擇結(jié)構(gòu)的方

4、法。4. 掌握 try 語(yǔ)句的使用。二 、實(shí)驗(yàn)內(nèi)容1. 求下列分段函數(shù)的值。2xx 6,x 0且 x 32y x25x 6,0 x 10,x 2且x 32 xx1其他要求：(1) 用 if 語(yǔ)句實(shí)現(xiàn)，分別輸出 x= , , , , , , 時(shí)的 y值。 提示： x 的值從鍵盤(pán)輸入，可以是向量。( 2) 仿照實(shí)驗(yàn)一第 1 題第 4 小題，用邏輯表達(dá)式實(shí)現(xiàn)，從而體會(huì)MATLAB邏輯表達(dá)式的一種應(yīng)用。clearx=input('input the value of x :');if x<0 & x=3y=x*x+x-6;elseif x>=0 & x<

5、;10 & x=2 &x=3 y=x*x-5*x+6;elsey=x*x-x-1;endy2輸入一個(gè)百分制成績(jī)，要求輸出成績(jī)等級(jí)A,B,C,D,E 。其中 90100分為A，8089分為 B，7079 分為 C， 60 69 分為 D，60 分以下為 E。要求：( 1) 用 switch 語(yǔ)句實(shí)現(xiàn)。( 2) 輸入百分制成績(jī)后要判斷該成績(jī)的合理性，對(duì)不合理的成績(jī)應(yīng)輸出錯(cuò)誤信息。if 語(yǔ)句實(shí)現(xiàn)：clearx=input('input score:');if x<0disp('Sorry,The fail score!')elseif x<

6、60disp('E')elseif x<70disp('D')elseif x<80disp('C')elseif x<90disp('B')elseif x<=100disp('A')elsedisp('Sorry,The fail score!')end switch 語(yǔ)句實(shí)現(xiàn)：clear score=input('input the score:');switch score case num2cell(0:59) disp('E')ca

7、se num2cell(60:69)disp('D')case num2cell(70:79)disp('C')case num2cell(80:89)disp('B')case num2cell(90:100) disp('A') otherwise disp('Sorry!It''s the fail score,input the score again') end 3. 建立 5×6 矩陣，要求輸出矩陣第 n 行元素。當(dāng) n 值超過(guò)矩陣的行數(shù)時(shí)，自動(dòng)轉(zhuǎn)為輸出矩陣最 后一行元素，并給

8、出出錯(cuò)信息。clear A=1,1,1,1,1,1;2,2,2,2,2,2;3,3,3,3,3,3;4,4,4,4,4,4;5,5,5,5,5,5;n=input(' 請(qǐng)輸入輸出元素行數(shù) n=');tryB=A(n,:);catchB=A(5,:);lasterrendB實(shí)驗(yàn)四 循環(huán)結(jié)構(gòu)程序設(shè)計(jì)一、實(shí)驗(yàn)?zāi)康?. 掌握利用 for 語(yǔ)句實(shí)現(xiàn)循環(huán)結(jié)構(gòu)的方法。2. 掌握利用 while 語(yǔ)句實(shí)現(xiàn)循環(huán)結(jié)構(gòu)的方法。3. 熟悉利用向量運(yùn)算來(lái)代替循環(huán)操作的方法。、實(shí)驗(yàn)內(nèi)容1.根據(jù)6122232 L12，求 的近似值。當(dāng) n分別取 100、1 000、10 000 時(shí)，結(jié)果是 n多少要求：分別

9、用循環(huán)結(jié)構(gòu)和向量運(yùn)算（使用 sum 函數(shù)）來(lái)實(shí)現(xiàn)。應(yīng)用 sum 的方法： clear n=100;i=1:n; y=1./i.2; y1=sum(y); pai=sqrt(6*y1)第二種方法： clear y=0;n=input(' 輸入數(shù)據(jù) n： ') for i=1:n y=y+1/(i*i);end pai=sqrt(6*y)1 1 12.根據(jù) . y 1 L ，求：3 5 2n 1(1) y<3時(shí)的最大 n 值。(2) 于( 1)的 n 值對(duì)應(yīng)的 y值。clear y=0;n=1000;for i=1:1000 y=y+1/(2*i-1);if y>=3b

10、reak;endi=i+1;endi3. 已知f1 1f2 0,n 3f3 1fn fn 1 2fn 2 fn 3求 f1 f100 中：（ 1） 最大值、最小值、各數(shù)之和。（ 2） 正數(shù)、零、負(fù)數(shù)的個(gè)數(shù)。f=1,0,1;for n=4:100f(n)=f(n-1)-2*f(n-2)+f(n-3); end k=length(find(f<0) k=length(find(f>0) k=length(find(f=0) 4若兩個(gè)連續(xù)自然數(shù)的乘積減 1 是素?cái)?shù)， 則稱這兩個(gè)連續(xù)自然數(shù)是親密數(shù)對(duì)， 該素?cái)?shù)是親密素?cái)?shù)。 求2，50區(qū)間內(nèi)：(1) 親密數(shù)對(duì)的對(duì)數(shù)；(2) 與上述親密數(shù)對(duì)對(duì)應(yīng)

11、的所有親密素?cái)?shù)之和。實(shí)驗(yàn)五 函數(shù)與文件一、實(shí)驗(yàn)?zāi)康?. 掌握定義和調(diào)用 ,MATLAB 函數(shù)的方法。2. 掌握 MATLAB文件的基本操作。二、實(shí)驗(yàn)內(nèi)容1. 定義一個(gè)函數(shù)文件，求給定復(fù)數(shù)的指數(shù)、對(duì)數(shù)、正弦和余弦，并在命令文件中調(diào)用該函數(shù)文件。function a,b,c,d=susb(x) a=exp(x);b=sin(x);c=log10(x);d=cos(x);調(diào)用：x=1+2ia,b,c,d=susb(x)2. 一個(gè)自然數(shù)是素?cái)?shù)，且它的各位數(shù)字位置經(jīng)過(guò)任意對(duì)換后仍為素?cái)?shù)，則稱是絕對(duì)素?cái)?shù)。 例如 13 是絕對(duì)素?cái)?shù)。試求所有兩位的絕對(duì)素?cái)?shù)。要求：定義一個(gè)判斷素?cái)?shù)的函數(shù)文件。functionq

12、=apnum(p); % 以下為函數(shù)文件，需單獨(dú)保存for i=2:9n=find(rem(p,i)=0);p(n)=;% 找出能被 i 整除的 p 的位置% 消去能被 i 整除的 pend q=10*rem(p,10)+fix(p/10);for j=2:9% p 的個(gè)位和十位調(diào)換，調(diào)換后為 qm=find(rem(q,j)=0);q(m)=;% 找出能被 j 整除的 q 的位置% 消去能被 j 整除的 qend調(diào)用： clear all p=10:99;q=apnum(p)3. 已知 yf(40)f (30) f (20)1)當(dāng) f (n) n 10ln(n2 5)時(shí)，求 y 的值。2)當(dāng)

13、 f (n) 1 2 2 3 3 4 n (n 1) 時(shí)，求 y的值。第三次實(shí)驗(yàn)內(nèi)容：實(shí)驗(yàn)六 高層繪圖操作一、實(shí)驗(yàn)?zāi)康?. 掌握繪制二維圖形的常用函數(shù)。2. 掌握繪制三維圖形的常用函數(shù)。二、實(shí)驗(yàn)內(nèi)容21. 已知 y1 x ,y2 cos(2x), y3 y1 y2 完成下列操作：( 1)在同一坐標(biāo)系下用不同的顏色和線型繪制三條曲線。(2) 以子圖形式繪制 3 條曲線。(3) 分別用條圖形、階梯圖、桿圖和填充圖繪制3 條曲線。12. 分別用 plot 和 fplot 函數(shù)繪制函數(shù) y sin 的曲線，分析兩曲線的差別。3. 繪制函數(shù)的曲面圖和等高線。 已知 z (x2 2x)e x y xy 提

14、示：繪制三維曲面圖，首先要選定一平面區(qū)域并在該區(qū)域產(chǎn)生網(wǎng)絡(luò)坐標(biāo)矩陣。在做本題前，先分 析并上機(jī)驗(yàn)證下列命令的執(zhí)行結(jié)果，從中體會(huì)產(chǎn)生網(wǎng)絡(luò)坐標(biāo)矩陣的方法。x,y=meshgrid(-1:2,1:5)3. 已知 6<x<30， 15<y<36，求不定方程 2x+5y=126 的整數(shù)解。x=7:29; y=16:35;x,y=meshgrid(x,y); %在7,29 ×16,35區(qū)域生成網(wǎng)格坐標(biāo)z=2*x+5*y; k=find(z=126); x(k),y(k)%找出解的位置%輸出對(duì)應(yīng)位置的x,y 即方程的解5. 裁掉表現(xiàn)函數(shù) z=sin(y)cos(x)三維曲面

15、圖中 z>部分。 x=0:2*pi;x,y=meshgrid(x);z=sin(y).*cos(x);I,J=find(z>for ii=1:length(I)z(I(ii),J(ii)=NaN;endsurf(x,y,z);實(shí)驗(yàn)七低層繪圖操作一、實(shí)驗(yàn)?zāi)康?. 掌握?qǐng)D形對(duì)象屬性的基本操作。2. 掌握利用圖形對(duì)象進(jìn)行繪圖操作的方法。二、實(shí)驗(yàn)內(nèi)容1. 建立一個(gè)圖形窗口，使之背景顏色為紅色，并在窗口上保留原有的菜單項(xiàng)，而且在按下鼠標(biāo)器 的左鍵之后顯示出 Left Button Pressed 字樣。hf=figure('Color',1,0,0,'WindowBu

16、ttonDownFcn','disp(''Left Button Pressed.'')');2 2x2. 先利用缺省屬性繪制曲線 y x2e2x ，然后通過(guò)圖形句柄操作來(lái)改變曲線的顏色、線型和線寬， 并利用文字對(duì)象給曲線添加文字標(biāo)注 y x2e2x 。實(shí)驗(yàn)十二 菜單設(shè)計(jì)一、實(shí)驗(yàn)?zāi)康?. 了解圖形用戶界面的特點(diǎn)。2. 掌握菜單設(shè)計(jì)的方法。3. 掌握對(duì)話框設(shè)計(jì)的方法。二、實(shí)驗(yàn)內(nèi)容1. 設(shè)計(jì)菜單。菜單條僅有 File菜單項(xiàng)， File下有 New 、Plot和Exit等 3個(gè)選項(xiàng)。選擇 New時(shí)利用 Edit命令建立一 個(gè)新的 M 文件。選擇

17、 Plot 將顯示下一級(jí)菜單，其中有 Sine Wave和 Cosine Wave兩個(gè)子菜單項(xiàng)，且 若選擇了其中的 Sine Wave子菜單項(xiàng)，則將打開(kāi)一個(gè)新的圖形窗口并顯示出正弦曲線。若選擇了其 中的 Cosine Wave 子菜單項(xiàng)，則將打開(kāi)一個(gè)新的圖形窗口并顯示出余弦曲線。 如果選擇 Exit 菜單項(xiàng)， 則將關(guān)閉窗口并推出用戶系統(tǒng)回到MATLAB命令窗口。2. 為圖形窗口建立快捷菜單，用以控制窗口的背景顏色和大小。4. 設(shè)計(jì)一個(gè)對(duì)話框，其中有一個(gè)編輯框和按鈕，當(dāng)單擊按鈕時(shí)，使編輯框的內(nèi)容加5。第四次實(shí)驗(yàn)內(nèi)容：實(shí)驗(yàn)二 矩陣分析一、實(shí)驗(yàn)?zāi)康?. 掌握生成特殊矩陣的方法。2. 掌握矩陣分析的方

18、法。二、實(shí)驗(yàn)內(nèi)容1. 產(chǎn)生 5階希爾伯特矩陣 H和5階帕斯卡矩陣 P，且求其行列式的值 Hh和 Hp以及它們的條件數(shù) Th和 Tp，判斷哪個(gè)矩陣性能更好，為什么H=hilb(5)%產(chǎn)生 55 階希爾波特矩陣P=pascal(5)%產(chǎn)生 55 階帕斯卡矩陣Hh=det(H);%計(jì)算行列式的值%計(jì)算條件數(shù)Ph=det(P);Th=cond(H);Tp=cond(P);HhPhThTp2. 建立一個(gè) 5×5 矩陣，求它的行列式值、跡、秩和范數(shù)。A=1 2 3 4 5;2 3 4 5 6;3 4 5 6 7;4 5 6 7 8;5 6 7 8 9;det(A)trace(A)rank(A)

19、norm(A)3. 已知29 6 18A 20 5 128 8 5求 A 的特征值及特征向量，并分析其數(shù)學(xué)意義。實(shí)驗(yàn)八 數(shù)據(jù)處理和多項(xiàng)式計(jì)算、實(shí)驗(yàn)?zāi)康?.3.掌握數(shù)據(jù)統(tǒng)計(jì)和分析的方法。 掌握數(shù)值插值與曲線擬合的方法及其應(yīng)用。 掌握多項(xiàng)式的常用運(yùn)算。、實(shí)驗(yàn)內(nèi)容1. 將 100 個(gè)學(xué)生 5 門(mén)功課的成績(jī)存入矩陣 P 中，進(jìn)行如下處理：(1)分別求每門(mén)課的最高分、最低分及相應(yīng)學(xué)生序號(hào)。 (2)分別求每門(mén)課的平均分合標(biāo)準(zhǔn)方差。(3) 5 門(mén)課總分的最高分、最低分及相應(yīng)學(xué)生序號(hào)。(4) 將 5 門(mén)課總分按從大到小順序存入 zcj 中，相應(yīng)學(xué)生序號(hào)存入 xsxh。 提示：上機(jī)調(diào)試時(shí)，為避免輸入學(xué)生成績(jī)的麻

20、煩，可用取值范圍在 45,95 之間的隨機(jī)矩陣來(lái)表示學(xué) 生成績(jī)。A=rand(5,100);P=100*A;% 求每門(mén)功課的最高分、序號(hào)% 求每門(mén)功課的最低分、序號(hào)% 求每門(mén)功課的平均值% 求每門(mén)功課的標(biāo)準(zhǔn)方差% 求五門(mén)功課的總成績(jī)% 求總成績(jī)的最高分% 求總成績(jī)的最低分Y,U1=max(P,2);X,U2=min(P,2);YU1XU2Z1=mean(P,2)Z2=std(P,1,2)B1=cumsum(P);B=B1(5,:);YZ,UZ1=max(B,2) XZ,UZ2=min(B,2) zcj=-sort(-B,2)2. 已知 lg(x)在1,101區(qū)間 11個(gè)整數(shù)采樣點(diǎn)的函數(shù)值如表

21、2所示。表 2 lg(x)在 10 個(gè)采樣點(diǎn)的函數(shù)值x1112131415161718191101Lg(x)0試求 lg(x)的 5次擬合多項(xiàng)式 p(x),并分別繪制出 lg(x)和 p(x)在1,101區(qū)間的函數(shù)曲線。 x=1:10:101;y=0,;P,S=polyfit(x,y,5)plot(x,y,'k*',x,polyval(P ,x),'k-')4 3 223. 有 3個(gè)多項(xiàng)式 p1(x)x42x34x25,p2(x)x 2, p3(x)x22x 3,試進(jìn)行下列操作：求 p(x) p1(x) p2(x) p3(x) 。1) 求 p(x) 的根。2)

22、當(dāng) x 取矩陣 A 的每一元素時(shí)，求 P( x)的值。其中1 1.2 1.4 A 0.75 23.5052.5(3) 當(dāng)以矩陣 A 為自變量時(shí)，求 p1=1,2,4,0,5 p22=1,2;p2=0,0,0,p22; p33=1,2,3;p3=0,0,p33; p=p1+p2.*p3 x=roots(p) A=-1,;,2,;0,5,;G=poly(A)第五次實(shí)驗(yàn)內(nèi)容：P（x）的值。其中 A 的值與（ 3）相同。% 求多項(xiàng)式的根% X=A 時(shí) P（X）的值實(shí)驗(yàn)九 符號(hào)計(jì)算基礎(chǔ)與符號(hào)微積分一、實(shí)驗(yàn)?zāi)康?. 掌握定義符號(hào)對(duì)象的方法。2. 掌握符號(hào)表達(dá)式的運(yùn)算法則以及符號(hào)矩陣運(yùn)算。3. 掌握求符號(hào)函

23、數(shù)極限及導(dǎo)數(shù)的方法。4. 掌握求符號(hào)函數(shù)定積分和不定積分的方法。、實(shí)驗(yàn)內(nèi)容1. 已知 x=6,y=5,利用符號(hào)表達(dá)式求 zx13 x y提示：定義符號(hào)常數(shù) x=sym( 6' ),y=sym(。 5') x=sym('6');y=sym('5');z=(x+1)/(sqrt(3+x)-sqrt(y)2. 用符號(hào)方法求下列極限或?qū)?shù)。1) lxim0sinx x(etanx1) 2(e 1)3sin xsyms x;f=(x*(exp(sin(x)+1)-2*(exp(tan(x)-1)/(sin(x)3; limit(f,x,0)2）已知 Atc

24、osxt3lnxdA，分別求 dA dxd2A,d2A2dt 2 dxdt%求 X 一階導(dǎo)數(shù)%求 X 二階導(dǎo)數(shù)syms a t x;A=ax,t3;t*cos(x),log(x); diff(A,x) diff(A,x,2)3. 用符號(hào)方法求下列積分。1)dx1 x4 x8x=sym('x'); f=1/(1+x4+x8); int(f)x=sym('x'); f=(x2+1)/(x4+1); int(f,0,inf)實(shí)驗(yàn)十 級(jí)數(shù)與方程符號(hào)求解一、實(shí)驗(yàn)?zāi)康?. 掌握級(jí)數(shù)求和的方法。2. 掌握將函數(shù)展開(kāi)為泰勒級(jí)數(shù)的方法。3. 掌握微分方程符號(hào)求解的方法。4. 掌握

25、代數(shù)方程符號(hào)求解的方法。、實(shí)驗(yàn)內(nèi)容1. 級(jí)數(shù)符號(hào)求和。10計(jì)算 Sn 1 2n 1n=sym('n');S=symsum(1/(2*n-1),n,1,10)21) 求級(jí)數(shù) n2xn 1 之和函數(shù)，并求nn 之和。nn 1n 1%求 N的從 1 到無(wú)窮的和syms n x;S2=symsum(n2)*x(n-1),n,1,inf)函數(shù)S3=symsum(n2)/5n)2. 將 ln(x) 在 x=1 處按 5 次多項(xiàng)式展開(kāi)為泰勒級(jí)數(shù)。x=sym('x'); f=log(x); taylor(f,x,6,1) 2 求微分方程的符號(hào)解。d2ydx2k2 y 0y(0)

26、 a y'(0) ba,b,k為任意常數(shù)syms x y a b k; m=dsolve('D2y+k2*y','y(0)=a','Dy(0)=b','x')4. 求下列方程和方程組的符號(hào)解。1) 3xex 5sin x 78.5 0x=solve('3*x*exp(x)+5*sin(x)','x')2)x2 y2 100 03x 5y 8 0x y=solve('sqrt(x2+y2)-100','3*x+5*y-8','x,y')53) l

27、n(1 x) 21 sinxx=solve('log(1+x)-5/(1+sin(x)=2','x')4) x2 9 x 1 1 0x=solve('x2+9*sqrt(x+1)-1','x')第六次實(shí)驗(yàn)內(nèi)容： 實(shí)驗(yàn)十三 Simulink 的應(yīng)用 一、實(shí)驗(yàn)?zāi)康?. 熟悉 Simulink 的操作環(huán)境并掌握構(gòu)建系統(tǒng)模型的方法。2. 掌握 Simulink 中子系統(tǒng)模塊的建立與封裝技術(shù)。3. 對(duì)簡(jiǎn)單系統(tǒng)所給出的數(shù)學(xué)模型能轉(zhuǎn)化為系統(tǒng)仿真模型并進(jìn)行仿真分析。二、實(shí)驗(yàn)內(nèi)容1. 假設(shè)從實(shí)際應(yīng)用領(lǐng)域(力學(xué)、電學(xué)、生態(tài)或社會(huì))中，抽象出有初始狀態(tài)

28、為 0 的二階微分方程x 0.2x 0.4 0.2u(t) ，u(t)是單位階躍函數(shù)。 用積分器直接構(gòu)造求解微分方程的模型并仿真。2. 利用傳遞函數(shù)模塊建模。仍以下面的二階微分方程為例 x 0.2x 0.4 0.2u(t)x 0.2x 0.4 0.2u(t)3. 利用狀態(tài)方程模塊建模。仍以下面的二階微分方程為例第七次實(shí)驗(yàn)內(nèi)容：實(shí)驗(yàn)十四 綜合實(shí)驗(yàn)（ 1）MATLAB解決控制系統(tǒng)仿真分析與設(shè)計(jì)的方法。一、實(shí)驗(yàn)?zāi)康?綜合運(yùn)用所學(xué)知識(shí)，掌握利用二、實(shí)驗(yàn)內(nèi)容1用 Matlab 表示傳遞函數(shù)為2s 3324s 3s 2s 1 的系統(tǒng)。2用 Matlab 表示傳遞函數(shù)為2(s 3)s(s 1)(s 2)的系統(tǒng)。3. 用 Matlab 表示狀態(tài)空間表達(dá)式為 x3010x10001x20u123x31 的系統(tǒng)。x1y 1 0 0 x2x3A=0 1 0;0 0 1;-1 -2 -3;B=0;0;1;C=1 0 0;D=0;Gf=ss(A,B,C,D)4. 已知某控制系統(tǒng)的傳遞函數(shù)為G(s)12 s2 3s 2，求Matlab 描述的傳遞函數(shù)模型及零極點(diǎn)增益模型。num=1; den=1 3 2; Gtf=tf(num,den) Gzpk=zpk(Gtf