24小时热门版块排行榜

返回列表

当前主题已经存档。

飞翔的海鸥

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 463.3
红花: 1
帖子: 80
在线: 5.9小时
虫号: 406142
注册: 2007-06-18
性别: GG
专业: 激光

[交流] 【求助】matlab程序问题急

function pdex4
clc
clear
m = 0;
x = linspace(-50,150,30);
t = linspace(0,250,30);
u = inline('u(x,t)');
sol = pdepe(m,@pdex4pde,@pdex4ic,@pdex4bc,x,t);
u1 = sol(:,:,1);
u2 = sol(:,:,2);
u3 = sol(:,:,3);
u4 = sol(:,:,4);
u5 = sol(:,:,5);
u6 = sol(:,:,6);
u7 = sol(:,:,7);

figure(1)
mesh(x,t,(abs(u7)).^2)
title('(abs(a)).^2')
xlabel('Distance x')
ylabel('Time t')
% A solution profile can also be illuminating.

figure(2)
surf(x,t,abs(u2),'FaceColor',[1 1 1])
title('u2(x,t)')
xlabel('Distance x')
ylabel('Time t')
% A solution profile can also be illuminating.
figure(3)
plot(t,abs(u7(end,

))
title('Solution at x =0 ')
xlabel('time')
ylabel('u7(0,

')

% --------------------------------------------------------------
function [c,f,s] = pdex4pde(x,t,u,DuDx,DuDt)
c = [1; 1; 1; 1; 1; 1; 1];
f = [0; 0; 0; 0; 0; 0; u];
ganmp=1.0;
gam=0;
g=1.0;
gam1=2.0;
gam2=2.0;
cottht=10.*(1-0.5.*tanh(0.1.*(t-65))+0.5.*tanh(0.1.*(t-185)));
y1=-(ganmp+i.*gam1).*u(1)-(i./2).*g.*(u(7).*(u(5)-u(4))-cottht.*u(2));
y2=i.*(gam2-gam1).*u(2)-(i./2).*g.*(u(7).*u(3)-cottht.*u(1));
y3=-(ganmp-i.*gam2).*u(3)-(i./2).*g.*(u(2).*conj(u(7))+cottht.*(u(6)-u(5)));
y4=u(5)-(i./2).*g.*(u(7).*conj(u(1))-u(1).*conj(u(7)));
y5=-(2+gam).*u(5)-(i./2).*g.*(u(1).*conj(u(7))-u(7).*conj(u(1))+cottht.*(conj(u(3))-u(3)));
y6=u(5)-(i./2).*g.*cottht.*(u(3)-conj(u(3)));
y7=(i./2).*g.*u(1);
s = [y1; y2; y3; y4; y5; y6; y7];
% --------------------------------------------------------------
function u0 = pdex4ic(x);
u0 = [0;0;0;1;0;0;0.01.*exp(-((x./15).^2))];
% --------------------------------------------------------------
function [pl,ql,pr,qr] = pdex4bc(xl,ul,xr,ur,t)
pl = [ul(1); ul(2); ul(3); ul(4); ul(5); ul(6); 0];
ql = [0; 0; 0; 0; 0; 0; 1];
pr = [ur(1); ur(2); ur(3); ur(4); ur(5); ur(6); 0];
qr = [0; 0; 0; 0; 0; 0; 1];

方程7 是 DuDt+DuDx=i=(i./2).*g.*u(1) ,是随空间变化的，但结果怎么不随空间变化了，谁能帮我改改这个程序。
感谢所有看到此程序的人，特别感谢帮我修改的人！！！

[ Last edited by lvzhu2007 on 2008-4-19 at 15:25 ]

回复此楼