声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1190|回复: 3

[综合讨论] 求助!关于优化命令Lsqcuevefit的问题。多谢

[复制链接]
发表于 2006-8-26 12:58 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
clc
clear
data=xlsread('D:\wjt\资料\其他\nsxyhjffzq');
xdata1=data(:,1);
ydata1=data(:,2);
xs=xdata1(1);ffs=ydata1(1);
x0=[1000,5];
fun1=inline('ffs+(xdata1-xs)/(x1(1)*(1+ffs/x1(1))-(xdata1-xs))*(x1(2)+ffs)','x1','xdata1','ffs','xs')
[x1,resnorm1]=lsqcurvefit(fun1,x0,xdata1,ydata1);  就是这个命令的错误。请各位大佬指教指教!多谢先!

出现的问题如下:
??? Error using ==> D:\matlab6.1\toolbox\optim\private\lsqncommon
User supplied expression or inline function ==> ffs+(xdata1-xs)/(x1(1)*(1+ffs/x1(1))-(xdata1-xs))*(x1(2)+ffs)
failed with the following error:

Error using ==> inline/feval
Not enough inputs to inline function.



Error in ==> D:\matlab6.1\toolbox\optim\lsqcurvefit.m
On line 129  ==> [x,Resnorm,FVAL,EXITFLAG,OUTPUT,LAMBDA,JACOB] = ...
回复
分享到:

使用道具 举报

发表于 2006-8-26 15:28 | 显示全部楼层
其实这个函数也不是很熟悉,看了一下,你试一下,行不行。
先写个函数
function F=fun1(x1,xdata1)
ffs=xdata1(end);
xs=xdata1(1);
xdata1=xdata1(1:length(xdata1)-1);
F=ffs+(xdata1-xs)./(x1(1)*(1+ffs/x1(1))-(xdata1-xs))*(x1(2)+ffs);
保存后,在命令窗口输入,
global xs;
global ffs;
data=xlsread('D:\wjt\资料\其他\nsxyhjffzq');
xdata1=data(:,1);
ydata1=data(:,2);
xdata1=[xdata1;ydata1(1)];
x0=[1000,5];
[x1,resnorm1]=lsqcurvefit(@fun1,x0,xdata1,ydata1)
这个指令中函数fun1,好像只能传入xdata1。
没有数据验证,不知对不对。

[ 本帖最后由 buaalqq 于 2006-8-26 15:48 编辑 ]

评分

1

查看全部评分

 楼主| 发表于 2006-8-27 00:22 | 显示全部楼层

多谢您的指点!只是还有一点点问题。

>> global xs;
global ffs;
data=xlsread('D:\wjt\资料\其他\nsxyhjffzq');
xdata1=data(:,1);
ydata1=data(:,2);
xdata1=[xdata1 ydata1(1)];
x0=[1000,5];
[x1,resnorm1]=lsqcurvefit(@fun1,x0,xdata1,ydata1)

出现问题如下:
??? Error using ==> horzcat
All matrices on a row in the bracketed expression must have the
same number of rows.

十分感谢!
 楼主| 发表于 2006-8-27 00:28 | 显示全部楼层

恩。谢了。可以了。

:lol
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-9-23 05:32 , Processed in 0.054558 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表