24小时热门版块排行榜    

查看: 2517  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料调剂 +4 爱擦汗的可乐冰 2026-02-28 4/200 2026-03-01 00:38 by 猫猫球alter
[基金申请] 面上模板改不了页边距吧? +5 ieewxg 2026-02-25 6/300 2026-03-01 00:10 by addressing
[考研] 307求调剂 +4 73372112 2026-02-28 6/300 2026-03-01 00:04 by ll247
[考研] 304求调剂 +3 52hz~~ 2026-02-28 5/250 2026-03-01 00:00 by 52hz~~
[考研] 317一志愿华南理工电气工程求调剂 +5 Soliloquy_Q 2026-02-28 8/400 2026-02-28 23:36 by xyx2012xyx
[考研] 272求调剂 +3 材紫有化 2026-02-28 3/150 2026-02-28 22:52 by ms629
[考研] 290求调剂 +5 材料专硕调剂; 2026-02-28 6/300 2026-02-28 21:40 by gaoxiaoniuma
[考研] 材料学调剂 +5 提神豆沙包 2026-02-28 5/250 2026-02-28 21:34 by gaoxiaoniuma
[考研] 264求调剂 +3 巴拉巴拉根556 2026-02-28 3/150 2026-02-28 21:31 by gaoxiaoniuma
[考研] 311求调剂 +8 南迦720 2026-02-28 8/400 2026-02-28 21:30 by gaoxiaoniuma
[考研] 求调剂 +4 repeatt?t 2026-02-28 4/200 2026-02-28 21:16 by gaoxiaoniuma
[考研] 298求调剂 +8 人间唯你是清欢 2026-02-28 11/550 2026-02-28 20:26 by L135790
[考研] 276求调剂 +3 路lyh123 2026-02-28 4/200 2026-02-28 19:45 by 路lyh123
[考研] 0856材料求调剂 +10 hyf hyf hyf 2026-02-28 11/550 2026-02-28 18:50 by 无际的草原
[考研] 265分求调剂不调专业和学校有行学上就 +4 礼堂丁真258 2026-02-28 6/300 2026-02-28 16:18 by 求调剂zz
[考研] 0856调剂 +3 刘梦微 2026-02-28 3/150 2026-02-28 13:22 by houyaoxu
[考研] 寻找调剂 +3 LYidhsjabdj 2026-02-28 3/150 2026-02-28 12:59 by miniwendy
[考研] 304求调剂 +5 曼殊2266 2026-02-28 6/300 2026-02-28 12:44 by 迷糊CCPs
[硕博家园] 博士自荐 +6 科研狗111 2026-02-26 9/450 2026-02-28 12:32 by seaskyy
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
请填处理意见