afdsfr 2022-08-20 18:10 采纳率: 33.3%
浏览 718
已结题

Matlab报错“function value and YDATA sizes are not equal”

在拟合过程中,Matlab报错“function value and YDATA sizes are not equal”,报错图如下

img

望能人解答,谢谢

代码附上
#main
clear all;clc;
%读取数据
data=xlsread('工作簿2.xls');
x0=data(:,1);
y0=data(:,2);
%拟合
beta0=[0.005,50,2];
% lb=[0,0,0];
% ub=[0.2,1000,5];
[p]=lsqcurvefit(@func,beta0,x0,y0);
%结果展示
xx=[1:10:360];
yy=func(p,xx);
figure;
plot(x0,y0,'ko');
hold on;
plot(xx,yy);
grid on;

#func

function yhat = func( beta0,X )
p00=beta0(1)
p10=beta0(2)
p20=beta0(3)

f=0.4-p00X/(1+X/p10);
yhat=p20
sqrt((52079166+10835766ff-62914932*f)1000000/(10417-3817f)/2020);

end

  • 写回答

4条回答 默认 最新

  • 斗迷飞鸟 2022-08-20 21:08
    关注

    (1)原函数代码存在的问题

    function yhat = original_func( beta0,X )
    p00=beta0(1)  //应该打分号减少不必要的屏幕输出
    p10=beta0(2)
    p20=beta0(3)
    
    f=0.4-p00*X/(1+X/p10);  //这里X是向量不是标量,因此"X/(1+X/p10)"应该修改为"X./(1+X/p10)"
    // (f是向量不是标量,"10835766/f"应修改为"10835766./f""1000000/(10417-3817*f)"应修改为"1000000./(10417-3817*f)")
    yhat=p20*sqrt((52079166+10835766/f-62914932*f)*1000000/(10417-3817*f)/2020); 
    end
    

    (2)修改后的函数代码

    function yhat = func( beta0,X )
    p00=beta0(1);
    p10=beta0(2);
    p20=beta0(3);
    
    f=0.4-p00*X./(1+X/p10);
    yhat=p20*sqrt((52079166+10835766./f-62914932*f)*1000000./(10417-3817*f)/2020);
    
    end
    

    (3)代码验证
    由于没有“工作簿2.xls”,修改主程序,通过模型数据加随机数得到待拟合数据,修改后的主程序如下:

    clear all;clc;
    %读取数据
    %data=xlsread('工作簿2.xls');
    % x0=data(:,1);
    % y0=data(:,2);
    x0=[1:10:360];
    y0=func([0.005,50,2],x0)+100*rand(size(x0));
    %拟合
    beta0=[0.005,50,2];
    % lb=[0,0,0];
    % ub=[0.2,1000,5];
    [p]=lsqcurvefit(@func,beta0,x0,y0);
    %结果展示
    xx=[1:10:360];
    yy=func(p,xx);
    figure;
    plot(x0,y0,'ko');
    hold on;
    plot(xx,yy);
    grid on;
    axis tight;
    

    程序运行结果如下:

    img

    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论 编辑记录
查看更多回答(3条)

报告相同问题?

问题事件

  • 系统已结题 9月2日
  • 已采纳回答 8月25日
  • 创建了问题 8月20日

悬赏问题

  • ¥30 自适应 LMS 算法实现 FIR 最佳维纳滤波器matlab方案
  • ¥15 lingo18勾选global solver求解使用的算法
  • ¥15 全部备份安卓app数据包括密码,可以复制到另一手机上运行
  • ¥15 Python3.5 相关代码写作
  • ¥20 测距传感器数据手册i2c
  • ¥15 RPA正常跑,cmd输入cookies跑不出来
  • ¥15 求帮我调试一下freefem代码
  • ¥15 matlab代码解决,怎么运行
  • ¥15 R语言Rstudio突然无法启动
  • ¥15 关于#matlab#的问题:提取2个图像的变量作为另外一个图像像元的移动量,计算新的位置创建新的图像并提取第二个图像的变量到新的图像