方法就是先用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;
没有评论:
发表评论