| 查看: 1102 | 回复: 1 | ||
[求助]
二阶微分方程画三维图的问题
|
设计要求画三维图，顿时崩溃。不会编程，求大神们指教！ 其中eta1包含B0z，要求B0z是x轴。画出B0z x 和  的三维图。我编的程序是 方程文件： function dx=Yqiangpo_F(t,x,flag,eta2,P0,omega,eta1,omegan) dx=[x(2);-eta1*x(2)-omegan^2*x(1)+eta2*P0*sin(omega*t)]; 求解画图： clear clc Omega=2*pi*150; n=60*Omega/(2*pi); rho=2670; sigma0=3.63*(10^7); mu=0.34; E=71*(10^9); h=0.01; P0=5000; B0z=0:0.04:4; R=0.5; DM=E*h.^3./(12*(1-mu^2)); omega=2*pi*50; C1=-(6+2*mu)/(5+mu); C2=(1+mu)/(5+mu) A1=rho.*h.*R.^2/2+((rho.*h.*R.^2)/2-(rho*h.^3)/6)*C1+(rho*h.*R.^2/3-rho*h.^3/3)*C2+(rho*h.*R.^2/6-rho*h.^3/12)*C1^2+(rho*h.*R.^2/10-rho*h.^3/6)/C2^2+(rho.*h.*R.^2/4-5*rho*h.^3/18)*C1*C2; A2=-sigma0*h.^3.*B0z.^2.*(2*C1+4*C2+C1^2+C2^2+10*C1*C2/3)/12; A3=DM.*(32*C2^2/3+12*C1*C2+32*C2)./R.^2+rho*h.^3*Omega^2*(2*C1+4*C2+C1^2+2*C2^2+10*C1*C2/3)/12; A5=(1/2+(C1)/4+(C2)/6).*R.^2; omegan=sqrt(A3./A1); qst=(A5*P0)./(A1.*omegan.^2) s=omega./omegan; epsilon=A2./(2*omegan.*A1); eta2=A5/A1; eta1=A2/A1; beta=1./sqrt((1-s.^2).^2+(2*epsilon.*s).^2); x0=[0.00001,0.1]; t=0:0.1:10; options=odeset('RelTol',1e-6,'AbsTol',1e-6); [t,x]=ode45('Yqiangpo_F',t,x0,options,eta2,P0,omega,eta1,omegan); [B,X1]=meshgrid(B0z,x(1)); surf(B,X1,x(2)) 程序前面大部分都是系数。我想画B x1 x2 的关系。但是图中ode45开始就出错了。请问大家应该怎么改？ |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有7人回复
-
各位基础老师们，看完中央《从严治党》的决定后，对学术界更加有信心了吧
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 二阶非线性偏微分方程如何解？
已经有4人回复
- 微分方程初值问题和边值问题的区别？
已经有5人回复
- 关于根据一个系统模型建立复杂微分方程解的问题，一般都是设其解为某种形式。
已经有3人回复
- 关于二阶偏微分方程求解
已经有21人回复
- 师兄师姐们，帮小弟看看这个 Laplace变换后的二阶常微分方程怎么解啊
已经有23人回复
- 求助如何用origin绘制三维图的问题
已经有3人回复
- 非线性二阶微分方程组求解，matlab
已经有12人回复
- 求解二阶变系数常微分方程
已经有7人回复
- COMSOL例题库中（shallow-water-equation）偏微分方程求解问题
已经有3人回复
- Matlab数值求解二阶常微分方程
已经有9人回复
- 一个二阶常微分方程的求解方法？
已经有5人回复
- 如何用 数值解法 求 二阶微分方程的本征值
已经有18人回复
- 微分方程组求参数问题，求高人指教，感谢
已经有11人回复
- MATLAB中画三维双对数图形的问题
已经有10人回复
- 求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复
- 请教mathematic问题？关于解微分方程组的和画图的一些问题。
已经有19人回复
- 如何用matlab求解矩阵系数的二阶微分方程
已经有13人回复
- 如何求解描述振动的二阶微分方程
已经有21人回复
- 关于函数作图和解微分方程的软件两个问题
已经有4人回复
- 【求助】先解常微分方程 然后画图【已完结】
已经有11人回复
2楼2015-04-21 15:39:24
50