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

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 有人能看一下我宿舍管理系统的报修功能该怎么改啊?链表那里总是越界
  • ¥15 cs loadimage运行不了,easyx也下了,没有用
  • ¥15 r包runway详细安装教程
  • ¥15 Html中读取Json文件中数据并制作表格
  • ¥15 谁有RH342练习环境
  • ¥15 STM32F407 DMA中断问题
  • ¥15 uniapp连接阿里云无法发布消息和订阅
  • ¥25 麦当劳点餐系统代码纠错
  • ¥15 轮班监督委员会问题。
  • ¥20 关于变压器的具体案例分析