24小时热门版块排行榜    

查看: 2527  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 317求调剂 +9 申子申申 2026-03-19 15/750 2026-03-21 17:31 by 学员8dgXkO
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 336求调剂 +5 rmc8866 2026-03-21 5/250 2026-03-21 17:24 by 学员8dgXkO
[考研] 298求调剂 +4 上岸6666@ 2026-03-20 4/200 2026-03-21 17:14 by 学员8dgXkO
[考研] 265求调剂 +12 梁梁校校 2026-03-19 14/700 2026-03-21 13:38 by lature00
[考研] 能源材料化学课题组招收硕士研究生8-10名 +5 脱颖而出 2026-03-16 15/750 2026-03-21 10:16 by 脱颖而出
[考研] 307求调剂 +3 wyyyqx 2026-03-17 3/150 2026-03-21 03:20 by JourneyLucky
[考研] 303求调剂 +5 睿08 2026-03-17 7/350 2026-03-21 03:11 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-18 3/150 2026-03-21 00:46 by JourneyLucky
[考研] 一志愿西南交大,求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 287求调剂 +7 晨昏线与星海 2026-03-19 8/400 2026-03-20 22:19 by JourneyLucky
[考博] 招收博士1-2人 +3 QGZDSYS 2026-03-18 3/150 2026-03-20 11:58 by 呱呱呱呱叫
[考研] 求调剂 +3 暗涌afhb 2026-03-16 3/150 2026-03-20 00:28 by 河南大学校友
[考研] 085601专硕,总分342求调剂,地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 考研化学学硕调剂,一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
信息提示
请填处理意见