東北林業(yè)大學(xué)機(jī)電工程學(xué)院電氣matlab實(shí)驗(yàn)內(nèi)容答案

實(shí)驗(yàn)一MATLAB運(yùn)算基礎(chǔ)先求下列表達(dá)式的值，然后顯示MATLAB工作空間的使用情況并保存全部變量。1)12sin85。

zl=1+e2clear；z1=2*sin(85/180*pi)/(1+exp(2))2)ln(x+v/l+x2)其中=—0.451+2i5clear；x=[2,1+2*i;-0.45,5];z2=log(x+sqrt(1+xA2))/2e0.3a—e—0.3a3)z3= ESin(a+0.3),a=—3.0,—2.9,—2.&…,2.8,2.9,3.03)2clear；a=(-3:0.1:3);z3=((exp(0.3*a)-exp(-0.3*a))/2).*sin(a+0.3)12 0<t<1(4)z4=<t2—1 1<t<2，其中t=0：0.5：2.512—2t—1 2<t<3clearfort=0:0.5:2.5ift>=0&t<1z4=tA2elseift>=1&t<2z4=tA2-1elsez4=tA2-2*t-1endend2.已知「1234—4"「13—1_A=34787，B=20336573—27求下列表達(dá)式的值：(1)A+6=B和A-B+I(其中I為單位矩陣)。(2)A*B和A.*B。(3)AA3和AA.3。(4)A/B和B\A。(5)[A,B]和[A([1,3],;);BA2]。clearA=[12,34,-4;34,7,87;3,65,7];B=[1,3,-1;2,0,3;3,-2,7];Z1=A+6*BI=eye(3,3)Z11=A-B+IAA3%ThethirdA.A3A/B%TheforthB\A[A,B]%Thefifth[&[1,3],:)田人2]3?設(shè)有矩陣A和B「1 2 3 4 5-「3 0 166 7 8 9 1017—6 9A=1112131415B=0 23—416171819209 7 021222324254 1311(1)求它們的乘積C。將矩陣C的右下角3X2子矩陣賦給D查看MATLAB工作空間使用情況。clearA=[1,2,3,4,5;6,7,8,9,10;11,12,13,14,15;16,17,18,19,20;21,22,23,24,25];B=[3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11];C=A*BD=C(3:5,2:3)whowhos4.完成下列操作：求[100，999]之間能被21整除的數(shù)的個(gè)數(shù)。cleara=100:999;k=find(rem(a,21)==0);length(k)建立一個(gè)字符串向量，刪除其中的大寫(xiě)字母。clearch='FlfkljFGDKJ';k=find(ch>='A'&ch<='Z');ch(k)=0實(shí)驗(yàn)二矩陣分析產(chǎn)生5階希爾伯特矩陣H和5階帕斯卡矩陣P,且求其行列式的值Hh和Hp以及它們的條件數(shù)Th和Tp,判斷哪個(gè)矩陣性能更好，為什么？H=hilb(5)P=pascal(5)Hh=det(H);Ph=det(P);Th=cond(H);Tp=cond(P);HhPhThTp建立一個(gè)5X5矩陣，求它的行列式值、跡、秩和范數(shù)。A=[12345;23456;34567;45678;56789];det(A)trace(A)rank(A)norm(A)3.已知—29618，求A的特征值及特征向量，并分析其數(shù)學(xué)意義。A=20 5 12，求A的特征值及特征向量，并分析其數(shù)學(xué)意義?！? 8 5clearA=[-29,6,18;20,5,12;-8,8,5];[V,D]=eig(A,'nobalance')實(shí)驗(yàn)三選擇結(jié)構(gòu)程序設(shè)計(jì)1.求下列分段函數(shù)的值。廠(chǎng)x2+x—6, x<0且x豐—3y=<x2—5x+6, 0<x<10,x豐2且x豐3x2—x—1 其他用if語(yǔ)句實(shí)現(xiàn)，分別輸出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0時(shí)的y值。clearx=input('inputthevalueofx:');ifx<0&x~=3y=x*x+x-6;elseifx>=0&x<10&x~=2&x~=3y=x*x-5*x+6;elsey=x*x-x-1;endy2?輸入一個(gè)百分制成績(jī)，要求輸出成績(jī)等級(jí)A,B,C,D,E。其中90?100分為A,80?89分為B,70?79分為C,60?69分為D,60分以下為E。(1) 用switch語(yǔ)句實(shí)現(xiàn)。輸入百分制成績(jī)后要判斷該成績(jī)的合理性，對(duì)不合理的成績(jī)應(yīng)輸出錯(cuò)誤信息。clearx=input('inputscore:');ifx<0disp('Sorry,Thefailscore!')elseifx<60disp('E')elseifx<70disp('D')elseifx<80disp('C')elseifx<90disp('B')elseifx<=100disp('A')elsedisp('Sorry,Thefailscore!')end建立5X6矩陣，要求輸出矩陣第n行元素。當(dāng)n值超過(guò)矩陣的行數(shù)時(shí)，自動(dòng)轉(zhuǎn)為輸出矩陣最后一行元素，并給出出錯(cuò)信息。clearA=[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ì)兀2111 11?根據(jù)=廠(chǎng)+ + +…+—，求兀的近似值。當(dāng)n分別取100、1000、10000時(shí)，結(jié)果是多少？6 12 22 32 n2clearn=100;i=1:n;y=1./i.A2;y1=sum(y);pai=sqrt(6*y1)112?根據(jù).y二1+-+5+…（1）y<3時(shí)的最大n值。（2）于（1）的n值對(duì)應(yīng)的y值。cleary=0;n=1000;fori=1:1000y=y+1/(2*i-1);ify>=3break;endi=i+1;endi3.已知ff=11f=0<2f=13[f=f-2f +f,n>3n n-1 n-2 n-3求f1-f100中：1100最大值、最小值、各數(shù)之和。正數(shù)、零、負(fù)數(shù)的個(gè)數(shù)。f=[1,0,1];forn=4:100f(n)=f(n-1)-2*f(n-2)+f(n-3);endk=length(find(f<0))k=length(find(f>0))k=length(find(f==0))4．若兩個(gè)連續(xù)自然數(shù)的乘積減1是素?cái)?shù)，則稱(chēng)這兩個(gè)連續(xù)自然數(shù)是親密數(shù)對(duì)，該素?cái)?shù)是親密素?cái)?shù)。求[2，50]區(qū)間內(nèi)：(1)親密數(shù)對(duì)的對(duì)數(shù)；與上述親密數(shù)對(duì)對(duì)應(yīng)的所有親密素?cái)?shù)之和。clearl=0;g=2;e=1;sum=2;f(1)=1;f(2)=0;f(3)=1;fori=4:100f(i)=f(i-1)-2*f(i-2)+f(i-3);iff(1)>f(i)small=f(i)endiff(1)<f(i)big=f(i)endiff(i)<0l=l+1;endiff(i)>0g=g+1endiff(i)==0e=e+1endsum=sum+f(i);endlgesmallbigsum實(shí)驗(yàn)五函數(shù)與文件定義一個(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+2i[a,b,c,d]=susb(x)一個(gè)自然數(shù)是素?cái)?shù)，且它的各位數(shù)字位置經(jīng)過(guò)任意對(duì)換后仍為素?cái)?shù)，則稱(chēng)是絕對(duì)素?cái)?shù)。例如13是絕對(duì)素?cái)?shù)。試求所有兩位的絕對(duì)素?cái)?shù)。function[q]=apnum(p);fori=2:9n=find(rem(p,i)==0);p(n)=[];endq=10*rem(p,10)+fix(p/10);forj=2:9m=find(rem(q,j)==0);q(m)=[];end調(diào)用：clearallp=10:99;[q]=apnum(p)y=f(40)3?已知V f(30)+f(20)當(dāng)/(n)=n+l°ln(n2+5)時(shí)，求y的值。functionf=f(n)f=n+10*log(n*n+5)當(dāng)/(n)-1X2+2X3+3X4+ +nX(n+1)時(shí)，求y的值。functionf=fa(n)s=0;fori=1:n;fa(i)=i*(i+1);s=s+fa(i);ends實(shí)驗(yàn)六高層繪圖操作1.已知y1=x2,y2=cos(2x),y3=ylxy2完成下列操作：在同一坐標(biāo)系下用不同的顏色和線(xiàn)型繪制三條曲線(xiàn)。x=0:100;y1=x.*x;y2=cos(2*x);y3=y1.*y2;plot(x,y1,'k',x,y2,'b--',x,y3,'rp');axis([0,100,-10000,10000]);以子圖形式繪制3條曲線(xiàn)。x=0:pi/100:2*pi;a=0:100;y1=x.*x;y2=cos(2*x);y3=y1.*y2;b1=a.*a;b2=cos(2*a);b3=b1.*b2;subplot(2,2,1);plot(x,y1);title('y1');axis([0,2*pi,0,40]);subplot(2,2,2);plot(x,y2);title('y2');axis([0,2*pi,-1,1]);subplot(2,2,3);plot(x,y3);title('y3');axis([0,2*pi,-40,40]);subplot(2,2,4);plot(a,b1,'k',a,b2,'b--',a,b3,'rp');title('complex');axis([0,100,-10000,10000]);分別用條圖形、階梯圖、桿圖和填充圖繪制3條曲線(xiàn)。x=0:pi/10:2*pi;a=0:2:100;y1=x.*x;y2=cos(2*x);y3=y1.*y2;b1=a.*a;b2=cos(2*a);b3=b1.*b2;subplot(2,2,1);bar(x,y1);title('y1');axis([0,2*pi,0,40]);subplot(2,2,2);stairs(x,y2);title('y2');axis([0,2*pi,-1,1]);subplot(2,2,3);stem(x,y3);title('y3');axis([0,2*pi,-40,40]);subplot(2,2,4);fill(a,b1,'k',a,b2,'b--',a,b3,'rp');title('complex');axis([0,100,-10000,10000]);1分別用plot和fplot函數(shù)繪制函數(shù)y二sin—的曲線(xiàn)，分析兩曲線(xiàn)的差別。xx=0.1:0.1:20;y=sin(1./x);plot(x,y,'r');pausefplot('sin(1./x)',[0.1,20],'*')繪制函數(shù)的曲面圖和等高線(xiàn)。已知z=(x2-2x)e-x2-y2-小[x,y]=meshgrid(-1:0.5:2,1:5);z=(x.*x-2*x).*exp(-x.*x-y.*y-x.*y);surf(x,y,z)已知6<x<30,15<y<36，求不定方程2x+5y=126的整數(shù)解。x=7:29;y=16:35;[x,y]=meshgrid(x,y);%在[7,29]X[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即方程的解裁掉表現(xiàn)函數(shù)z=sin(y)cos(x)三維曲面圖中z>0.25部分。x=0:0.1:2*pi;[x,y]=meshgrid(x);z=sin(y).*cos(x);[I,J]=find(z>0.25);forii=1:length(I)z(I(ii),J(ii))=NaN;endsurf(x,y,z);實(shí)驗(yàn)七低層繪圖操作建立一個(gè)圖形窗口，使之背景顏色為紅色，并在窗口上保留原有的菜單項(xiàng)，而且在按下鼠標(biāo)器的左鍵之后顯示出LeftButtonPressed字樣。hf=figure('Color',[1,0,0],'WindowButtonDownFcn','disp(''LeftButtonPressed.'')');先利用缺省屬性繪制曲線(xiàn)y=x2e2.，然后通過(guò)圖形句柄操作來(lái)改變曲線(xiàn)的顏色、線(xiàn)型和線(xiàn)寬，并利用文字對(duì)象給曲線(xiàn)添加文字標(biāo)注y=x2e2x。x=0:0.1:10;y=x.A2.*exp(2*x);h=plot(x,y);text(9.0,81*exp(18),'\leftarrowy=xA2*eA2x','Fontsize',12);>>set(h,'Color','r','LineStyle',':','LineWidth',3);實(shí)驗(yàn)八數(shù)據(jù)處理和多項(xiàng)式計(jì)算1.將100個(gè)學(xué)生5門(mén)功課的成績(jī)存入矩陣P中，進(jìn)行如下處理：分別求每門(mén)課的最高分、最低分及相應(yīng)學(xué)生序號(hào)。分別求每門(mén)課的平均分合標(biāo)準(zhǔn)方差。5門(mén)課總分的最高分、最低分及相應(yīng)學(xué)生序號(hào)。將5門(mén)課總分按從大到小順序存入zcj中，相應(yīng)學(xué)生序號(hào)存入xsxh。A=rand(5,100);P=100*A;[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ù)值如表2所示。表2lg(x)在10個(gè)采樣點(diǎn)的函數(shù)值x1112131415161718191101Lg(x)01.04141.32221.49141.61281.70761.78531.85131.90851.95902.0043試求lg(x)的5次擬合多項(xiàng)式p(x)，并分別繪制出lg(x)和p(x)在[1,101]區(qū)間的函數(shù)曲線(xiàn)。x=[1:10:101];y=[0,1.0414,1.3222,1.4914,1.6128,1.7076,1.7853,1.8513,1.9085,1.9590,2.0043];[P,S]=polyfit(x,y,5)plot(x,y,'k*',x,polyval(P,x),'k-')3.有3個(gè)多項(xiàng)式P1(x)二X4+2X3+4x2+5,P2(x)二x+2,p?(x)二x2+2x+3,試進(jìn)行下列操作:p(x)p(x) p(x)p(x)TOC\o"1-5"\h\z求 =1 + 2 3 。1)求p(x)的根。(2)當(dāng)x取矩陣A的每一元素時(shí)，求P(x)的值。其中"-1 1.2-1.4\o"CurrentDocument"A二0.752 3.50 5 2.5當(dāng)以矩陣A為自變量時(shí)，求P(x)的值。其中A的值與(3)相同。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.*p3x=roots(p)A=[-1,1.2,-1.4;0.75,2,3.5;0,5,2.5];G=poly(A)

