2019年matalab實(shí)踐任務(wù)與書寫規(guī)范11034

1、說明:報(bào)告嚴(yán)格按照以下各項(xiàng)獨(dú)立完成，若有抄襲，以不及格處理，書寫工整浦江學(xué)院實(shí)習(xí)報(bào)告（2011-2012學(xué)年第二學(xué)期）課程名稱：Matlab實(shí)踐班級(jí)：浦通信1103-4學(xué)號(hào)：01姓名：張三 指導(dǎo)教師：王琛譚紅一、目的和要求MATLA映踐是通信工程專業(yè)技術(shù)基礎(chǔ)類別的必修環(huán)節(jié)，也是通信工程專業(yè)的重要實(shí) 踐環(huán)節(jié)。本次實(shí)踐旨在培養(yǎng)通信工程專業(yè)的學(xué)生理論聯(lián)系實(shí)際的能力，為通信工程的素 質(zhì)教育提供Matlab科學(xué)計(jì)算工具平臺(tái)。將學(xué)習(xí) Matlab這一現(xiàn)代的科學(xué)計(jì)算和系統(tǒng)仿真 語言的基本編程思想和方法，并利用Matlab對(duì)所學(xué)基礎(chǔ)課程進(jìn)行上機(jī)模擬實(shí)驗(yàn)和數(shù)值計(jì) 算，從而通過Matlab編程實(shí)驗(yàn)來驗(yàn)證和鞏固所學(xué)

2、的數(shù)學(xué)和工程理論，為進(jìn)一步學(xué)習(xí)通信仿真技術(shù)等課程打下扎實(shí)的計(jì)算機(jī)應(yīng)用能力的基礎(chǔ)。具體要求如下：1）熟悉Matlab運(yùn)行環(huán)境，了解Matlab的搜索路徑、菜單欄和工具欄；掌握 Matlab 各種表達(dá)式的書寫規(guī)則；熟悉矩陣的產(chǎn)生、操作和運(yùn)算，利用數(shù)學(xué)函數(shù)和邏輯函數(shù)實(shí)現(xiàn) 簡(jiǎn)單功能。2）熟悉曲線和坐標(biāo)軸的控制；了解高級(jí)圖形的繪制及圖形注釋。3）掌握函數(shù)文件和腳本文件的運(yùn)行方法，函數(shù)中流程控制的設(shè)計(jì)，用戶交互輸入的方式和程序調(diào)試的步驟。4）熟悉線性方程組的求法和矩陣函數(shù)；掌握 Matlab在電路分析和物理等方面的應(yīng) 用。5）能夠利用基本指令實(shí)現(xiàn)簡(jiǎn)單的設(shè)計(jì)工作。二、內(nèi)容1） Matlab的軟件環(huán)境Matl

3、ab的啟動(dòng)和退出，在命令窗口中計(jì)算數(shù)學(xué)表達(dá)式、輸入簡(jiǎn)單矩陣，查看工作空 問中的變量，命令窗口顯示清除，語句、變量及數(shù)據(jù)的表示，搜索路徑設(shè)置，簡(jiǎn)單程序 的編寫等。2） Matlab數(shù)值計(jì)算矩陣生成及計(jì)算、數(shù)組運(yùn)算、線性方程組求解、微積分問題求解等。3） Matlab程序設(shè)計(jì)建立和執(zhí)行M文件；程序選擇結(jié)構(gòu)、分支結(jié)構(gòu)、循環(huán)結(jié)構(gòu)、；函數(shù)文件的實(shí)現(xiàn)，定義 和調(diào)用Matlab函數(shù)；直流穩(wěn)態(tài)電路、正弦穩(wěn)態(tài)電路和一階動(dòng)態(tài)電路的求解；雙縫干涉和 衍射的模擬。4） Matlab繪圖功能二維圖形的繪制、三維圖形的繪制、圖形窗口的建立、圖形動(dòng)畫的實(shí)現(xiàn)等。5） Matlab符號(hào)計(jì)算與字符申利用符號(hào)計(jì)算實(shí)現(xiàn)“代數(shù)演算”，

