24小时热门版块排行榜    

查看: 407  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 情弦风 的主题更新
信息提示
请填处理意见