24小时热门版块排行榜    

查看: 691  |  回复: 4

卡萨丁

新虫 (小有名气)

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

我想拟合一个 三自变量 11参数的方程  已经用matlab试过了但结果不太好
有matlab程序,但不知道1stopt里怎么写
求代跑,价钱好商量
回复此楼
人的一切感情都是脑内的化学物质作用
已阅   回复此楼   关注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

» 本帖已获得的红花(最新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的回帖
相关版块跳转 我要订阅楼主 卡萨丁 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 【考研调剂】化学专业 281分,一志愿四川大学,诚心求调剂 +12 吃吃吃才有意义 2026-03-19 13/650 2026-03-24 10:06 by NCUT_Yu
[考研] 336求调剂 +4 收到VS 2026-03-20 4/200 2026-03-23 19:02 by macy2011
[考研] 328求调剂,英语六级551,有科研经历 +7 生物工程调剂 2026-03-17 12/600 2026-03-23 18:18 by YMU施老师
[考研] 工科0856求调剂 +5 沐析汀汀 2026-03-21 5/250 2026-03-23 17:56 by 海瑟薇-
[考研] 307求调剂 +3 余意卿 2026-03-21 3/150 2026-03-23 10:32 by Iveryant
[考研] 291求调剂 +5 孅華 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[考研] 317求调剂 +12 申子申申 2026-03-19 18/900 2026-03-22 22:23 by luoyongfeng
[考研] 293求调剂 +3 涛涛Wjt 2026-03-22 5/250 2026-03-22 22:21 by jiangpengfei
[考研] 石河子大学(211、双一流)硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 306求调剂 +5 来好运来来来 2026-03-22 5/250 2026-03-22 16:17 by BruceLiu320
[考研] 生物学调剂 +5 Surekei 2026-03-21 5/250 2026-03-22 14:39 by tcx007
[考研] 318求调剂 +4 plum李子 2026-03-21 7/350 2026-03-22 14:17 by ColorlessPI
[考研] 化学调剂 +5 yzysaa 2026-03-21 5/250 2026-03-21 22:12 by peike
[考研] 一志愿深大,0703化学,总分302,求调剂 +4 七月-七七 2026-03-21 4/200 2026-03-21 18:20 by 学员8dgXkO
[考研] 302求调剂 +12 呼呼呼。。。。 2026-03-17 12/600 2026-03-21 17:29 by ColorlessPI
[考研] 0805材料320求调剂 +3 深海物语 2026-03-20 3/150 2026-03-21 15:46 by 无际的草原
[考研] 材料与化工(0856)304求 B区 调剂 +3 邱gl 2026-03-21 3/150 2026-03-21 13:47 by lature00
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
[考研] 求调剂 +3 eation27 2026-03-20 3/150 2026-03-20 19:32 by JourneyLucky
信息提示
请填处理意见