| 查看: 719 | 回复: 3 | |||
[求助]
matlab拟合参数 已有1人参与
|
|
已知PA0 T x t dx/dt=k*(PA0*(1-x))^n 根据已知的数拟合k 和n 程序在附件中，就一直出现 ‘too many input arguments’ 求解决 |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有199人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何用Matlab曲线拟合一个多输入参数，而且其中一个输入参数与另一个参数有关呢？
已经有6人回复
- 拜求！matlab拟合动力学参数
已经有14人回复
- 求助用matlab拟合动力学方程
已经有9人回复
- matlab 拟合方法
已经有4人回复
- matlab非线性拟合模型参数，大神你在哪！
已经有73人回复
- matlab 拟合求助
已经有8人回复
- 关于matlab拟合
已经有4人回复
- 用matlab里的nlinfit如何拟合分段函数
已经有8人回复
- 用Matlab做多变量非线性拟合求参数值，求各位大侠鼎力相助！
已经有3人回复
- 用Matlab线性拟合求解两个参数，十万火急！
已经有9人回复
- 怎样用一组参数同时拟合两个曲线--matlab
已经有5人回复
- MATLAB微分方程参数拟合问题，求大神
已经有7人回复
- Matlab 含有两个参数的公式拟合
已经有4人回复
- matlab数据拟合问题求助
已经有7人回复
- Matlab多项式拟合中参数的不确定度及相关系数
已经有3人回复
- matlab 拟合反应动力学参数结果很差。大家帮忙看一下
已经有14人回复
- matlab如何编写共享参数拟合程序
已经有8人回复
- matlab参数拟合
已经有6人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- matlab拟合方程参数时初值的选择
已经有15人回复
- matlab拟合模型参数
已经有5人回复
- matlab非线性参数拟合问题
已经有7人回复
- 求助matlab隐函数参数拟合
已经有16人回复
- 【求助】用matlab最优化方法进行参数拟合
已经有17人回复
- 【求助】使用MATLAB拟合出按照公式的系数
已经有19人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2014-07-06 18:17:41
|
function biaoguan clear all;clc format long; %global T PA0 tspan x x0 T=[613 633 663]; F0=[15 17.5 20 22.5 25]; %t=0.03/0.78*134/0.004./F0; t=[85.8974 73.6264 64.4231 57.2650 51.5385]; x1=[0.632 0.581 0.540 0.496 0.468]; x2=[0.600 0.551 0.515 0.480 0.451]; x3=[0.551 0.492 0.459 0.427 0.402]; x=[x1;x2;x3]; x0=[0.632 0.600 0.551]; PA0=[0.001865374 0.002289233 0.00239448]; beta0=[1.1e5 100000 1]; lb=[0 0 0]; ub=[inf inf inf]; tspan =t; [beta,resnorm,resid,exitflag,output,lambda,jacobian] = ... lsqnonlin(@Opt,beta0,[],[],[],tspan,x0,T,PA0,x) ci = nlparci(beta,resid,jacobian); beta function f = Opt(beta,tspan,x0,T,PA0,x) [t xc1] = ode45(@Kin,tspan,x0(1),[],beta,T(1),PA0(1),x(1,  );[t xc2] = ode45(@Kin,tspan,x0(2),[],beta,T(2),PA0(2),x(2, );[t xc3] = ode45(@Kin,tspan,x0(3),[],beta,T(3),PA0(3),x(3, );f1= xc1' - x(1, ;f2= xc2' - x(2, ;f3= xc3' - x(3, ;f=[f1;f2;f3] % ------------------------------------------------------------------ function dxdt = Kin(t,beta,T,PA0,x) dxdt=beta(1).*exp(-beta(2)./(8.3145*T)).*(PA0.*(1-x)).^beta(3); |
» 本帖已获得的红花（最新10朵）
sara_ecust3楼2014-07-06 20:30:46
4楼2014-07-06 20:31:52
5