马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
function [yout,xcout,tcout]=xdkeuler(f,x0,t0,td,h,c0)<BR>m=(td-t0)/h;<BR>xcout{1}=x0;<BR>tcout(1)=t0;<BR>for i=1:m;<BR> yout{i}=c0*xcout{i}<BR> tcout(i+1)=tcout(i)+h;<BR> temp=subs(f,{'t','y'},{tcout(i),xcout{i}});<BR> xcout{i+1}=xcout{i}+h*temp;<BR>end<BR>我这里f函数是f=a*x+u a是方阵,u是向量,变量x也是向量,a和u的值在workspace里有当我把函数‘a*x+u’输进去后表达式会不会自动去把a和u附值附进去?? |