24小时热门版块排行榜    

查看: 1737  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

fengxingmin

新虫 (初入文坛)

[求助] 用Matlab编程,大家看看我编的程序,谁知道哪里错了啊?

我用matlab编程,来回归方程里面的两个参数,由于是初学者,不知道这个程序哪里错了,大家看看帮我改一下吧,,谢谢了。!
function jscs;
clear all; clc
%实验数据
x1=[0.00000 0.0098 0.0495 0.0811 0.1020 0.1497 0.2494 0.3213 0.4501 0.6905 0.9112 0.9656 0.9818 0.9999];
x2=1-x1;
y1e=[0.00000 0.01480 0.06030 0.09510 0.11020 0.14540 0.19880 0.22710 0.26250 0.32400 0.55510 0.79530 0.89990 0.99990];
y2e=1-y1e;
p=[101.325 101.325 101.325 101.325 101.325 101.325 101.325 101.325 101.325 101.325 101.325 101.325 101.325 101.325];
ps1=[43.2649 42.6417 41.7209 41.4176 41.2667 41.2667 41.5690 42.0260 43.4219 49.7470 71.2944 85.8818 92.1282 101.2401];
ps2=[101.7462 100.1531 97.8025 97.0291 96.6444 96.6444 97.4152 98.5809 102.1477 118.4094 174.7990 213.6398 230.4018 254.975];
a0=[1 1];
[a,resnorm]=lsqnonlin(@objFunc,a0,[],[]);
ci=nlparci(a,residual,jacobian);
fprintf('\nEstimated Parameters:\n')
fprintf('\t%.3f±%.3f\n',ci(1,2)-a(1))
fprintf('\t%.3f±%.3f\n',ci(2,2)-a(2))
fprintf('\tThe sum of the squares is:%3f’,resnorm)
% -------------------
function f=ObjFunc(a,x1,x2,ps1,ps2,p)
y1c=(ps1.*x1.*exp(a(1)*x2./(x1.+a(1)*x2.)))/(p.*(x1.+a(1)*x2.)*exp(a(2)*x2./(x2.+a(2)*x1.)));
y2c=(ps2.*x2.*exp(a(2)*x1./(x2.+a(2)*x1.)))/(p.*(x2.+a(2)*x1.)*exp(a(1)*x1./(x1.+a(1)*x2.)));
f=sqrt((y1e-y1c)^2+(y2e-y2c)^2);
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengxingmin

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by lurencyj at 2012-03-17 22:23:47:
首先,不知道你的程序是什么东东;
第二,程序里面所有的分号要从全角改成半角;
第三,程序里面一些单引号,错写成反引号。
再次,一个忠告,如果看不懂matlab,不知道它的语法,建议楼主慎用,因为你自己都不 ...

我是刚开始学Matlab,由于要回归参数,得自己学着编程序,是不是我编的离正确的程序差好远啊??
4楼2012-03-18 14:39:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

lurencyj

木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
首先,不知道你的程序是什么东东;
第二,程序里面所有的分号要从全角改成半角;
第三,程序里面一些单引号,错写成反引号。
再次,一个忠告,如果看不懂matlab,不知道它的语法,建议楼主慎用,因为你自己都不知道程序在做什么。
很女子很弓虽大
2楼2012-03-17 22:23:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emanlee

木虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
你运行的时候,在Command窗口会有错误提示,指出哪行有错误,是什么错误,这样可以很好调试程序。
3楼2012-03-18 07:51:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengxingmin

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by emanlee at 2012-03-18 07:51:53:
你运行的时候,在Command窗口会有错误提示,指出哪行有错误,是什么错误,这样可以很好调试程序。

恩好的,我试试,谢谢了
5楼2012-03-18 14:50:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +8 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:40 by l0VvVHGBGRLv
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 11:18 by l0VvVHGBGRLv
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-30 10:33 by winsaint
[基金申请] 29号明天会评吗 +3 笨笨唐 2026-08-28 3/150 2026-08-30 09:29 by cww8181
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 08:33 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:17 by ZPa0EcMwuECS
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 06:35 by ZPa0EcMwuECS
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 ASdOkHsho7FD 2026-08-28 10/500 2026-08-30 06:11 by ZPa0EcMwuECS
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 06:10 by ZPa0EcMwuECS
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 05:48 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 13/650 2026-08-30 05:38 by ZPa0EcMwuECS
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 02:40 by ZPa0EcMwuECS
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
信息提示
请填处理意见