永发信息网

matlab y-(a(2)*(1.1*sqrt(x)+a(1))./(0.5*x+52)) Matrix dimensions must agree????

答案:1  悬赏:30  手机版
解决时间 2021-11-30 19:00
matlab y-(a(2)*(1.1*sqrt(x)+a(1))./(0.5*x+52)) Matrix dimensions must agree????
最佳答案
你太可爱了,请你把参数的位置放对好不好,[a,r]=lsqnonlin(f,x,y,ab)追问
追答然后??你想表达啥???像你那样放置会报错的,不信你可以试试。追问ab就是 待求参数a与b的初始值,那help中的x0不是待求参数的初始值吗,xy是往函数里传递的变量,不是上下边界追答。。。。你的函数里的a在哪??追问a(1),a(2)就是待求的参数,y-(a(2)*(1.1*sqrt(x)+a(1))./(0.5*x+52))追答那你的x,y矩阵里的值难道不是初始值吗??追问

x y 对应的数据关系就是提问中一堆一一对应的 数据。
我就是想用lsqnonlin函数进行回归,求解ab,但是老是出现矩阵维数不一致问题。有的例题就可以解。有的不行。       
追答你这题有问题吧,两个方程两个未知数,x,y只要两组数据就能求a b 了。追问这个问题领域属于拟合参数,不是解方程。
General model:
f(x) = b*(1.1*sqrt(x)+a)/(0.5*x+52)
Coefficients (with 95% confidence bounds):
a = -0.2423 (-0.3988, -0.08585)
b = 92.99 (75.48, 110.5)
Goodness of fit: R-square: 0.9482

这个是用Matlab 的cftool命令做的。是交互式.
我想知道用程序哪里出错,怎么解决?????? !!!!
有人用1stopt软件也算了,结果相同追答不好意思哈,刚刚出去了下。你可以这样,代码我写好了,你试试
function Q=myfun(a,x,y)
x =[
1.1200
1.2400
1.3000
1.1700
1.1600
0.9760
0.8090
0.8100
0.4440
0.4290];
y =[
1.5900
1.6300
1.8600
1.7700
1.7300
1.4500
1.3600
1.2000
0.9520
0.8280];
Y=(a(2)*(1.1*sqrt(x)+a(1))./(0.5*x+52));
Q=Y-y;

主程序里:
clear;clc;
x =[
1.1200
1.2400
1.3000
1.1700
1.1600
0.9760
0.8090
0.8100
0.4440
0.4290];
y =[
1.5900
1.6300
1.8600
1.7700
1.7300
1.4500
1.3600
1.2000
0.9520
0.8280];
ab=[1 2];
[a,r]=lsqnonlin(@myfun,ab,x,y)

就这样。结果是:
a =
1.5900 1.6300

r =
19.4266

满意请采纳。
我要举报
如以上问答信息为低俗、色情、不良、暴力、侵权、涉及违法等信息,可以点下面链接进行举报!
大家都在看
求好听的校斗群名字与头像(多点)
十日观察法是准确的吗?为什么我听说人死了狗
如图,在平面直角坐标系中,正△OAB的顶点A的坐
3200米 800米等于几千米
测量农田面积运用什么地理信息技术
红茶洗头好还是绿色
银河W6地板配什么胶皮和海绵最好,硬度该是多
999分之77,9999分之777,999+9999分之77+77
as a r he was lucky to get a present
简述EU NAFTA APEC 主要区别
我家住在六楼,负一层安装着暖气加压泵,一天
三义学校这个地址在什么地方,我要处理点事
o.25xo.28简便计算
借款人为什么愿意从P2P借钱
家长网上抱怨为何导致孩子被开除?
推荐资讯
小波浪卷发短发怎么做
有没有可以自己导入自己整理的生词库的背单词
抖音我希望在你的爱情里融化的糖果歌词是什么
tmd现在东莞公交真tm乱,都不到那个地方却写
快手给别人双击怎么取消,一年前双击的
目前短期理财利率是多少
咸肉芋头菜饭怎么做如何做好吃
欧洲的重要国际组织是A.欧佩克B.欧盟C.东盟D.
做海鱼放什么配菜好,
北京哪里卖香港优良品话梅加应子
欠人几万元,没钱还,该怎么办,如何跟父母开
整日无所事事却未感半分轻松什么意思
正方形一边上任一点到这个正方形两条对角线的
阴历怎么看 ?