| 查看: 649 | 回复: 4 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】画三维图【已解决】
|
|||
|
function equ1=funb(t,y,flag,A,B,C,D,E,F,K,M) equ1=[y(2);(-E*y(2)-F*y(1)+((A*t+B)./M)*(cosh(C*D*t.^2./2)-cos(C*D*t.^2./2)-K*(sinh(C*D*t.^2./2)-sin(C*D*t.^2./2))))]; function ctn A=1;B=1;C=1;D=1;E=0.1;F=1;K=1;M=2; y0=[0;0]; tspan=[0:0.1:40]; [t,y]=ode45('funb',tspan,y0,[],A,B,C,D,E,F,K,M); 以上我是求的一个微分方程。 接下来，我要画一个三维图，Z(x,t)=Y(t)*W(x) W(x)表达式已知，现在我的问题是我不知道该如何把解微分方程得到的数值解Y(t)加入到画图程序中，我试着用过surf and plot3但总是出错，，，，望高手不吝赐教！ W(x)=(cosh(D*x)-cos(D.*x)-K.*(sinh(D.*x)-sin(D.*x))) [ Last edited by nono2009 on 2010-3-30 at 09:21 ] |
回复此楼» 猜你喜欢
应该是93bebmhtak前后十一个字符比较关键
已经有4人回复
-
这样的filecode谁见过
已经有16人回复
-
关于代码变化问题，想知道的进来
已经有16人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
nono2009(金币+2):鼓励应助！ 2010-03-26 18:25
小木虫(金币+0.5):给个红包，谢谢回帖交流
nono2009(金币+2):鼓励应助！ 2010-03-26 18:25
|
function ctn global A B C D E F K M A=1; B=1; C=1; D=1; E=0.1; F=1; K=1; M=2; y0=[0;0]; tspan=0:0.1:40; [t,y]=ode45(@funb,tspan,y0); x=0:0.1:40; Wx=(cosh(D*x)-cos(D.*x)-K.*(sinh(D.*x)-sin(D.*x))); Z1=y(:,1)*Wx; Z2=y(:,2)*Wx; [T,X]=meshgrid(0:0.1:40,0:0.1:40); figure(1) mesh(T,X,Z1) title('y1') figure(2) mesh(T,X,Z2) title('y2') function equ1=funb(t,y) global A B C D E F K M equ1=[y(2); (-E*y(2)-F*y(1)+((A*t+B)./M)*(cosh(C*D*t.^2./2)-cos(C*D*t.^2./2)-K*(sinh(C*D*t.^2./2)-sin(C*D*t.^2./2))))]; |
2楼2010-03-26 16:58:06
|
谢谢，我后来也编出来了，用的是个循环。但是图不如你的好看。。。 function ctn A=1;B=1;C=1;D=1;E=0.1;F=1;K=1;M=2; y0=[0;0]; tspan=[0:1:40]; [t,y]=ode45('funb',tspan,y0,[],A,B,C,D,E,F,K,M); x=0:0.025:1; a=length(t); b=length(x); k=1; for i=1:a for j=1:b z(j,i)=(cosh(D*x(j))-cos(D.*x(j))-K.*(sinh(D.*x(j))-sin(D.*x(j)))).*y(i,1); end end mesh(t,x,z); xlabel('t'),ylabel('x'),zlabel('z') colorbar function equ1=funb(t,y,flag,A,B,C,D,E,F,K,M) equ1=[y(2);(-E*y(2)-F*y(1)+((A*t+B)./M)*(cosh(C*D*t.^2./2)-cos(C*D*t.^2./2)-K*(sinh(C*D*t.^2./2)-sin(C*D*t.^2./2))))]; 还有一个问题，我在保存图片的时候，总是保存不上，打开之后就是一片黑的，没有图。有的说是内存不够了，但是我觉得这个图没有多大。PS：我换了好几种格式还是保存不上。 |
3楼2010-03-29 09:22:44
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
4楼2010-03-29 12:10:18
5楼2010-03-30 09:06:48
