求助非线性拟合

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

data aa;
do t=216 to 360 by 24;
input AD@@;
output; end;
cards;
1.15 1.30 1.49 1.78 1.66 1.68 1.71
1.20 1.43 1.59 2.05 2.01 2.02 1.99
1.23 1.29 1.40 1.52 1.67 1.82 1.65
1.28 1.37 1.72 2.04 2.22 2.04 2.11
1.29 1.46 1.81 1.86 2.15 2.10 2.09
1.31 1.56 1.89 1.95 2.12 1.99 2.04
1.35 1.35 1.70 1.79 1.87 1.96 1.87
1.39 1.53 1.70 2.09 2.22 2.13 2.12
1.40 1.53 1.76 1.95 1.98 1.97 1.88
1.47 1.57 1.58 1.96 2.01 1.95 2.02

;

proc nlin best=100;
parameters kin=0,c=0,k=0,R=0;
kout=(R-R*exp(216-t/k)-3*exp(216*k-t*k-1))/k;

model AD=exp((kin-kout)*t+c)+1;

run;
结果是不能运行，请问高手有方法解决吗？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：非线性 parameters Parameter paramete output 动力学

proc nlin ;
parameters kin=0 c=0 k=0 R=0;/*这里不需要逗号，但是初值取的也有问题*/
kout=(R-R*exp(216-t/k)-3*exp(216*k-t*k-1))/k;
model AD=exp((kin-kout)*t+c)+1;
run;

经改变可以运行了，只是初值的改变对迭代的结果影响还是很大的，请问初值的选择有探索的方法吗？谢谢ziyenano ！

sleepsoul 发表于 2013-1-13 00:04
经改变可以运行了，只是初值的改变对迭代的结果影响还是很大的，请问初值的选择有探索的方法吗？谢谢ziyena ...

一元非线性看看散点图吧，看你要要拟合的函数的参数大概在什么位置，
也可以设置初值的形式为 parameters a=0 to 1 初值会在这个范围内搜索
不过设置的范围不要太大

谢谢指导！