matlab物理例子,matlab在物理上的应用举例.doc

matlab在物理上的应用举例.doc

1单列波单列波T0000110A 振幅A W 频率W A A YASINWTAPLOTT,YPAUSE1,SOUNDYYLABEL Y ,XLABEL T 2光的单缝衍射现象LAMBDA500E9A A 可取02E3,1E3,2E3三种情况Z1YMAX3LAMBDAZ/ANY51YSLINSPACEYMAX,YMAX,NYNPOINTS51YPOINTLINSPACEA/2,A/2,NPOINTSFORJ1NYLSQRTYSJYPOINT2Z2PHI2PILZ/LAMBDASUMCOSSUMCOSPHISUMSINSUMSINPHIBJSUMCOS2SUMSIN2/NPOINTS2ENDCLF,PLOTYS,B, ,YS,BGRID3用毕奥－沙伐尔定律计算电流环产生的磁场MU04PI1E7I050RH1C0MU0/4PII0NGX21NGY21XLINSPACERH,RH,NGXYLINSPACE3,3,20YXNH20THETA0LINSPACE0,2PI,NH1THETA1THETA01NHY1RHCOSTHETA1Z1RHSINTHETA1THETA2THETA02NH1Y2RHCOSTHETA2Z2RHSINTHETA2DLX0DLYY2Y1DLZZ2Z1XC0YCY2Y1/2ZCZ2Z1/2FORI1NGYFORJ1NGXRXXJXCRYYIYCRZ0ZCR3SQRTRX2RY2RZ23DLXR_XDLYRZDLZRYDLXR_YDLYRXDLXRZBXI,JSUMC0DLXR_X/R3BYI,JSUMC0DLXR_Y/R3ENDENDCLFQUIVERX,Y,BX,BY4多普勒效应X0500V50Y020C330W1000T0000130RSQRTX0VT2Y02T1TR/CUSINWTSIN11WTU1SINWT1SIN11WT1SOUNDUPAUSE5SOUNDU15亥姆霍兹线圈CLEARALLMU04PI1E7I050RH1C0MU0/4PII0NGX21NGY21XLINSPACERH,RH,NGXYLINSPACERH,RH,NGYNH20THETA0LINSPACE0,2PI,NH1THETA1THETA01NHY1RHCOSTHETA1Z1RHSINTHETA1THETA2THETA02NH1Y2RHCOSTHETA2Z2RHSINTHETA2DLX0DLYY2Y1DLZZ2Z1XC0YCY2Y1/2ZCZ2Z1/2FORI1NGYFORJ1NGXRXXJXCRYYIYCRZ0ZCR3SQRTRX2RY2RZ23DLXR_XDLYRZDLZRYDLXR_YDLYRXDLXRZBXI,JSUMC0DLXR_X/R3BYI,JSUMC0DLXR_Y/R3ENDENDBAXBX,1121BX,111BAYBY,1121BY,111SUBPLOT1,2,1MESHX1121,Y,BAXXLABEL X YLABEL Y SUBPLOT1,2,2,PLOTY,BAX,GRID,XLABEL Y YLABEL BX ６．库仑引力CLEARALLN 电荷数目N FORIC1NFPRINTF \N对电荷G\N ,ICRC 电荷位置XY(米) XICRC1YICRC2QIC 输入电荷量(库仑) ENDE0885E12C01/4PIE0FORIC1NFX00FY00FORJC1NIFICJCXIJXICXJCYIJYICYJCRIJSQRTXIJ2YIJ2FXFXC0QICXIJ/RIJ3FYFYC0QICYIJ/RIJ3ENDENDFPRINTF 其它电荷作用在电荷G上的合力为\N ,ICFPRINTF X－分量GN\N ,FXFPRINTF Y－分量GN\N ,FYEND７．李萨如图形LISARUMSYMSTA1A2W1W2XCOSW1TA1YSINW2TA2A1 A1 A2 A2 W1 W1 W2 W2 TF10NS1000TLINSPACE0,TF,NSDTTF/NS1分NS个点,求出时间增量DTXPLOTXYPLOTY计算NS个点的位置XT,YTFIGUREGCFSUBPLOT1,2,1,FORI1750PLOTYPLOT1I,XPLOT1I画点的轨迹图AXIS EQUAL GRID使两轴比例相同PAUSE001END８．耦合振子M12M22K116K24X010XD066TF10MM1,00,M2KK1K2,K2K2,K1K2U,S�U,LEIGK,MTLINSPACE0,TF,101XZEROS2,101FORS12ALFASQRTU,S MU,SU,SU,S/ALFAWJSQRTLJ,JXTU,JU,J MX0COSWJTU,J MXD0/WJSINWJTXXXTENDFORR12X1,X1SUBPLOT2,1,RPLOTTXR,,GRIDXLABEL XXX YLABEL YYY ,NUM2STRREND９．拍频T0000110A1 1 W1 1 A2 2 W2 2 Y1A1SINW1TY2A2SINW2TYY1Y2SUBPLOT3,1,1,PLOTT,Y1,YLABEL Y1 SUBPLOT3,1,2,PLOTT,Y2,YLABEL Y2 SUBPLOT3,1,3,PLOTT,Y,YLABEL Y ,XLABEL T PAUSE,SOUNDY1PAUSE5,SOUNDY2PAUSE5,SOUNDY,PAUSESUBPLOT1,1,1１０。物块下滑M1 M1 M2 M2 THETA THETA度 THETATHETAPI/180G981AM1COSTHETA,M1,SINTHETA,0M1SINTHETA,0,COSTHETA,00,M2,SINTHETA,00,0,COSTHETA,1B0,M1G,0,M2G XA\BA1X1,A2X2,N1X3,N2X4１１。循环XUNHUANR831GAMA14注由于在MATLAB中无拉丁文,所以用GAMA代替NMOLES05P11E5V10012WTOTAL0QTOTAL0IPOINT1NCURVE100PPLOTP1VPLOTV1变量和图形初始化，输入气体的摩尔数NMOLES，初始压力P(1)，初始容积V(1)，气体常数R＝8314；给定起始总功WTOTAL0点序号IPOINT1画等温线用的点数NCURVE100P－V图第一点坐标PPLOTP1VPLOTV1T1P1V1/NMOLESR算出初始温度为了进入循环，先要设两个不相等的PATHTYPE和QUITTYPE值QUITTYPE5PATHTYPE0WHILEPATHTYPEQUITTYPE在菜单上选择‘退出’之前不断循环，选择路径类型或退出IPOINTIPOINT1下一点FPRINTF 对过程＃％G\N ,IPOINT1PATHTYPEMENUSPRINTF 过程％G选择下一路径 ,IPOINT1, 等压 , 等容 , 等温 , 绝热 , 退出 图形界面菜单生成语句SWITC