UC知道怎样编写一个matlab的程序解决一个微分方程问题（含有延时反馈信号）
来源：百度知道 编辑：UC知道 时间：2024/06/20 11:11:49
dA1/dt=-(1-5i)A1+A1*A2 (其中A1*为A1的共轭项)
dA2/dt=10.8-(1+i)A2-A1A1
其中：A1=x(1)+ix(2);
A2=x(3)+ix(4); %A1，A2是复数,x(1) x(2) x(3) x(4) 为原偏微分方程的解,其随时间t变化%
设： E(t)=/A1(t)/.^2 %E为A1模的平方，也随时间t变化%
既： /A1/.^2=x(1).^2+x(2).^2
下面是我所整理的M文件：（含E处不知道写法对否）%就是对上方方程组求解表达式%
function yx=qwhe(t,x)
yx=zeros(4,1);
yx(1)=-x(1)-5.*x(2)+x(1).*x(3)+x(2).*x(4);
yx(2)=-x(2)+5.*x(1)+x(1).*x(4)-x(2).*x(3);
yx(3)=1.08.*(E(t)-E(t-0.02))-x(3)+x(4)-x(1).*x(1)+x(2).*x(2);
yx(4)=-x(4)-x(3)-2.*x(1).*x(2);
%%%上面含E处是延时信号的差，就是此处不知道怎么表示%%%
执行程序是：
tspan=0:0.01:500;
x0=[1 1 1 1];
[t,x]=ode45('qwhe',tspan,x0);
plot(x(3000:48000,2),x(3000:48000,1))
我想问下怎样修改程序或重新编写使其输出以方程解x(1) ，x(2)为横纵坐标的图像？
希望大家能帮帮我，由于刚接触MATLAB很多基础都没有所以很难修改程序，如果有什么地方叙述不全或者不对 希望大家能知道我的想法，谢谢大家了！！：）
x随时间变化的方程就是原始偏微分方程为：
dA1/dt=-(1-5i)A1+A1*A2 (其中A1*为A1的共轭项)
dA2
x随时间t变化的方程呢?
表关注
PDE algorithms.
adaptmesh - Adaptive mesh generation and PDE solution.
assema - Assemble area integral contributions.
assemb - Assemble boundary condition contributions.
assempde - Assemble a PDE problem.
hyperbolic - Solve hyperbolic problem.
parabolic - Solve parabolic problem.
pdeeig - Solve eigenvalue PDE problem.
pdenonlin - Solve nonlinear PDE problem.
poisolv - Fast solution of Poisson's equation on a rectangular grid.
User interface algorithms and utilities.
pdecirc - Draw circle.
pdeellip - Draw ellipse.
pdemdlcv - Convert MATLAB 4.2c Model M-files for use with MATLAB 5.
pdepoly - Draw polygon.
pderect - Draw rectangle.
pdetool - PDE Toolbox graphical user interface (GUI).
Geometry algorithms.
csgchk - Check validity