89456418961 2022-02-28 16:33 采纳率: 84.6%
浏览 674
已结题

MATLAB fsolve:串联的矩阵的维度不一致,如何解决?

a=9;
d=0.6;
xr1=1;
xg1=1;
xb1=1.2;

x0=[2,2,2,2,2,2,0.6];      % 初始参考点

fun=@(x)[

(x(1)+x(2)+x(1)+2*x(7))-a ,...

(x(2)*3^0.5+d+(x(2)*3^0.5+d)/2)-a   ,...
   
(3^0.5*(x(1)*x(2)+x(2)*x(2)/2)/a/a)/(3^0.5*(x(3)*x(4)+x(4)*x(4)/2)/a/a)-xr1/xg1 ,...
(3^0.5*(x(3)*x(4)+x(4)*x(4)/2)/a/a)/(3^0.5*(x(5)*x(6)+x(6)*x(6)/2)/a/a)-xg1/xb1 ,...
(3^0.5*(x(1)*x(2)+x(2)*x(2)/2)/a/a)/(3^0.5*(x(5)*x(6)+x(6)*x(6)/2)/a/a)-xr1/xb1 ,...
    ];
x=fsolve(fun,x0);   %调用函数

Lr1=x(1)
Lr2=x(2)
Lg1=x(3)
Lg2=x(4)
Lb1=x(5)
Lb2=x(6)
Ld =x(7)

错误使用 vertcat
串联的矩阵的维度不一致。

出错 fsolve (line 230)
fuser = feval(funfcn{3},x,varargin{:});

出错 SIX (line 41)
x=fsolve(fun,x0); %调用函数

原因:
Failure in initial objective function evaluation.
FSOLVE cannot continue.

  • 写回答

1条回答 默认 最新

  • joel_1993 2022-02-28 19:11
    关注

    你好,把函数句柄里面的逗号改成分号就行

    a=9;
    d=0.6;
    xr1=1;
    xg1=1;
    xb1=1.2;
     
    x0=[2,2,2,2,2,2,0.6];      % 初始参考点
     
    fun=@(x)[
    (x(1)+x(2)+x(1)+2*x(7))-a ;...
    (x(2)*3^0.5+d+(x(2)*3^0.5+d)/2)-a   ;...
    (3^0.5*(x(1)*x(2)+x(2)*x(2)/2)/a/a)/(3^0.5*(x(3)*x(4)+x(4)*x(4)/2)/a/a)-xr1/xg1 ;...
    (3^0.5*(x(3)*x(4)+x(4)*x(4)/2)/a/a)/(3^0.5*(x(5)*x(6)+x(6)*x(6)/2)/a/a)-xg1/xb1 ;...
    (3^0.5*(x(1)*x(2)+x(2)*x(2)/2)/a/a)/(3^0.5*(x(5)*x(6)+x(6)*x(6)/2)/a/a)-xr1/xb1 ;...
        ]; % 把函数句柄里面的逗号改成分号就行
    x=fsolve(fun,x0);   %调用函数
     
    Lr1=x(1)
    Lr2=x(2)
    Lg1=x(3)
    Lg2=x(4)
    Lb1=x(5)
    Lb2=x(6)
    Ld =x(7)
    
    

    结果

    
    Lr1 =
    
       1.722973255909702
    
    
    Lr2 =
    
       3.117691453619087
    
    
    Lg1 =
    
       2.367819755812734
    
    
    Lg2 =
    
       2.738056112727723
    
    
    Lb1 =
    
       2.522158681936574
    
    
    Lb2 =
    
       3.038179617322688
    
    
    Ld =
    
       1.218181017277643
    

    有帮助望采纳哟,谢谢啦

    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论

    报告相同问题?

    问题事件

    • 系统已结题 3月9日
    • 已采纳回答 3月1日
    • 创建了问题 2月28日

    悬赏问题

    • ¥15 关于网上一个easyx制作的见缝插针小游戏(c++)
    • ¥15 开地址法双散列函数处理碰撞
    • ¥15 想问一下这个是什么情况 虚拟机Linux打不开了
    • ¥15 联通光猫掉注册了怎么重新注册上去
    • ¥15 关于unity开发steamvr程序遇到的问题
    • ¥60 求tc downloader的下载方式
    • ¥15 华为 快捷方式 手电筒 接口
    • ¥15 Qt6.5支不支持Android13开发啊
    • ¥20 网络只能跑一半,应该如何设置
    • ¥20 Python调用百度开发者平台人脸识别接口