【正文】 tr,ts]=mdlInitializeSizes(A,B,C,D)。C=[1,0]。B=[0。 A=[0,1。202。S186。181。205。207。208。193。210。L=laplace(F)第七章function[sys,x0,str,ts]=ww(t,x,u,flag)%182。39。)=maple(39。simplify(sin(x)^2+cos(x)^2)。Dy(0)=039。y(0)=139。D2y+2*Dy+2*y=039。2*xyz=139。xy+z=239。x+y+z=139。a*x^2+b*x+c39。syms sinv(s*eye(2)a)。) xF=ztrans(x*exp(x*10))=[0 1。),39。,39。dg=diff(g)=int(int(39。g=sin(y)。*39。:O39。Y=sin(X)。Y=sin(X)。b=[13,9,6,0]39。0,2,1,1。 A=[2,1,5,1。)。temps=interp2(x,h,T,xi,hi,39。xi=[0::10]。88,48,32,12,6。h=[0:30:60]39。*39。:O39。Y=sin(X)。Y=sin(X)。fxyz39。b=[13,9,6,0]39。0,2,1,1。=[2,1,5,1。z=u(3)。end i k第五章function f=fxyz(u)x=u(1)。k=0。i=i+1。k=0。end。while sumsum=sum+i。end s sum=0。end y s=0。n=100。m3=rem(m,10)。while sum(63)jcsum1(63) for m=100:999m1=fix(m/100)。end s=0。end y=0。n=100。m3=rem(m,10)。end m=100:999m1=fix(m/100)。while ifunction k=jcsum(n)k=0。function k=jcsum1(n)k=0。endfunction[s,p]=fcircle(r)s=pi*r*r。1,0,0])。colormap([0,1,0。c2(:,1:4)=3*ones(size(c2(:,1:4)))。hold on z2=z。c1=ones(size(z1))。z1=z。)。polar(theta,rho,39。theta=0::2*pi。)。),zlable(39。),ylable(39。xlabel(39。origin39。)。title(39。p39。z=4*sqrt(2)*sin(t)。x=8*cos(t)。plot(x,y),grid on。) x=0:pi/100:2*pi。,20)axis equal grid on xlabel(39。,39。plot(n,y,39。) n=(10:10)39。PA=poly(A)PPA=poly2str(PA,39。14 15 16。s39。s39。余多項式為39。A39。A(:)=4:5L=abs(A)3 islogical(L)X=A(L) A=[4,15,45,10,6。盡管如此，但是在實驗中依然容易把for循環(huán)跟C語言中的for語句弄混，最后經(jīng)過不懈努力下，終于弄明白了兩者之間的差別，使我能更好的運用這些指令語句。i=i+1。 elsey=1。 elseif n==139。 elseif n==039。if n39。請輸入一個正數(shù)n=39。!的函數(shù)程序，取n分別為80、10驗證其正確性(輸入n為負(fù)數(shù)時輸出出錯信息)。 [65] 39。39。39。[100][40]39。39。39。39。39。[85][76][95]39。39。39。end A A =39。E39。D39。C39。B39。A39。A+39。A(2,:)={85,76,95,100,40,65}。f39。e39。d39。c39。b39。a39。end A=cell(3,6)。c(k)={69+k}。clear,clc for k=1:10a(k)={89+k}。要求：（1）用switch語句實現(xiàn)。 end，要求輸出成績等級A+、A、B、C、D、E。if m0x1=(b+sqrt(m))/(2*a)x2=(bsqrt(m))/(2*a)elseif m==0x=(b)/(2*a)else39。為恒不等式39。m=b^24*a*c。請輸入一個數(shù)c=39。)。b=input(39。請輸入一個數(shù)a=39。并輸入幾組典型值加以檢驗。二、實驗設(shè)備及條件計算機一臺（）。實驗二 MATLAB程序設(shè)計一、實驗?zāi)康?。x2y2=1繪制平面曲線，并分析參數(shù)a對其形狀的影響。2在－54，繪制函數(shù)的曲面圖和等高線z=cosxcosyex2+y24其中x的21個值均勻分布在[5,5]范圍，y的31個值均勻分布在[0,10]，要求使用subplot（2，1，1）和subplot（2，1，2）將產(chǎn)生的曲面圖和登高圖畫在同一個窗口上。1In(x+1+x2),x0239。e y=237。x163。3，已知236。五、實驗內(nèi)容及步驟：1，設(shè)y=[+3sinx]cosx，在x＝0～2π區(qū)間取101點，繪制函數(shù)曲線。四、預(yù)習(xí)要求：1．復(fù)習(xí)二維與三維圖形的繪圖函數(shù)。為了保證圖象的顯示效果，一般還應(yīng)使用colormap函數(shù)設(shè)置圖象色圖。6，圖像處理1）imread和imwrite函數(shù) 這兩個函數(shù)分別用于將圖象文件讀入matlab工作空間，以及將圖象數(shù)據(jù)和色圖數(shù)據(jù)一起寫入一定格式的圖象文件。一般情況下，x，y，z是維數(shù)相同的矩陣。當(dāng)x，y是同維矩陣時，則以x，y對應(yīng)列元素為橫、縱坐標(biāo)分別繪制曲線，曲線條數(shù)等于矩陣的列數(shù)。二、實驗原理：1，二維數(shù)據(jù)曲線圖（1）繪制單根二維曲線 plot(x,y)。盡管學(xué)習(xí)上遇到了很多困難，結(jié)果也差強人意?？刂葡到y(tǒng)的 Simulink 仿真按圖建立系統(tǒng)的 Simulink 模型，對不同的輸入信號進行仿真，改變參數(shù)，觀察不同的仿真結(jié)果。第一篇：基于 Matlab 的離散控制系統(tǒng)仿真2014 / 2015 學(xué)年第 1 學(xué)期計算機控制技術(shù) 實班 級 學(xué) 生 指 導(dǎo) 驗 報 告學(xué) 號 1108030301 姓 名 蔡 夢 教 師 張 坤 鰲實驗二 基于 Matlab 的離散控制系統(tǒng)仿真一、實驗?zāi)康暮鸵螅簩W(xué)習(xí)使用 Matlab 的命令對控制系統(tǒng)進行仿真的方法學(xué)習(xí)使用 Matlab 中的 Simulink 工具箱進行系統(tǒng)仿真的方法二、實驗環(huán)境X86系列兼容型計算機，Matlab軟件三、實驗原理控制系統(tǒng)命令行仿真1）建立如圖所示一階系統(tǒng)控制模型并進行系統(tǒng)仿真：一階系統(tǒng)閉環(huán)傳遞函數(shù)為G（S）=s1+333s=s+3，轉(zhuǎn)換為離散系統(tǒng)脈沖傳遞函數(shù)并仿真。2）建立如圖所示二階系統(tǒng)控制模型并進行系統(tǒng)仿真：52s(s+2**5)25251+s(s+2**5)=s2+2**5s+52，二階系統(tǒng)閉環(huán)傳遞函數(shù)為G（S）=轉(zhuǎn)換為離散系統(tǒng)脈沖傳遞函數(shù)并仿真，改變參數(shù)，觀察不同的系統(tǒng)的仿真結(jié)果。將上述系統(tǒng)離散化并仿真，觀察仿真結(jié)果四、實驗步驟根據(jù)實驗原理對控制系統(tǒng)進行軟件仿真觀察記錄輸出的結(jié)果，與理論計算值相比較自行選擇參數(shù)，練習(xí)仿真方法，觀察不同的仿真結(jié)果5252s(s+2**5)s(s+2**5)=52521+1+s(s+2**5)s(s+2**5)進行軟二階系統(tǒng)閉環(huán)傳遞函數(shù)為G（S）=件仿真如下圖：分別進行離散仿真：五、實驗心得針對這次實驗設(shè)計，我通過各種渠道，上課認(rèn)真學(xué)習(xí)，請教老師、上網(wǎng)搜索，圖書館查閱，詢問同學(xué)等學(xué)習(xí)到了很多知識，一步步了解最少拍控制系統(tǒng)設(shè)計，鍛煉了自我學(xué)習(xí)能力。但我們在不斷處理困難的過程中磨練了處理事物的能力和耐心，也讓同學(xué)間學(xué)會了互相學(xué)習(xí)，共享資源第二篇：MATLAB與控制系統(tǒng)仿真實驗報告《MATLAB與控制系統(tǒng)仿真》實驗報告20132014學(xué)年 第 1 學(xué)期專業(yè)： 班級： 學(xué)號： 姓名：實驗三 MATLAB圖形系統(tǒng)一、實驗?zāi)康模?。?）繪制多根二維曲線 plot(x,y)當(dāng)x是向量，y是有一維與x同維的矩陣時，則繪制多根不同顏色的曲線。（3）含有多個輸入?yún)?shù)的plot函數(shù) plot(x1,y1,x2,y2,…,xn,yn)（4）具有兩個縱坐標(biāo)標(biāo)度的圖形 plotyy(x1,y1,x2,y2)2,圖形標(biāo)注與坐標(biāo)控制 1）title(圖形名稱)； 2）xlabel（x軸說明）3）ylabel（y軸說明）4）text（x，y圖形說明）5）legend（圖例1，圖例2，…）6）axis（[xmin xmax ymin ymax zmin zmax]）3, 圖形窗口的分割 subplot（m,n,p）4,三維曲線plot3（x1,y1,z1,選項1，x2,y2，選項2，…,xn,yn,zn,選項n）5，三維曲面mesh(x,y,z,c)與surf(x,y,z,c)。X，y是網(wǎng)格坐標(biāo)矩陣，z是網(wǎng)格點上的高度矩陣，c用于指定在不同高度下的顏色范圍。2）image和imagesc函數(shù) 這兩個函數(shù)用于圖象顯示。三、實驗儀器和設(shè)備：計算機一臺（）。2．復(fù)習(xí)圖形輔助操作。21+x2，已知y1=x2,y2=cos(2x),y3=y1*y2,完成下列操作：（1）在同一坐標(biāo)系下用不同的顏色和線型繪制三條曲線；（2）分別用條形圖、階梯圖、桿圖和填充圖繪制三條曲線。x+p,239。02239。239。238。z=x2+y2+sin(xy)的曲面及等高線圖。+a25a2四、心得體會：通過這次實驗我能熟練掌握二維和三維圖以及其他特殊圖形的制作，弄清楚了基本的圖形操作規(guī)則，大大加深了我對matlab的興趣。三、實驗內(nèi)容+bx+c=0的根的函數(shù)（這個方程不一定為一元二次方程，因a、b、c的不同取值而定），這里應(yīng)根據(jù)a、b、c的不同取值分別處理，有輸入?yún)?shù)提示，當(dāng)a=0,b=0,c~=0時應(yīng)提示“為恒不等式!”。clear,clc a=input(39。)。請輸入一個數(shù)b=39。c=input(39。)。if a==0if b==039。end endm=b^24*a*c。不存在正實根39。其中100分為A+，90分～99分為A，80分～89分為B，70分～79分為C，60分～69分為D，60分以下為E。（2）輸入百分制成績后要判斷該成績的合理性，對不合理的成績應(yīng)輸出出錯信息。b(k)={79+k}。d(k)={59+k}。A(1,:)={39。,39。,39。,39。,39。,39。}。for k=1:6switch A{2,k}case 100r=39。case ar=39。case br=39。case cr=39。case dr=39。otherwiser=39。endA(3,k)={r}。a39。b39。c39。B39。C39。A39。d39。e39。A+39。E39。f39。D39。clear,clc n=input(39。)。輸入錯誤39。n!=039。n!=139。for i=1:1:ny=y*i。endy end 請輸入一個正數(shù)n=89ans =輸入錯誤 請輸入一個正數(shù)n=0ans =n!=0 請輸入一個正數(shù)n=1ans =n!=1 請輸入一個正數(shù)n=3y =6 請輸入一個正數(shù)n=10y =3628800四、實驗心得體會：通過本次實驗課，我能熟練運用for循環(huán)語句，switch條件語句以及if條件語句的新用法，和在C中的區(qū)別。第三篇：控制系統(tǒng)的Matlab仿真與設(shè)計課后答案MATLAB課后習(xí)題答案 x=[15 22 33 94 85 77 60] x(6)x([1 3 5])x(4:end)x(find(x70)) A=zeros(2,5)。56,0,17,45,0] find(A=10amp。cr=39。disp([cq,poly2str(q,39。)]),disp([cr,poly2str(r,39。)]) A=[11 12 13。17 18 19]。s39。y=abs(n)。r.39。MarkerSize39。n39。y=2*exp(*x).*sin(2*pi*x)。 t=0:pi/50:2*pi。y=4*sqrt(2)*sin(t)。plot3(x,y,z,39。)。Line in 3D Space39。text(0,0,0,39。)。X39。Y39。Z39。grid。rho=sin(2*theta).*cos(2*theta)。k39。[x,y,z]=sphere(20)。z1(:,1:4)=NaN。surf(3*x,3*y,3*z1,c1)。c2=2*ones(size(z2))。surf(*x,*y,*z2,c2)。,0,0。grid on hold off 第四章function f=factor(n)if nf=factor(n1)*n。p=2*pi*r。i=0。for i=0:n k=k+2^i。m2=rem(fix(m/10),10)。ifm==m1*m1*m1+m2*m2*m2+m3*m3*m3 disp(m)end end[s,p]=fcircle(10)=0。for i=1:ny=y+1/i/i。for i=1:5s=s+factor(i)。i=1。m2=rem(fix(m/10),10)。if m==m1*m1*m1+m2*m2*m2+m3*m3*m3disp(m)end end y=0。for i=1:ny=y+1/i/i。for i=1:5s=s+factor(i)。i=1。i=i+1。n=i2 i=0。while ik=k+2^i。end k ii=0。for i=0:63k=k+2^i。y=u(2)。f=x+y.^2./x/4+z.^2./y+2./z。1,5,0,7。1,6,1,4]。x=Ab[U,fmin]=fminsearch(39。,[,])=linspace(0,2*pi,50)。P=polyfit(X,Y,3)AX=linspace(0