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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » matlab解方程
232/3返回列表上一页123下一页
查看: 1871  |  回复: 22
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

津津爱睡觉

金虫 (小有名气)
引用回帖:
9楼: Originally posted by 月只蓝 at 2015-07-20 14:23:09
1stOpt软件代码：
Variable  x,y ;
Function x=1/(x2*49.005^2)*(ln((y/x1)*((x1-49.005)/(y-49.005)))-49.005*(1/y-1./x1));
Data;
         0   83.9499
    0.2500   69.1806
    0.5000   63.7801
     ...

clear,clc
x=[0,0.25,0.5,0.75,1,1.25,1.5,1.75,2,2.5,3,3.5,4,4.5,5,5.5,6,6.5]';
y=[83.94986152,69.18056238,63.78006779,62.16487246,60.65377495,59.89822619,58.46520206,57.7096533,56.95410455,55.769572,54.95441885,53.89832959,53.14278084,52.14041949,51.64595481,51.34373531,51.11791018,50.85514711]';
plot(x,y,'k+');
hold on
X=[x,y];n=length(X);
fx=@(b,X)(1/(b(2)*49.005^2)).*[log((X(:,2)./b(1)).*((b(1)-49.005)./(X(:,2)-49.005)))-49.005.*(1./X(:,2)-1./b(1))]-X(:,1);
b=[83.94986152,0.000124];
for l=1:1
    b=lsqcurvefit(fx,b,X,y);
   
end
b
fx=@(b,x,y)(1/(b(2)*49.005^2)).*[log((y./b(1)).*((b(1)-49.005)./(y-49.005)))-49.005.*(1./y-1./b(1))]-x
y1=fx(b,x,y);
plot(x,y1,'r');
hold on
c=corrcoef(y,y1)
版主能帮我看这个拟合代码哪里有错误吗？你发的那个链接我有些地方看不懂
赞一下 回复此楼
有梦一起追
11楼2015-07-20 14:39:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

津津爱睡觉

金虫 (小有名气)
引用回帖:
7楼: Originally posted by 月只蓝 at 2015-07-20 12:42:34
先拟合出：
x2   =              0.000124354608284474
x1     =            71.0611809557734

相关系数(R): 0.99453642822085
相关系数之平方(R^2): 0.989102707058286
决定系数(DC): 0.989102707058287

代入之后用MA ...

function fit_nonl
clear all;clc
format long

data=[ 83.94986152    0
       69.18056238    0.25   
       63.78006779    0.5
       62.16487246    0.75
       60.65377495    1
       59.89822619    1.25
       58.46520206    1.5
       57.7096533     1.75
       56.95410455    2
       55.769572      2.5
       54.95441885    3
       53.89832959    3.5
       53.14278084    4
       52.14041949    4.5
       51.64595481    5
       51.34373531    5.5
       51.11791018    6
       50.85514711    6.5
];
xspan=data(:,2);    %x的数据，在此输入
Texp=data(:,1);    %T的数据，在此输入

k0=[83.94986152 0.0001];
lb=[];
ub=[];
OPTIONS=optimset('MaxFunEvals',1000);
[k,resnorm,residual,exitflag,output,lambda,jacobian] = ...
    lsqnonlin(@ObjFunc,k0,lb,ub,OPTIONS,xspan,Texp);

ci = nlparci(k,residual,jacobian);
%residual;
fprintf('\n\n拟合结果:\n')
fprintf('\n\t参数 E1 = %.16f',k(1))
fprintf('\n\t参数 E2 = %.16f',k(2))

y=KineticsEqs(xspan,k);
R2=1-sum((Texp-y).^2)./sum((Texp-mean(y)).^2);
fprintf('\n\tR^2 = %.16f',R2);
figure
plot(xspan,KineticsEqs(xspan,k),'b',xspan,Texp,'or'),legend('计算值','实验值','Location','Best')
function f = ObjFunc(k,xspan,Texp)
f=KineticsEqs(xspan,k)-Texp;

