北京交通大學(xué)matlab復(fù)習(xí)資料

1、北京交通大學(xué)matlab復(fù)習(xí)資料P1051. 主程序：s1,kk=quad(f1,0,2);s2,kk=quad(f2,2,4);s=s1+s2 函數(shù)：function y=f1(x)y=log(x.2);function y=f2(x)y=log(16)./(2+sin(x+1).*pi);f1=inline(log(x.2),x);f2=inline(log(16)./(2+sin(x+1).*pi),x);s1,kk=quad(f1,0,2);s2,kk=quad(f2,2,4);s=s1+s2 3. syms x a taus=int(log(x.2),0.5,tau)solve(s=

2、10)4. 求解方程x5+6x4-3x2=10的5個根，并將其位置用五角星符號標(biāo)記在復(fù)平面上，要求橫縱坐標(biāo)軸刻度等長，注明虛軸和實軸，在title位置上寫出方程。clc;close all;clear all; factor=1 6 0 -3 0 -10; ro=roots(factor);disp(ro);scatter(real(ro),imag(ro),p);axis equal;title(x5+6x4-3x2=10);xlabel(real axis);ylabel(imagine axis(j); 6. 某班同學(xué)成績已經(jīng)存放在矩陣A中，每行為某一位同學(xué)的數(shù)據(jù)，第1列為學(xué)號，第2列至

3、第4列為其三門課程的成績，試編程按照3門課平均成績由小到大的順序重新排列成績表，并放在矩陣B中。clc;clear all;close all;A=,95,73,88; p873. 在同一坐標(biāo)系下繪制t2，- t2，t2sint在t0,2pi內(nèi)的曲線圖；選擇合適的“曲線線型”，“曲線顏色”、“標(biāo)記符號”選項，并在圖形上加注坐標(biāo)軸名和圖名。 close all;clear all;clc; t=linspace(0,2*pi,1024); y1=t.2; y2=-t.2; y3=(t.2).*sin(t); figure(2); plot(t,y1,:y,t,y2,-g,t,y3,-r); le

4、gend(y=t2,y=-t2,y=t2*sin(t); title(三個函數(shù)圖形);xlabel(t);ylabel(f(t); 4. theta=linspace(0,16*pi,2048); subplot(2,2,1);polar(theta,1.0013*theta.2); subplot(2,2,2);polar(theta,cos(3.5*theta); subplot(2,2,3);polar(theta,sin(theta)./theta); subplot(2,2,4);polar(theta,1-(cos(7*theta).3); 5. ttl=主視圖,左視圖,俯視圖,三

5、維圖;angle=0,0,-90,0,0 90,-37.5,30;for i=1:4 subplot(2,2,i); x=0:0.01:1;y=0:0.02:2;x,y=meshgrid(x,y);z=1./sqrt(1-x.2)+y.*2)+1./sqrt(1+x.2)+y.2);mesh(x,y,z) view(anglei);title(ttli);end6. 用鼠標(biāo)左鍵在圖形窗口上取5個點，在每個點的位置處寫出一個字符串來顯式鼠標(biāo)點的橫坐標(biāo)值，然后將這些點連成折線。 axis(0,5,0,5); hold on; box on; x=;y=; gtext(); for i=1:5 x1

6、,y1,button=ginput(1); if(button=1) break; end plot(x1,y1,o); x=x,x1;y=y,y1; text(x1,y1,num2str(x1); line(x,y); end hold off老師： figure (1)u,v,b=ginput(5);line(u,v);for i=1:5;str=(,num2str(u(i),);text(u(i),v(i),str);end油桶clc;clear;syms x y h;w=ezplot(x2/(2.52)+(y-2)2/4-1=0); %繪制橢圓.set(w,Color,k); hold

7、 on;axis(-2.55,2.55,0,4.05); %設(shè)置坐標(biāo)值.set(gca,XTick,XColor,w); %將坐標(biāo)隱去.set(gca,YTick,YColor,w);line(0,0,0,4,Color,k,LineStyle,-) %畫出中間刻度線.hold on;s=int(100-25*(y-2).2)/16).0.5,0,h); %用int函數(shù)求不定積分. for i=1:25 H(i)=solve(s-0.3125*i,h); %求得各個刻度的高度. end for i=1:25 if mod(i,4)=0 plot(-0.3,0.3,H(i),H(i),b) ;

8、%每四格為一長刻度. else plot(-0.1,0.1,H(i),H(i),k) ; end end%標(biāo)注長刻度分別為20,40,60,80,100,120.text(0.4,0.86,20); text(0.4,1.42,40);text(0.4,1.93,60);text(0.4,2.43,80);text(0.35,2.97,100);text(0.35,3.64,120);用循環(huán)語句形成一個有20個分量的數(shù)組，滿足fibonacci規(guī)則，令第K+2個元素滿足a（k+2）=a（k）+a（k+1）a(1)=1;a(2)=1;for k=1:18a(k+2)=a(k)+a(k+1);en

9、da在同一圖形窗口上分成四個繪圖區(qū)域，試分別繪制主視圖，左視圖，俯視圖和默認視角的三位表面圖。ttl=主視圖,左視圖,俯視圖,三維圖;angle=0,0,-90,0,0 90,-37.5,30;for i=1:4 subplot(2,2,i); ezmesh(關(guān)于x,y,z的關(guān)系表達式,取值范圍); view(anglei);title(ttli);end ttl=主視圖,左視圖,俯視圖,三維圖;angle=0,0,-90,0,0 90,-37.5,30;for i=1:4 subplot(2,2,i); ezmesh(1/(1-X)2+Y2)0.5+1/(1+X)2+Y2)0.5,-2 2)

10、; view(anglei);title(ttli);End試用文字和數(shù)學(xué)公式描述下面指令要完成的計算任務(wù)roots(11 0 -6 2 0 -3)roots為求多項式的根，其向量表達多項式的系數(shù)，對應(yīng)降冪排列，那么就有本題對應(yīng)求解11x5-6x3+2x2-3的根。試舉例說明matlab語言中三位數(shù)組是如何定義的。在Matlab中習(xí)慣性的會將二維數(shù)組中的第一維稱為“行”第二維稱為“列”，而對于三維數(shù)組的第三位則是習(xí)慣性的稱為“頁”。在Matlab中將三維及三維以上的數(shù)組統(tǒng)稱為高維數(shù)組，三維數(shù)組也是高級運算的基礎(chǔ)。舉例如課本31頁。在程序段中插入語句A=zeros(10000,1);能顯著縮短程

11、序運行時間。因為插入本語句后，在循環(huán)語句前預(yù)先為變量A分配了足夠的單元，循環(huán)時只需向各元素內(nèi)填數(shù)，如果不插入的話則未預(yù)先定義A的大小，每次執(zhí)行一次循環(huán)體便要擴充一次y的長度，增加了內(nèi)部工作量。Matlab語言的基本編程單元為 復(fù)數(shù)矩陣攝氏度轉(zhuǎn)華氏度clc;clear; C=input(c=); f=32+C*1.8; disp(sprintf(f=%2.2f,f); 華氏度轉(zhuǎn)攝氏度clc;clear; f=input(f=); c=(f-32)/1.8; disp(sprintf(c=%2.2f,c); 實例二： 數(shù)據(jù)的多項式曲線擬合例： 已知某壓力傳感器的標(biāo)定數(shù)據(jù)見下表，p為電壓值，u為電壓

12、值。p0.01.12.12.84.25.06.16.98.19.09.9u1011131417182224293439試用多項式擬合其特性函數(shù)，求出和，把擬合曲線和各個標(biāo)定點畫在同一幅圖上。解：將壓力視為自變量，電壓視為函數(shù)，編程如下：p=0,1.1,2.1,2.8,4.2,5,6.1,6.9,8.1,9,9.9;u=10,11,13,14,17,18,22,24,29,34,39;A=polyfit(p,u,3);a=A(1),b=A(2),c=A(3),d=A(4),p1=0:0.01:10;u1=polyval(A,p1);plot(p1,u1,p,u,o)1.如何求未保存數(shù)據(jù)(上次ma

13、tlab執(zhí)行了 save x )：執(zhí)行l(wèi)oad x y即可得到y(tǒng)的數(shù)據(jù)2.畫視圖 P87頁上3.roots的作用相當(dāng)于求11s5+0s4-6s3+2s2+0s-3=0的根。4.三維數(shù)組的建立 P31下5.A=zeros（，1）的意義預(yù)先定維可以節(jié)約時間。1.求定積分函數(shù)程序1：函數(shù)名為jifen1.mfunction y=jifen1（x）y=exp(pi-x)函數(shù)程序2：函數(shù)名為 jifen2.mfunction y=jifen2（x）y=log10(16)/(2+sin(x+1)*pi)主程序：y1=quad(jifen1,0,3.5);y2=quad(jifen2,3.5,7);y=y1+y12.求解微分方程函數(shù)程序：函數(shù)名difen.mfunction xdot=difeq(t,x)xdot=zeros(4,1)xdot(1)=x(2)xdot(2)=2*x(4)+x(1)xdot(3)=x(4)xdot(4)=-2*x(2)+x(3)主程序：t0=0;tf=20;x0=1.1 0 0 -1;t,x=ode23(difeq,t0,tf,x0)xx=x(:,3)yy=x(:,3)plot(yy,xx)3.求采樣點的幅頻、相頻、條形、波形圖（數(shù)據(jù)y和h存儲在signaldata.mat文件中）。load signal