3/26/2008

matlab中解析解和数值解一起求

方法就是先用syms 指令将其符号化,然后解方程得到解析解
然后将参数赋值,用eval 指令将其数值解求出。
例子

%hw6prob3

clear;
syms gamma omega delta;
diary hw6prob3.txt;
A1=[-gamma, i*omega,-i*omega,0;i*omega,-(gamma/2+i*delta),0,-i*omega;-i*omega,0,-(gamma/2-i*delta),i*omega;1,0,0,1],

B=[0;0;0;1],

X1=A1\B;
fprintf('The analytical solution for this linear equation is: \n ');
X1,

gamma=1;
omega=2;
delta=0;

fprintf('The numerical solution for this linear equation is: \n ');
X2=eval(X1),

diary off;

没有评论: