应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » matlab数据拟合问题求助
51/1返回列表
查看: 2445  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

hudie012

荣誉版主 (文坛精英)

小兔爱吃萝卜

优秀版主优秀版主优秀版主

[求助] matlab数据拟合问题求助

已知一组数据y，x。根据这组数据要拟合一个公式，确定其中的参数，拟合公式的形式为：y=A*x^n，其中A，n是要拟合的参数。
没有接触过matlab，所以也不懂咋弄，有没有帮忙拟合一下，写一下程序。[ Last edited by hudie012 on 2012-9-6 at 15:04 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : 数据.xlsx
    • 2012-09-06 15:04:09, 27.62 K

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
todayisagift
1楼 2012-09-06 10:31:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengxing24

新虫 (初入文坛)
我先按你的建了一个curvefun3.m的函数，然后打开curve11.m可以进行拟合，x0是初始的向量，可以改动，看看效果

我另外建了一个指数形式的拟合curve22，希望能帮到你，可以自己去修改，形式很简单
matlab数据拟合问题求助
赞一下 回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : curve11.m
    • 2015-07-11 15:56:03, 14.81 K
  • 附件 2 : curve22.m
    • 2015-07-11 15:56:04, 14.82 K
  • 附件 3 : curvefun2.m
    • 2015-07-11 15:56:04, 63 bytes
  • 附件 4 : curvefun3.m
    • 2015-07-11 15:56:04, 60 bytes
8楼2015-07-11 15:59:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
hudie012: 金币+5 2012-09-10 17:38:31
xdata=1:10;
ydata=0.5*xdata.^3;   %xdata,ydata是待拟合数据，这里做个例子
p0=[1,5];   %p0是给出的拟合初值，为了与正确结果相区别，随便写了一个
p = lsqcurvefit(@(p,xdata)p(1)*xdata.^p(2),p0,xdata,ydata)  %拟合出来的p就是实际结果，实际结果为0.5，和3，与第2行的一致
赞一下(1人) 回复此楼
showmethemoney
2楼2012-09-06 10:50:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

z萝莉

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
csgt0: 金币+1, 多谢应助 2012-09-06 13:20:48
hudie012: 金币+5 2012-09-10 17:38:40
Sample TextSample Text
这个其实并不难，在研一的课上会学数值分析方法，看上去这是一个非线性函数的拟合，其实可以通过对左右取自然对数将其转化为线性函数的拟合。取对数得到：ln(y)=n*ln(a*x)=n*（ln(a)+ln(x)）=n*ln(a)+n*ln(x)这样就可以把原来的式子化为拟合数据y=a*x+b
方法如下：x1=log(x);
y1=log(y)；
m=ployfit(x1,y1,1)
然后会出现两个数（这里写成m1,m2)，则式子就是ln(y)=m1*ln(x)+m2
所以n=m1,n*ln(a)=m2.
赞一下 回复此楼
做好自己该做的事
3楼2012-09-06 11:36:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

★ ★ ★ ★ ★
hudie012: 金币+5 2012-09-10 17:38:53
拟合效果并不好，建议修改公式
CODE:
data=xlsread('数据.xlsx');
xdata=data(:,1);
ydata=data(:,2);
p0=[0.1,0.5];
options=optimset('MaxFunEvals',1000);
[p,r] = lsqcurvefit(@(p,xdata)p(1)*xdata.^p(2),p0,xdata,ydata,[],[],options);
format long
p
r
y=p(1)*xdata.^p(2);
plot(xdata,ydata,'r',xdata,y,'b')
赞一下 回复此楼
showmethemoney
4楼2012-09-06 16:26:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