4、字符串的生成和比較等。三、步驟基本操作1 .試分別生成3階的單位陣和4階的零矩陣eye zeroseye(3)>>ans =1 00010001zeros(4) ans =00000000000000002 .生產(chǎn)列向量 x=1,3, 5, 7, 9,，15冒號(hào)x=1:2:15 x =135791113153 .生成以x的前5個(gè)元素為對(duì)角線的矩陣 AdiagA=diag(x(1:5)10000030000050000070000094 .生成一個(gè)與A同階的隨機(jī)矩陣Brand sizeB=rand(size(A)B =5 .生成由B與其同階魔方矩陣點(diǎn)乘得到的矩陣DMagicD=B.*

5、magicze(B)D =6 .生成由D的第2、4、5行和第4、1列組成的子矩陣 EE=D(2 4 5,4 1)E =2967. A 205881812 ,求A的特征值和特征向量。5eigA=-29 6 18;20 5 12;-8 8 5V,D=eig(A)A =-2961820512-885V =D =00>>符號(hào)運(yùn)算(symS) n1、計(jì)算 L lim lg(x h) 1g(x)和 M lim 1 - h 0 hn nsyms , limitsyms h xf=(log(x+h)-log(x)/hlimit(f,h,0)>>f =(log(x+h)-log(x)/h

6、ans =1/xsyms n xm=(1-x/n)Anlimit(m,n,inf)m =(1-x/n)Anans =exp(-x)2、化簡(jiǎn)函數(shù) f(x) 3 139 12 8,x x xsyms , simplifysyms c alpha betaf=(1/xA3+6/xA2+12/x+8)A(1/3)y1=simplify(f)g1=simple(f)g2=simple(g1)f =y1 =(2*x+1)A3/xA3)A(1/3)gi =(2*x+1)/xg2 =2+1/x3、指出下面各條語句的輸出結(jié)果f=sym('2*u');subs(f,'u',2)f2

7、=subs(f,'u,u+2,)a=3;subs(f2,'u',a+2)subs(f2,'u','a+2')syms x yf3=subs(f,'u',x+y)subs(f3,x,y,1,2)f =2*uans =4f2 =2*(u+2)a =3ans =14ans =2*(a+2)+2)f3 =2*x+2*yans =6程序的編寫1、隨機(jī)投擲骰子，驗(yàn)證各點(diǎn)出現(xiàn)的概率是否為1/6 。 randperm switch casen=6000m1=0;m2=0;m3=0;m4=0;m5=0;m6=0for i=1:nx=rand

8、perm(6)y=x(1)switch ycase 1m1=m1+1;case 2m2=m2+1;case 3m3=m3+1;case 4m4=m4+1;case 5m5=m5+1;otherwisem6=m6+1;endendb1=m1/n,b2=m2/n,b3=m3/n,b4=m4/n,b5=m5/n,b6=m6/nb1 =b2 =b3 =b4 =b5 =b6 =2、兩個(gè)連續(xù)自然數(shù)的乘積減一是素?cái)?shù)，則稱這兩個(gè)連續(xù)自然數(shù)是一對(duì)親密數(shù)素?cái)?shù),求2, 50之間的親密素?cái)?shù)的對(duì)數(shù)for if3、某公司員工的工資計(jì)算方法如下:(1)工作時(shí)間超過120小時(shí)，超過部分加發(fā)15%(2)工作時(shí)間低于60小時(shí)，扣