實(shí)驗(yàn)九符號(hào)計(jì)算基礎(chǔ)與符號(hào)微積分1.已知x=6,y=5,利用符號(hào)表達(dá)式求z=x=sym('6');y=sym('5');z=(x+1)/(sqrt(3+x)-sqrt(y))用符號(hào)方法求下列極限或?qū)?shù)。x(esinx+1)一2(etanx_1)(1)limx—0 sm3Xsymsx;axax(2)已知A=tcosxt3 dAd2Ad2Alnx'刀別求dx'dt2'dxdtsymsatx;A=[aAx,tA3;t*cos(x),log(x)];diff(A,x)diff(A,x,2)diff(diff(A,x),t)用符號(hào)方法求下列積分。1)dx1+x4+x8x=sym('x');f=1/(1+xA4+xA8);int(f)2)x1)dx1+x4+x8x=sym('x');f=1/(1+xA4+xA8);int(f)2)x2+1,dxx4+1x=sym('x');f=(xA2+1)/(xA4+1);int(f,0,inf)實(shí)驗(yàn)十級(jí)數(shù)與方程符號(hào)求解1.級(jí)數(shù)符號(hào)求和。(1)計(jì)算S=牙丄2n一1n=1n=sym('n');S=symsum(1/(2*n-1),n,1,10)n2n2xn_1之和函數(shù)，并求乙L之和。5nn=1n=1symsnx;S2=symsum((nA2)*xA(n-1),n,1,inf)S3=symsum((nA2)/5An)將ln(x)在x=1處按5次多項(xiàng)式展開(kāi)為泰勒級(jí)數(shù)。x=sym('x');f=log(x);taylor(f,x,6,1)求微分方程的符號(hào)解。?空+k2y=0dx2<y(0)=ay'(0)=ba,b,k為任意常數(shù)symsxyabk;m=dsolve('D2y+kA2*y','y(0)=a','Dy(0)=b','x')求下列方程和方程組的符號(hào)解。(1)3xex+5sinx-78.5=0x=solve('3*x*exp(x)+5*sin(x)-78.5','x')Jx2+y2一100=0(2)$*3x+5y一8=0[xy]=solve('sqrt(xA2+yA2)-100','3*x+5*y-8','x,y')x=solve('log(1+x)-5/(1+sin(x))=2','x')3)3)ln(1+x)一51+sinx(4)x2+91：x+1—1=0x=solve('xA2+9*sqrt(x+1)-1','x')實(shí)驗(yàn)十三Simulink的應(yīng)用二、實(shí)驗(yàn)內(nèi)容1.假設(shè)從實(shí)際應(yīng)用領(lǐng)域(力學(xué)、電學(xué)、生態(tài)或社會(huì))中，抽象出有初始狀態(tài)為0的二階微分方程x〃+0.2x'+0.4=°?2t(u(t)是單位階躍函數(shù)。用積分器直接構(gòu)造求解微分方程的模型exm1.mdl并仿真。2.利用傳遞函數(shù)模塊建模。仍以下面的二階微分方程為例x+°2x'+°?4=°.2u(t)實(shí)驗(yàn)十四綜合實(shí)驗(yàn)(1)2s+31?用Matlab表示傳遞函數(shù)為4S3+3S2+2S+1的系統(tǒng)。num=[23];den=[4321];

