24小时热门版块排行榜    

查看: 2526  |  回复: 5

星尘2516

木虫 (小有名气)

[求助] Matlab非线性拟合(y=sqrt(1-(x-a).^2))

x=-5:0.01:5;
y=sqrt(1-(x).^2);
plot(x,y,'bo','Markersize',5);
x=x';y=y';

s=fitoptions('Method','NonlinearLeastSquares','StartPoint',[20 0 0.8]);
f=fittype('abs(a0.*(1-(x-a1).^2).^0.5+a2)','options',s);
cfun=fit(x,y,f);
yf=abs(cfun.a0.*(1-(x-cfun.a1).^2).^0.5+cfun.a2);
hold on;
plot(x,yf,'r','LineWidth',2);
以上是我的程序,我打算对y=a*sqrt(1-(x-b).^2)+c,这个类型的函数进行拟合,但是总是不能完成。我用上面的程序拟合高斯、洛伦兹等函数就可以

[ Last edited by 星尘2516 on 2013-5-22 at 22:12 ]
回复此楼

» 收录本帖的淘帖专辑推荐

matlab学习

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
是否是公式中“1-(x-b).^2”可能小于0而导致开方出错?
2楼2013-05-23 08:44:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

星尘2516

木虫 (小有名气)

引用回帖:
2楼: Originally posted by dingd at 2013-05-23 08:44:51
是否是公式中“1-(x-b).^2”可能小于0而导致开方出错?

有可能,不知道是不是由于二次根号的存在,我试着用cftool也不行
3楼2013-05-23 09:06:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

加个“abs()”试试。
4楼2013-05-23 09:26:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

星尘2516

木虫 (小有名气)

引用回帖:
4楼: Originally posted by dingd at 2013-05-23 09:26:39
加个“abs()”试试。

不行的,还是不能拟合
5楼2013-05-24 11:07:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
星尘2516: 金币+10, ★★★很有帮助, 我的这个拟合有点问题,修改后已经好了,谢谢回复 2013-05-26 21:26:44
产生数据的代码:“y=sqrt(1-(x).^2);”都不对吧,x的范围在[-5,5],如何计算“sqrt(1-(x).^2)”
6楼2013-05-24 14:34:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 星尘2516 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 336求调剂 +3 rmc8866 2026-03-21 3/150 2026-03-21 11:16 by 求调剂zz
[考研] 070300化学319求调剂 +7 锦鲤0909 2026-03-17 7/350 2026-03-21 03:46 by JourneyLucky
[考研] 324分 085600材料化工求调剂 +4 llllkkkhh 2026-03-18 4/200 2026-03-21 01:24 by JourneyLucky
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +4 葵梓卫队 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[考研] 材料与化工 322求调剂 +4 然11 2026-03-19 4/200 2026-03-20 22:12 by luoyongfeng
[考研] 一志愿中南化学(0703)总分337求调剂 +8 niko- 2026-03-19 9/450 2026-03-20 21:57 by luoyongfeng
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] 295复试调剂 +8 简木ChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
[考研] 319求调剂 +3 小力气珂珂 2026-03-20 3/150 2026-03-20 19:47 by JourneyLucky
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 281求调剂(0805) +14 烟汐忆海 2026-03-16 25/1250 2026-03-20 15:47 by yuncha
[考研] 0856调剂,是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
[考研] 298-一志愿中国农业大学-求调剂 +9 手机用户 2026-03-17 9/450 2026-03-20 14:24 by 无懈可击111
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考研] 本科郑州大学物理学院,一志愿华科070200学硕,346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 材料,纺织,生物(0856、0710),化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
信息提示
请填处理意见