9、發(fā)15%(3)其余按每小時(shí)80元計(jì)發(fā)。編程實(shí)現(xiàn)輸入工作時(shí)間，輸出應(yīng)發(fā)工資。gh=input('請(qǐng)輸入工號(hào)gh'),gs=input('請(qǐng)輸入工時(shí)gs');if gs>120%gz=120*80+(gs-120)*(1+*80elseif gs<60%gz=gs*80-700else gz=gs*80%end1112344、求解線性方程組111345111456inv normA=1/2 1/3 1/4;1/3 1/4 1/5;1/4 1/5 1/6;b1= '%x1=inv(A)*b1%求工時(shí)超過120小時(shí)的工資求工時(shí)低于60小時(shí)的工資求工時(shí)

10、為其他值的工資X10.95x20.67,并計(jì)算系數(shù)矩陣A的條件數(shù)。X30.52產(chǎn)生矩陣產(chǎn)生矩陣Ab的秩求方程的解C=norm(A)結(jié)果分析x1 = 6/53/53/5C = 9473/75、模擬楊氏雙縫干涉:兩相干光源到接收屏的距離r1D22D21d 2 2y -,相位差22 r2-1-,光強(qiáng) I 4I0 cos2 一 。取25 107, d 2 103, D=1, I0=iplotlinspace for sqrt sum image口 1234clear alllam=500*10A(-9);% 輸入波長(zhǎng)a=2*10A(-3);D=1;ym=5*lam*D/a;xs=ym;%設(shè)定光屏的范圍

11、n=101;ys=linspace(-ym,ym,n);%把光屏的 y 方向分成 101 點(diǎn)for i=1:nr1=sqrt(ys(i)-a/2).A2+DA2);r2=sqrt(ys(i)+a/2)A2+DA2);phi=2*pi*(r2-r1)/lam;B(i,:)=4*cos(phi/2).A2;endN=255;%確定用的灰度等級(jí)為255級(jí)Br=(B/*N;%使最大光強(qiáng)對(duì)應(yīng)于最大灰度級(jí)(白色)subplot(1,2,1)image(xs,ys,Br);%畫干涉條紋colormap(gray(N);subplot(1,2,2)plot(B,ys)%畫出光強(qiáng)變化曲線x(0) 0y(0) 1

12、 , t 0,12的數(shù)值解z(0) 1(一 一 x(t) yz丁 ，一、 一一 , 一,6、求解方程組 y (t) xz'一一 _ .z (t)0.51xyode45 functionx,y,z=dsolve(,Dx=y*z,Dy=-x*z,Dz=*x*y',x(0)=0, 'y(0)=1','z(0)=1','t')電路分析1、電路如圖，已知R=5歐姆，wL=3歐姆，wC=ffln, Uc 10 30 ,求Ir, %,Il , Ul , Us ,并畫出相量圖abs angle2、下圖中，開關(guān)S在位置1已很久，t=0時(shí)開關(guān)合向位置

13、2,求換路后的電流i和若在端口 ab接負(fù)載ZL, ZL為何值時(shí)它能獲得最大功率？最大功率為多大?繪圖1、繪制w x2 y2 z2的圖形。meshgrid surfacex,y,z=meshgrid(-1:1,-1:1,-1:1);w=x.A2+y.A2+z.A2;isosurface(x,y,z,w)1 、2、繪制雙葉雙曲面，a=b=3, c=5o< P< 2fx = 2 tan©cos，y - ft - tansin£?z = c * see中Ezsurfezsurf('3*tan(u)*cos(v)','3.*tan(u)*sin(v

14、)','5.*sec(u)',-pi./2,3*pi./2,0,2*pi) axis equalgrid onaxis squarexlabel('x 軸');ylabel('y 軸');zlabel('z 軸')；title('雙葉雙曲面,)雙葉雙曲面3、繪制軸截面的曲邊為拋物線的螺線（取 a=2, b=2, c=1/3, 0 t 50x-a-l' cost，y = b*tsintEzplotezplot3（'2*t*cos','2*t*sin'， '1/3*t*t', 0,50）;grid onaxis squarexlabel（'x 軸'）;ylabel（'y 軸'）;zlabel（'z 軸'）；title（'軸截面的曲邊為拋物線的螺線,）4、實(shí)現(xiàn)霓虹燈閃爍效果，圖形用 z 17x2 16| x| y 17y2 225inspace meshgrid pcolor shading interp pa