应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 资源共享 » 求助 1stopt代跑
51/1返回列表
查看: 654  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

卡萨丁

新虫 (小有名气)

[交流] 求助 1stopt代跑 已有1人参与

我想拟合一个 三自变量 11参数的方程  已经用matlab试过了但结果不太好
有matlab程序，但不知道1stopt里怎么写
求代跑，价钱好商量
回复此楼

» 猜你喜欢
人的一切感情都是脑内的化学物质作用
1楼 2015-12-07 16:23:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡萨丁

新虫 (小有名气)
clc
clear all
I=[
-4        5        25        4986
-4        10        25        5332
-4        16        25        6549
-4        24        25        7495
-4        30        25        8774
-3.5        5        25        5001
-3.5        10        25        5162
-3.5        16        25        6309
-3.5        24        25        7119
-3.5        30        25        8539

-3        10        25        1243
-3        16        25        823
-3        24        25        920
-3        30        25        1195
-2.5        5        25        206
-2.5        10        25        237
-2.5        16        25        328
-2.5        24        25        465
-2.5        30        25        656
-2        5        25        141
-2        10        25        204
-2        16        25        310
-4        5        85        3244
-4        10        85        3452
-4        16        85        4147
-4        24        85        4772
-4        30        85        5484
-3.5        5        85        3224
-3.5        10        85        3330
-3.5        16        85        4077
-3.5        24        85        4699
-3.5        30        85        5384

-3        10        85        1796
-3        16        85        1048
-3        24        85        1122
-3        30        85        1497
-2.5        5        85        271
-2.5        10        85        274
-2.5        16        85        373
-2.5        24        85        522
-2.5        30        85        704
-2        5        85        163
-2        10        85        223
-2        16        85        339
-4        5        125        2713
-4        10        125        2810
-4        16        125        3278
-4        24        125        3812
-4        30        125        4259
-3.5        5        125        2705
-3.5        10        125        2766
-3.5        16        125        3219
-3.5        24        125        3793
-3.5        30        125        4242

-3        10        125        1891
-3        16        125        1713
-3        24        125        1865
-3        30        125        2207
-2.5        5        125        297
-2.5        10        125        324
-2.5        16        125        429
-2.5        24        125        585
-2.5        30        125        743
-2        5        125        169
-2        10        125        240
-2        16        125        354
];
%取自变量　（Vd Vg T） 和因变量
V=I(:,1:1:3);
Y=I(:,4);
%c0=rand(1,11)
%c0=[100 450 0.02 0.5 30 50 -3.01 33 6 88 0.5]                 %也就是k=0 初值
c0=[ 3.0895    0.1222    0.0108    0.6112   15.4216   99.0000    3.0595   33.0000    6.0000   88.0000    0.5000]
Cgs=inline('220+(4500.*(0.0108.*xd(:,2)+0.6112)-15.4216.*xd(:,3)).*(1-tanh(3.0895.*(xd(:,1)+3.0595+0.1222.*tanh(xd(:,2)-6).*tanh(100.*(xd(:,1)+3.0595)))))-xd(:,3).*0.5','x','xd');      %Cgs是函数名。。。
%Cgs=inline('x(1)+(x(2).*(x(3).*xd(:,2)+x(4))-x(5).*xd(:,3)).*(1-tanh(x(6).*(xd(:,1)-x(7)+x(8).*tanh(xd(:,2)-x(9)).*tanh(x(10).*(xd(:,1)-x(7))))))-xd(:,3).*x(11)','x','xd');  
[c,resonorm]=lsqcurvefit(Cgs,c0,V,Y)
%画图
i=[1:1:66]';
r=Cgs(c,V);
plot(i,Y,'r')
hold on
plot(i,r)
赞一下 回复此楼
人的一切感情都是脑内的化学物质作用
2楼2015-12-07 16:28:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
CODE:
Variable x(3),y;
Function  y=p1+(p2*(p3*x2+p4)-p5*x3)*(1-tanh(p6*(x1-p7+p8*tanh(x2-p9)*tanh(p10*(x1-p7)))))-x3*p11;
Data;
-4        5        25        4986
-4        10        25        5332
-4        16        25        6549
-4        24        25        7495
-4        30        25        8774
-3.5        5        25        5001
-3.5        10        25        5162
。。。。。

均方差(RMSE):373.352928646608
残差平方和(SSE):9199899.01571395
相关系数(R): 0.986884864274678
相关系数之平方(R^2): 0.973941735334449
确定系数(DC): 0.973941735334448
卡方系数(Chi-Square): 3438.75695684732
F统计(F-Statistic): 211.06547464271

参数        最佳估算
----------        -------------
p1        -0.267100989415004
p2        6.41932445252313
p3        7.62308716702917
p4        439.158225823665
p5        18.6933809575902
p6        4.51190114503422
p7        -3.1393790750151
p8        0.106414567514015
p9        13.1675584815892
p10        2.5650382236529
p11        -5.17475748709202
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

卡萨丁
3楼2015-12-07 17:19:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡萨丁

新虫 (小有名气)
送红花一朵
引用回帖:
3楼: Originally posted by dingd at 2015-12-07 17:19:16
Variable x(3),y;
Function  y=p1+(p2*(p3*x2+p4)-p5*x3)*(1-tanh(p6*(x1-p7+p8*tanh(x2-p9)*tanh(p10*(x1-p7)))))-x3*p11;
Data;
-4        5        25        4986
-4        10        25        5332
- ...

多谢  
只是结果和matlab跑的差不多， 后面较小的值基本是平的，请问怎么解决
赞一下 回复此楼
人的一切感情都是脑内的化学物质作用
4楼2015-12-07 19:23:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡萨丁

新虫 (小有名气)
求高精度的 来代跑一次
赞一下 回复此楼
人的一切感情都是脑内的化学物质作用
5楼2015-12-09 09:47:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 卡萨丁 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