由买买提看人间百态

boards

本页内容为未名空间相应帖子的节选和存档,一周内的贴子最多显示50字,超过一周显示500字 访问原贴
Physics版 - 一个关于用matlab解微分方程的小问题 (转载)
相关主题
请教Matlab高手求指教:MATLAB的优化工具可以用于常微分方程组吗?
请教初级问题。谢谢!数值计算含参变量的Fourier积分 (转载)
请教一个空气降温的问题 (转载)请问一个关于解Schrodinger方程的问题
请大家推荐一本经典的统计物理学书 请教KK变换问题
问个问题:知道 偏微分方程 和 matlab 的请帮忙 (转载)国内往弦论上投入那么多钱
如何对隐函数做data fitting?对Grassmann viariable的积分?
一个/组哈密顿运动方程的求解看Polchinski的书需要什么基础啊
有人知道浮点数的压缩算法吗?有损也没关系 (转载)咆哮体:不难证明,学物理的人你伤不起啊!!!!!(转载)
相关话题的讨论汇总
话题: xdata话题: y0话题: tspan话题: ode23s话题: numcal
进入Physics版参与讨论
1 (共1页)
k******o
发帖数: 61
1
【 以下文字转载自 Mathematics 讨论区 】
发信人: kafeimao (咖啡毛), 信区: Mathematics
标 题: 一个关于用matlab解微分方程的小问题
发信站: BBS 未名空间站 (Sun May 8 17:29:15 2011, 美东)
大家好,我要用matlab做变系数微分方程组的参数拟合,参考其他论坛上的教程,我修
改了自己的程序如下:
function dy=dydt(t,y,k)
dy=zeros(2,1)
dy(1)=0.0321*k(1)*(k(2)-y(1))-k(3)*y(1)-y(2)
dy(2)=0.25*k(4)*exp(-k(4)*t)*k(2);
function y=numcal(k,x)
global y0
tspan=[0 max(x)];
[m,n]=size(x);
[tt yy] = ode23s(@dydt,tspan,y0,[],k);
yc=spline(tt',yy',x);
y=yc;
主程序
clc;clear;
global y0
xdata=[1,2,3,4];ydata=data[3,4,5,6];
k0=[0,0,0,0,];%要识别参数的初始值
lb=[0,0,0,0];%要识别参数的下限
ub=[6.5,5,6];%要识别参数的上限;课根据参数的范围自己设定。
options=optimset('TolFun',1e-20,'TolX',1e-20,'MaxFunEvals',100,'Algorithm','
trust-region-reflective','Display','iter');
beta=lsqcurvefit(@numcal,k0,xdata,ydata,lb,ub,options);
Jc=numcal(k,xdata);
plot(xdata,ydata,'o',xdata,Jc);%画图查看参数识别后目标值和实验目标值之间的差
异;
结果还是出现错误:“??? Error using ==> odearguments at 19
When the first argument to ode23s is a function handle, the tspan and y0
arguments must be supplied.
Error in ==> ode23s at 188
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...
Error in ==> numcal at 8
[tt yy] = ode23s(@dydt,tspan,y0,[],k);
Error in ==> lsqcurvefit at 209
initVals.F = feval(funfcn_x_xdata{3},xCurrent,XDATA,varargin{:});
Caused by:
Failure in initial user-supplied objective function evaluation.
LSQCURVEFIT cannot continue.”
请问我的错误在哪里,纠结了几天了,还是查不出来,谢谢大家。
1 (共1页)
进入Physics版参与讨论
相关主题
咆哮体:不难证明,学物理的人你伤不起啊!!!!!(转载)问个问题:知道 偏微分方程 和 matlab 的请帮忙 (转载)
電磁場能動量張量的計算問題,廣義相對論下的如何对隐函数做data fitting?
求助!!!陳省身《微分幾何講義》題目!!!一个/组哈密顿运动方程的求解
物理系的学姐不好追(zz) (转载)有人知道浮点数的压缩算法吗?有损也没关系 (转载)
请教Matlab高手求指教:MATLAB的优化工具可以用于常微分方程组吗?
请教初级问题。谢谢!数值计算含参变量的Fourier积分 (转载)
请教一个空气降温的问题 (转载)请问一个关于解Schrodinger方程的问题
请大家推荐一本经典的统计物理学书 请教KK变换问题
相关话题的讨论汇总
话题: xdata话题: y0话题: tspan话题: ode23s话题: numcal