24小时热门版块排行榜    

查看: 519  |  回复: 1

情弦风

新虫 (小有名气)

[求助] 非线性方程组求解

function F=myfun(X)
I0=dicomread(\'C:\\Users\\Desktop\\0.dcm\');
I500=dicomread(\'C:\\Users\\Desktop\\500.dcm\');
I800=dicomread(\'C:\\Users\\Desktop\\800.dcm\');
I1000=dicomread(\'C:\\Users\\Desktop\\1000.dcm\');
I0=double(I0);
I500=double(I500);
I800=double(I800);
I1000=double(I1000);
I500=I500./I0;
I800=I800./I0;
I1000=I1000./I0;
A=zeros(208,256);
B=zeros(208,256);
C=zeros(208,256);

for(i=1:208)
    for(j=2:256)
      F(1)=(1-X(1))*exp(-500*X(2))+X(1)*exp(-500*X(3))-I500(i,j);
      F(2)=(1-X(1))*exp(-800*X(2))+X(1)*exp(-800*X(3))-I800(i,j);
      F(3)=(1-X(1))*exp(-1000*X(2))+X(1)*exp(-1000*X(3))-I1000(i,j);
      A(i,j)=X(1);
      B(i,j)=X(2);
      C(i,j)=X(3);
      j=j+1;
    end
    i=i+1;
end

%  X=fsolve(\'myfun\',[0,0,0],optimset(\'Display\',\'off\'))  写在Command窗口的
实际上I1000,I500,I800为208*256 的矩阵
为什么我在解的时候只有一组解????
回复此楼

» 猜你喜欢

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

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

情弦风

新虫 (小有名气)

自己顶!d=====( ̄▽ ̄*)b  有木有解决一下
2楼2015-04-07 17:05:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 情弦风 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 DI1qFHFK0KbE 2026-10-24 5/250 2026-10-26 01:39 by cLOojpOxkX0G
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 DI1qFHFK0KbE 2026-10-24 3/150 2026-10-26 01:37 by cLOojpOxkX0G
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 DI1qFHFK0KbE 2026-10-24 4/200 2026-10-26 01:24 by cLOojpOxkX0G
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 7xHjyB3bAiwg 2026-10-24 5/250 2026-10-26 01:20 by cLOojpOxkX0G
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 7xHjyB3bAiwg 2026-10-23 12/600 2026-10-26 00:54 by cLOojpOxkX0G
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 7xHjyB3bAiwg 2026-10-23 5/250 2026-10-26 00:41 by cLOojpOxkX0G
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 GVGRDgNfv7P4 2026-10-23 10/500 2026-10-26 00:15 by cLOojpOxkX0G
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 GVGRDgNfv7P4 2026-10-23 5/250 2026-10-26 00:08 by cLOojpOxkX0G
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 WqtG0Aem8Ufl 2026-10-24 3/150 2026-10-25 20:24 by cLOojpOxkX0G
[硕博家园] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 WqtG0Aem8Ufl 2026-10-24 5/250 2026-10-25 20:07 by cLOojpOxkX0G
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 3xN9SualEUPt 2026-10-24 3/150 2026-10-25 14:23 by XiYP04It3Evb
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-25 14:20 by XiYP04It3Evb
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-25 13:54 by XiYP04It3Evb
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 DI1qFHFK0KbE 2026-10-24 4/200 2026-10-25 13:18 by XiYP04It3Evb
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 7xHjyB3bAiwg 2026-10-23 8/400 2026-10-25 12:59 by XiYP04It3Evb
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 7xHjyB3bAiwg 2026-10-23 10/500 2026-10-25 12:51 by XiYP04It3Evb
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-24 4/200 2026-10-25 01:22 by XcOYG53d3x4p
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 7xHjyB3bAiwg 2026-10-23 8/400 2026-10-25 00:15 by XcOYG53d3x4p
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 7xHjyB3bAiwg 2026-10-24 3/150 2026-10-24 15:41 by Equinoxhua
[考博] 郑州大学申博 +5 stars~ 2026-10-19 8/400 2026-10-22 15:06 by stars~
信息提示
请填处理意见