大學(xué)畢業(yè)設(shè)計MATLAB語言程序設(shè)計

1、MATLAB語言程序設(shè)計 系別：電子電氣工程系 班級：08級自動化班 姓名：XXX 學(xué)號：1. 編寫一個求圓的面積的函數(shù)文件。>> f(1)s = 3.1416>> f(2)s = 12.5664通過此題掌握了獨立文件與函數(shù)文件的區(qū)別。2. 三次拋物線的方程為：y=a*x3+b*x2+c*x試探討參數(shù)a，b和c對其圖形的影響。>> clear all;>> subplot(1,3,1)>> fplot('(-2:2)*x.3+x.2+x',-2 2) %繪制變量a=-2,-1,0,1,2時的圖形>> gri

2、d,axis('equal'),axis(-2 2 -4 4) %顯示網(wǎng)格，縱橫坐標軸保持一致圖>> %顯示區(qū)間>> subplot(1,3,2)>> fplot('x.3+(-2:2)x.2+x',-2 2) %繪制變量b=-2,-1,0,1,2時的圖形>> fplot('x.3+(-2:2)*x.2+x',-2 2) %繪制變量b=-2,-1,0,1,2時的圖形>> grid,axis('equal'),axis(-2 2 -4 4)>> subplot(

3、1,3,3)>> fplot('x.3+x.2+(-2:2)*x',-2 2) %繪制變量c=-2,-1,0,1,2時的圖形>> grid,axis('equal'),axis(-2 2 -4 4)>> gtext('a=2'),gtext('a=-2'),gtext('a=0') %用gtext命令在其他兩個子圖上標注字符程序運行結(jié)果如下：分析：由圖可得a=0時，此函數(shù)為拋物線，a>0時，隨著a的增大，其圖形越來越陡，a<0時，隨著a的減小，其圖形越來越陡，且a&l

4、t;0時，圖形的極值點增多。隨著b，c的增大，其圖形越來越陡。3設(shè)計GUI，通過調(diào)節(jié)滑塊可以畫出不同頻率的三角波形，同時學(xué)習(xí)法對Push Button、Checkbox、Slider、Axes、Popup Menu、Static Text控件的使用和操作。布局好的控件如下：添加的代碼如下：運行結(jié)果如下：通過此題，掌握了對Push Button、Checkbox、Slider、Axes、Popup Menu、Static Text控件的使用和操作，知道了圖形用戶界面設(shè)計工具GUIDE的強大功能。4 運算放大電路如圖所示，試分析放大器開環(huán)增益和頻率響應(yīng)對整個電路閉環(huán)頻率響應(yīng)的影響，并繪出曲線。程序

5、如下：>> z2=20,100,500*1000;z1=2000; %設(shè)置元件參數(shù)>> A0=2e6;w1=500;w2=2e6;w3=5e7;>> w=logspace(2,8); %設(shè)定頻率數(shù)組>> b=A0*w1*w2*w3;>> a=poly(-w1,-w2,-w3); %列出運算放大器分子分母系數(shù)向量>> A=polyval(b,j*w)./polyval(a,j*w); %求放大器開環(huán)頻率響應(yīng)>> for i=1:3z12(i)=z2(i)/z1;H(i,:)=-z12(i)./(1+(1+z12(

6、i)./A); %求放大器閉環(huán)響應(yīng)semilogx(w,abs(H(i,:),hold on %畫出頻率-增益曲線end>> v=axis;axis(v); %保持w坐標>> semilogx(w,abs(A)>> hold off運行結(jié)果如下：>> z2=20,100,500*1000;z1=2000; %設(shè)置元件參數(shù) >> A0=2e6;w1=50;w2=2e6;w3=5e7;>> w=logspace(2,8); %設(shè)定頻率數(shù)組>> b=A0*w1*w2*w3;>> a=poly(-w1,-w

7、2,-w3); %列出運算放大器分子分母系數(shù)向量>> A=polyval(b,j*w)./polyval(a,j*w); %求放大器開環(huán)頻率響應(yīng)>> for i=1:3z12(i)=z2(i)/z1;H(i,:)=-z12(i)./(1+(1+z12(i)./A); %求放大器閉環(huán)響應(yīng)semilogx(w,abs(H(i,:), hold on %畫出頻率-增益曲線end>> v=axis;axis(v); %保持w坐標>> semilogx(w,abs(A)>> hold off將w1減小10倍的H(w)的圖如下：由題可知，Z2越大

8、，越容易造成運算放大器的自激現(xiàn)象。消除自激可以通過減小w1，或增大w2，w3，在放大器已選定的情況下通常只能用加消振電容的方法減小w1。5分別利用命令和simulink模型求y=cos(t)dt的結(jié)果，其中初值分別為y1(0)=0和y2(0)=1。（1）simulink模型求解如下：y1（0）=0時：其結(jié)果如下：y2（0）=1時：其結(jié)果如下：（2）命令求解如下：>> clear>> f='cos(x)' %定義被積函數(shù)>> F=int(f) %求積分 F = sin(x)6許多工業(yè)控制過程都可以等效成二階環(huán)節(jié)，設(shè)計典型二階環(huán)節(jié)H(s)=20/

9、(1.6s2+4.4s+1)的模糊控制器，使系統(tǒng)盡快跟隨系統(tǒng)輸入。eudeNBNSZRPSPBNBPBPBPSPSZRNSPBPSPSZRZRZRPSPSZRZRNSPSPSZRZRNSNSPBZRZRNSNSNBu為輸出信號，e與de為輸入信號。u，e，de的取值范圍均在-6 6之間。在命令窗口下輸入命令fuzzy彈出如下界面：在Edit菜單下點擊Add Variable下的input，分別點擊input1和input2，在Name中改其名稱分別為e和de，將output1名字改為u得到如下圖形：雙擊輸入信號e，進入Membership Function Editor窗口界面，將其Range

10、與Display Range都改為-6 6,然后在membership function plots下刪掉原有的三個輸入波形，得圖如下：在Edit菜單下，打開Add Mfs窗口，修改Number of MFs為5，點擊OK按鈕。然后依次選中五個波形在Name中分別對其進行命名，對de與u做相同的操作，得出如下圖形：然后，返回到FIS Editor界面，雙擊Untitled進入Rule Editor界面，按照表格對其輸入輸出進行添加，完成后如下圖所示：然后點擊Close，在File菜單下選中Export中的To workplace，修改它的workspace variable 為mohu，如下所示：點擊OK按鈕，就可以在如下的仿真電路中進行仿真得出的仿真波形如下：若運行有錯誤提示，則在FIS Editor窗口下的Defuzzificat