24小时热门版块排行榜    

查看: 2525  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 299求调剂 +6 △小透明* 2026-03-17 6/300 2026-03-21 02:42 by JourneyLucky
[考研] 一志愿西南交大,求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 321求调剂 +9 何润采123 2026-03-18 11/550 2026-03-20 23:19 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +4 @taotao 2026-03-19 4/200 2026-03-20 22:14 by JourneyLucky
[考研] 350求调剂 +5 weudhdk 2026-03-19 5/250 2026-03-20 22:04 by luoyongfeng
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +22 rare12345 2026-03-18 22/1100 2026-03-20 20:39 by zhukairuo
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
[考研] 广西大学家禽遗传育种课题组2026年硕士招生(接收计算机专业调剂) +3 123阿标 2026-03-17 3/150 2026-03-20 15:58 by 飞行琦
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 6/300 2026-03-19 17:18 by fei626-918
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 10/500 2026-03-19 15:08 by 阳阳哇塞
[考研] 286求调剂 +6 lemonzzn 2026-03-16 10/500 2026-03-19 14:31 by lemonzzn
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
信息提示
请填处理意见