%------------------------------------------------------------------------
function xt = KineticsEqs(x,k)
E1=k(1);E2=k(2);
xt=1/(E2*49.005^2).*[log((x./E1).*((E1-49.005)./(x-49.005))-49.005.*(1./x-1/E1)];
版主版主，照你告诉的我的网址改了下代码，还是有错。。帮我看看哪里错了好不好
赞一下 回复此楼
有梦一起追
12楼2015-07-20 15:13:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

津津爱睡觉

金虫 (小有名气)
引用回帖:
9楼: Originally posted by 月只蓝 at 2015-07-20 14:23:09
1stOpt软件代码：
Variable  x,y ;
Function x=1/(x2*49.005^2)*(ln((y/x1)*((x1-49.005)/(y-49.005)))-49.005*(1/y-1./x1));
Data;
         0   83.9499
    0.2500   69.1806
    0.5000   63.7801
     ...

为啥我照着你那个在1sOpt中输的，得到完全不同的结果啊。。拟合参数都是零点几。。都要改哪些条件吗。。求版主帮忙
赞一下 回复此楼
有梦一起追
13楼2015-07-20 15:53:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
13楼: Originally posted by 津津爱睡觉 at 2015-07-20 15:53:10
为啥我照着你那个在1sOpt中输的，得到完全不同的结果啊。。拟合参数都是零点几。。都要改哪些条件吗。。求版主帮忙...

把拟合结果贴出来。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
14楼2015-07-20 16:13:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

津津爱睡觉

金虫 (小有名气)
引用回帖:
9楼: Originally posted by 月只蓝 at 2015-07-20 14:23:09
1stOpt软件代码：
Variable  x,y ;
Function x=1/(x2*49.005^2)*(ln((y/x1)*((x1-49.005)/(y-49.005)))-49.005*(1/y-1./x1));
Data;
         0   83.9499
    0.2500   69.1806
    0.5000   63.7801
     ...

Variable  x,y ;
Function x=1/(x2*89.005^2)*(ln((y/x1)*((x1-89.005)/(y-89.005)))-89.005*(1/y-1./x1));
Data;
         0   241.2453915
    0.2500   183.3464975
    0.5000   155.1400864
    0.7500    144.0572731
    1.0000   130.1157019
    1.2500   123.8167847
    1.5000   119.6818387
    1.7500   114.7628452
    2.0000   112.2128814
    2.2500   109.1828393
    2.5000   107.8005032
    3.0000   103.5545843
    3.5000   101.0962938
    4.0000   97.44383853
    4.5000   96.25691121
    5.0000   93.11348376
    5.5000   91.75768466
    6.0000   91.33309277
    6.5000   90.457372
月只蓝，帮我看看换个方程怎么做不出来了，是什么样的原因呢？麻烦你啦
赞一下 回复此楼
有梦一起追
15楼2015-07-20 17:29:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
15楼: Originally posted by 津津爱睡觉 at 2015-07-20 17:29:16
Variable  x,y ;
Function x=1/(x2*89.005^2)*(ln((y/x1)*((x1-89.005)/(y-89.005)))-89.005*(1/y-1./x1));
Data;
         0   241.2453915
    0.2500   183.3464975
    0.5000   155.1400864
    0.7500 ...

方程 形式不对 不用点除

[ 发自小木虫客户端 ]
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
16楼2015-07-20 17:45:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

津津爱睡觉

金虫 (小有名气)
引用回帖:
16楼: Originally posted by 月只蓝 at 2015-07-20 17:45:21
方程 形式不对 不用点除
...

把点删了也不对的，只有一条线。。R2什么都没有
赞一下 回复此楼
有梦一起追
17楼2015-07-20 19:19:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

津津爱睡觉

金虫 (小有名气)
引用回帖:
14楼: Originally posted by 月只蓝 at 2015-07-20 16:13:22
把拟合结果贴出来。...

??? Error: File: fit_nonl.m Line: 51 Column: 82
Unbalanced or unexpected parenthesis or bracket.
只有这一句。。因为我方程是隐函数，所以不太会弄。。之前的显函数都可以做出来。。这个怎么也不行
赞一下 回复此楼
有梦一起追
18楼2015-07-20 19:24:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

津津爱睡觉

金虫 (小有名气)
引用回帖:
14楼: Originally posted by 月只蓝 at 2015-07-20 16:13:22
把拟合结果贴出来。...

迭代数: 16
计算用时(时:分:秒:微秒): 00:00:00:431
优化算法: 麦夸特法(Levenberg-Marquardt) + 通用全局优化法
计算结束原因: 达到收敛判断标准
均方差(RMSE): 3.36438388822804
残差平方和(SSE): 215.0625
相关系数(R): 0
相关系数之平方(R^2): 0
决定系数(DC): -1.90186418109188
卡方系数(Chi-Square): 51.75
F统计(F-Statistic): 0

参数        最佳估算
----------        -------------
x2                 0.063293183920905
x1                 0.320298498496413

====== 结果输出 =====

No        实测值x        计算值x
1        0        0.0000000
2        0.25        0.0000000
3        0.5        0.0000000
4        0.75        0.0000000
5        1        0.0000000
6        1.25        0.0000000
7        1.5        0.0000000
8        1.75        0.0000000
9        2        0.0000000
10        2.25        0.0000000
11        2.5        0.0000000
12        3        0.0000000
13        3.5        0.0000000
14        4        0.0000000
15        4.5        0.0000000
16        5        0.0000000
17        5.5        0.0000000
18        6        0.0000000
19        6.5        0.0000000
赞一下 回复此楼
有梦一起追
19楼2015-07-20 19:25:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

津津爱睡觉

金虫 (小有名气)
引用回帖:
14楼: Originally posted by 月只蓝 at 2015-07-20 16:13:22
把拟合结果贴出来。...

之前不知道要改算法。。改完算法可以了。。谢谢你。。麻烦你了一下午。。感谢感谢。。帮了我大忙
赞一下 回复此楼
有梦一起追
20楼2015-07-20 19:32:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 津津爱睡觉 的主题更新
232/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +17 Lanmanbaby 2026-08-09 29/1450 2026-08-12 00:24 by zhanghaozhu
[基金申请] 是这周出结果还是下周出结果？ +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +8 Tide man 2026-08-10 9/450 2026-08-11 20:39 by beefly
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] FileCode能看出啥？ +7 要乐观耀哥 2026-08-10 17/850 2026-08-11 10:50 by 要乐观耀哥
[基金申请] filecode +9 documentary 2026-08-10 9/450 2026-08-11 10:08 by wk7465
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题，想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读？ +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 固定端突然变了，今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 好消息？这个有何含义？？？ +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
信息提示
关闭
请填处理意见
关闭