24小时热门版块排行榜

返回列表

当前主题已经存档。

岸芷汀兰-

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 163
帖子: 200
在线: 8.5小时
虫号: 905122
注册: 2009-11-17
专业: 力学

[交流] 【求助】画三维图【已解决】

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 ]

回复此楼

» 猜你喜欢

面上项目申报已经有3人回复
有时候真觉得大城市人没有县城人甚至个体户幸福已经有5人回复
酰胺脱乙酰基已经有9人回复
CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优已经有5人回复
有70后还继续奋斗在职场上的吗？已经有6人回复
博士延得我，科研能力直往上蹿已经有7人回复
退学或坚持读已经有27人回复
面上基金申报没有其他的参与者成吗已经有5人回复
遇见不省心的家人很难过已经有22人回复

1楼 2010-03-25 16:52:33

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师

应助: 44 (小学生)
金币: 17724.5
红花: 17
帖子: 2413
在线: 546.7小时
虫号: 496517
注册: 2008-01-19
专业: 化学反应工程

★ ★ ★
小木虫(金币+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人)

回复此楼

2楼2010-03-26 16:58:06

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

岸芷汀兰-

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 163
帖子: 200
在线: 8.5小时
虫号: 905122
注册: 2009-11-17
专业: 力学

引用回帖:

Originally posted by change0618 at 2010-03-26 16:58:06:
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) ...

谢谢，我后来也编出来了，用的是个循环。但是图不如你的好看。。。
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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