作者wmin (Gary)
看板MATLAB
标题[问题] 隐函数拟合後 求相对应y值
时间Fri Feb 10 03:04:24 2012
发上一篇问的问题後
发觉自己的函数是个隐函数Orz...
参考google大神上对於隐函数的拟合问题
终於解决了~
不过我还有一个问题...
xlsFile='data.xls';
t(1,:)=xlsread(xlsFile,'sheet1','A1:A600')';
y(1,:)=xlsread(xlsFile,'sheet1','B1:B600')';
tt=t';
yy=y';
F=@(p,x) ((log((x(:,2)*(p(2)-p(3)))./(p(2)*(x(:,2)-p(3))))
-p(3)*(x(:,2).^-1-p(2)^-1))./(p(1)*p(3)^2))-x(:,1);
x=[tt yy];
p0=[100000 0.00012 0.00004];
p=nlinfit(x,zeros(size(x,1),1),F,p0)
plot(x(:,1),x(:,2),'ro');
hold on;
ezplot(@(t,y)F(p,[t y]),[0,43200,0.00018,0.000039]);
format long;
就是到这里一切都没问题
但是我想求 拟合得到参数p後
带入x值 求出相对应的y值
这边不知道有人可以给我建议该如何求吗?
我有google到一个例子
就是假设f(x,y) = exp(x/y)-x^sin(y)=0
给定x值 求相对应y可利用
Y = @(x) fzero(@(y) exp(x/y)-x^sin(y),1);
x=1.5时之y值即为
Y(1.5)
但是我的例子要怎麽套用上去
我就搞了好久...
有人可以给我一点hint吗?
先感谢各位了Q_Q
--
年轻就是要挥霍青春热血一下!
什麽事都要精打细算过後才去做的话
也许都忘了一开始是为了什麽...
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.115.150.220
1F：推 maplefog:其实图都画出来了 可以直接从图上把所有点的座标读取出来 02/11 00:22
2F：→ maplefog:再用内插法就可以了 02/11 00:23
3F：→ maplefog:座标取出来和内插法就靠自己研究了 02/11 00:28
谢谢a大提供的意见 不过我用fzero其他方法解出来了XD
※ 编辑: wmin 来自: 140.115.150.220 (02/19 18:36)