sys=tf(num,den)2(s+3)2-用Matlab表示傳遞函數(shù)為s(s+l)(s+2)的系統(tǒng)。z=-3;p=[0-1-2];k=2;y=[1y=[10]X「010X011X二001X+022X-1-2-3X131—33.用Matlab表示狀態(tài)空間表達(dá)式為u1的系統(tǒng)。x2X3A=[010;001;-1-2-3];B=[0;0;1];C=[100];D=0;Gf=ss(A,B,C,D)4.已知某控制系統(tǒng)的傳遞函數(shù)為G4.已知某控制系統(tǒng)的傳遞函數(shù)為G(s)=1s2+3s+2求Matlab描述的傳遞函數(shù)模型及零極點(diǎn)增益模型。num=[1];den=[132];Gtf=tf(num,den)Gzpk=zpk(Gtf)num=[1];den=[132];Gtf=tf(num,den)Gzpk=zpk(Gtf)用MATLAB繪制典型二階系統(tǒng)的單位階躍響應(yīng)曲線(xiàn)。step(Gtf)6.已知單位負(fù)反饋系統(tǒng)開(kāi)環(huán)傳遞函數(shù)為G(6.已知單位負(fù)反饋系統(tǒng)開(kāi)環(huán)傳遞函數(shù)為G(s)=s+2

s(s+1)(s+3)試判斷系統(tǒng)閉環(huán)穩(wěn)定性。k=1;z=[-2];p=[0-1-3];k=1;z=[-2];p=[0-1-3];G=zpk(z,p,k);Gtf=feedback(G,1)pzmap(Gtf)5(s+0.2)7.已知一單位負(fù)反饋系統(tǒng)的開(kāi)環(huán)傳遞函數(shù)為G(s)=s(s_05)(s+15)，試求單位階躍信號(hào)作為參考輸入時(shí)產(chǎn)生的穩(wěn)態(tài)誤差。k=5;z=[-0.2];p=[0z