matlab画波动图像,【基于Matlab的波动方程的可视化实现最终版材料】

基于Matlab的波动方程的可视化实现(最终版)

《基于Matlab的波动方程的可视化实现.doc》由会员分享，可免费在线阅读全文，更多与《基于Matlab的波动方程的可视化实现(最终版)》相关文档资源请在帮帮文库(www.woc88.com)数亿文档库存里搜索。

1、()解析解可以用以下程序作动画，程序中取,,,rvka；clearr=；v=；k=；a=；theta=linsace(,*i,)；rho=::；[Th,Rh]=meshgrid(theta,rho)；[X,Y]=olcart(Th,Rh)；rh=sqrt(X^+Y^)；th=atan(YX)；fort=::u=real(v*r^*(rh^+i*krh)*cos(th)*ex(k*(rh+*t)*i))；surf(X,Y,u)view(,)ause()end图是动画中的几幅画面，这是由球面向外传播的球面波。虽然从远场近似解的表达式可以看出远场的极化现象，也就是声波传播具有明显的两极的方向性，但是从近场的表达式是看不出来的，将近场的解析解画出的图形以后，它的极化现象其实更显著。图解析解的表面图图也可以用等值线来作动画演示，只需将指令中下列两句surf(X,Y,u)view(,)改为contour(X,Y,u)axis([])axissquare即可。图就是动画中的几幅画面。图解析解等值线图四级声源的研究半径为r的球面径向速度分布为(cos。

2、选择Color,Height(Dlot)和Animation，在外为了看得更清，把函数值适当放大，做法在Proerty栏目下与Height(Dlot)对齐的位置，选择Userentry栏目下，与它对齐的位置在空白栏中填入*u最后单击按钮Plot即可。这里画的过柱轴的截面的运动，不难想象，面上每个点都应该围绕平衡位置来回振动，这真是动画图所表现的图像。图柱体内的振动柱体外的振动问题研究半径为的长圆柱面，其径向速度分布为cosvvt，试求解这个长圆柱面在空气中辐射出去的声场的中的速度势。设远小于声波的波长。[]所求的速度势满足二维波动方程，取平面极坐标系，极点在柱轴上，则定解问题是()costtuauuvt()问题的解析解是下式的实部：()()itvuiHea()在远场区即大的区域，渐近解为sin()auvta()这是振幅按减小的柱面波。下面PDETOOL求解数值解。画两个同心圆C和C，圆心都在原点，半径分别为和两个值。C代表柱体的横切面,C表示求解区域的外边界。问题的求解区。

3、最终的解为()()(,,)cos()()sinsin()kAckkyuxytatxkbccb()从这个解中可以看出，矩形膜的本征函数和本征值分别为()()sinsin()()mmmnnxmyXxYycbmnbc()圆膜的振动边缘固定，半径为的圆形膜，初始形状是旋转抛物面()()utu，初始速度为零，求膜的振动情况。[]定解问题是()()()ttuauuutu()它的解是()()()()(,)()cos[]()nnnnnxxuxtuJatxJx()下面直接用偏微分方程工具箱来解决这个问题。计算中可取,,au。先画一个半径为圆心在原点圆，边界条件为周边固定，所以不必改变默认设置。方程取Hyerbolic(双曲型)，方程的系数取,,,cafd。方程求解的时间范围可取::，初速取零而在初始位移一栏中输入*(x^y^)作图时选择contour,Height(Dlot),Animat。

4、域是由这两个同心圆组成的环形区域CC。在外边界C上，取狄里克利边界条件u表示无穷远函数为零，即,hr，在内边界C上去纽曼边界条件，为cos,gtq。方程取Hyerbolic(双曲型)，在对话框中取,,,cafd。方程求解的时间范围可取:，初始位移和初速都取零。作图时选择contour和Animation，最后单击Plot即可。图是画出的函数在某个时刻的等值线，在动画中，犹如水波一样往外传播。图向外传播的柱面波偶极声源的研究半径为r的球面，径向速度分布为coscosvvt，试求解这球面所发射的稳恒声振动的速度势u，设r远小于声波的波长。[]用球坐标，极点取在球心，定解问题是()(cos)ttitruauurrvPe()在球面的边界条件是cost即为itRe，上面写成了ite，这要求在计算结果中也取实部。问题的解析解是()()(cos)iratavruierar()上式的实部就是所要求的解。在远场取渐近公式近似，并取实部，得到的解为()(cos)iratavrue。

5、ion,最后单击按钮Plot即可。所得图形如图所示。图圆形膜的振动用Matlab软件研究三维振动问题柱体内的振动研究匀质圆柱，半径为，高L，上下底面固定，侧面自由，初始位移为零，初始速度为uz，求柱体内各处的振动情况。[]定解问题是()()()()()tttuauuuzuzLututuz()问题的解析解是()()()()()()()()()sin()()sin()nnnnnnuLuaxJxxLxxatzJLL()下面用PDETOOL求数值解，以x表示柱体的高度，y表示柱体宽度，画一个柱体的纵切面。为了较好的演示效果，求解的区域是,xy。在柱体的上下底，取狄里克利边界条件即在对话框中取,hr，在侧边界上诺伊曼边界条件，即在对话框中取为,gq。方程取Hyerbolic(双曲型)，在对话框中取,,,cafd。方程求解的时间范围可取::，在初始位移取零，初速y*y*x，为了获得更高的精度，将区域划分网络以后，要作两次细分。在作图对话框中，。

6、选项Otion中，Animationrate一栏中填入，Numberofreeat一栏填入。另外在Colorma一栏中选ink,在动画时间(Timeforlot)一栏中填入最后单击按钮Plot即可。可以看到，是初始位移形成的波在膜上来回传播，如图所示。图矩形膜的振动这个问题是有解析解的，解法如下，将定解问题分离变量得三个方程为''()()()()()XxXxXxXxc()''''()()()()()YyYyYyYybTtaT()分别求得()sin,(),,,()sin,(),,,()cos()sin()mmmmnmmnmnmnmnmnmymYymbbmxnXxnbcTtCatDat定解问题的解的形式,,(,,)()()()(cos()sin())sinsinmnmmmnmnmnmnmnmnuxytTtXxYynxmyCatDatcb()由初始条件定出,,,,,,,,,,(),,,mnmnnDmnCmAcnkkCnnkk。