24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1687  |  回复: 6

ww汪红

新虫 (小有名气)

[求助] 求大神解释~~F=@(p,x)p(1)*x(:,1).^2+p(2)*x(:,1).*x(:,2)+p(3)*x(:,2).^2+p(4); 已有1人参与

F=@(p,x)p(1)*x(:,1).^2+p(2)*x(:,1).*x(:,2)+p(3)*x(:,2).^2+p(4);
% 离散数据点
x=load('G:\Users\Daniel\Desktop\Total.txt');
p0=[1 1 1 1];
warning off
% 拟合系数,最小二乘方法
p=nlinfit(x,zeros(size(x,1),1),F,p0);
plot(x(:,1),x(:,2),'r.');
axis square;
hold on;
xmin=min(x(:,1));
xmax=max(x(:,1));
ymin=min(x(:,2));
ymax=max(x(:,2));
% 作图
ezplot(@(x,y)F(p,[x,y]),[-1+xmin,1+xmax,-1+ymin,1+ymax]);
title('曲线拟合');
legend('样本点','拟合曲线')
set(gca,'xtick',[-1:0.2:1]);
set(gca,'ytick',[-1:0.2:1]);
xlim([-1 1]);
ylim([-1 1]);
line([-1 1],[0 0]);
line([0 0],[-1 1]);

麻烦请问能不能设置拟合出来的椭圆长轴沿y方向,短轴沿x方向?
还有这个公式F=@(p,x)p(1)*x(:,1).^2+p(2)*x(:,1).*x(:,2)+p(3)*x(:,2).^2+p(4)

哪个代表长轴哪个代表短轴呀~?小白一个跪谢大神
回复此楼

» 猜你喜欢

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

guo_zy_123

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
F=@(p,x)p(1)*x(:,1).^2+p(2)*x(:,2).^2+p(3)%匿名函数写成这个

%P的初值只要
p0 = [1 1 1];

其他部分不变

上次给你程序了,长短轴,根据算出来的长短轴,比较大小就行了

见你自己上一个帖子
老古董
2楼2015-12-27 09:42:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guo_zy_123

木虫 (正式写手)

【答案】应助回帖

长短轴公式我就不列了吧

这个高中生的东西,自己查一下下

x^2/(a^2)+y^2/(b^2) = 1

把求解出的p转换成上面的标准格式

ab都能求,然后比较大小就行了。。。。
老古董
3楼2015-12-27 09:44:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guo_zy_123

木虫 (正式写手)

【答案】应助回帖

这次看了你的数据,
以前没有看过

感觉你的数据比较适合分段拟合

比较明显的是几个环状数据

绿色的线是所有数据用最小二乘法拟合的,

蓝色的线是分了几个段拟合的

蓝色拟合的基础是数据根据到原点的具体把数据手动分几类,然后就能做了
求大神解释~~F=@(p,x)p(1)*x(:,1).^2+p(2)*x(:,1).*x(:,2)+p(3)*x(:,2).^2+p(4);
分段拟合.jpg

老古董
4楼2015-12-27 09:58:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guo_zy_123

木虫 (正式写手)

金币谢谢
金币谢谢
金币谢谢

重要的事情说3遍
老古董
5楼2015-12-27 10:00:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ww汪红

新虫 (小有名气)

引用回帖:
3楼: Originally posted by guo_zy_123 at 2015-12-27 09:44:23
长短轴公式我就不列了吧

这个高中生的东西,自己查一下下

x^2/(a^2)+y^2/(b^2) = 1

把求解出的p转换成上面的标准格式

ab都能求,然后比较大小就行了。。。。

p的解在哪可以看到?
麻烦举例说一下,怎么样利用p值求解 ab
谢谢~!
6楼2015-12-27 15:28:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guo_zy_123

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
ww汪红: 金币+30, ★★★★★最佳答案, 谢谢谢谢~!! 2015-12-27 19:41:56
引用回帖:
6楼: Originally posted by ww汪红 at 2015-12-27 15:28:29
p的解在哪可以看到?
麻烦举例说一下,怎么样利用p值求解 ab
谢谢~!...

p=nlinfit(x,zeros(size(x,1),1),F,p0);

程序里这一行是,还有你这也太不动脑了吧。。。。。。

F=@(p,x)p(1)*x(:,1).^2+p(2)*x(:,2).^2+p(3)%匿名函数写成这个

则椭圆方程就是

(x^2)*p(1)+(y^2)*p(2) = - p(3)

长短轴自己对应分析吧。。。。。。。
老古董
7楼2015-12-27 19:24:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ww汪红 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿085502,267分求调剂 +14 再忙也要吃饭啊 2026-04-08 15/750 2026-04-12 22:33 by fqwang
[考研] 电气专硕320求调剂 +6 小麻子111 2026-04-10 6/300 2026-04-12 10:54 by lemon6009
[考研] 305求调剂 +6 77Qi 2026-04-07 6/300 2026-04-12 02:30 by 秋豆菜芽
[考研] 316求调剂 +5 想读研究生( ?∵ 2026-04-07 5/250 2026-04-12 00:43 by 蓝云思雨
[考研] 电子信息279求调剂,有书读就行 +8 wwwooden 2026-04-08 11/550 2026-04-11 20:22 by cq2548
[考研] 085501机械专硕 302分 不挑专业求调剂 +7 汪某. 2026-04-09 7/350 2026-04-11 14:37 by luhong1990
[考研] 302分求调剂 +9 凡语祈愿 2026-04-08 10/500 2026-04-10 23:26 by 314126402
[考研] 083200 305分 求二轮调剂 不接受跨专业 +9 Claireyyyy 2026-04-09 10/500 2026-04-10 21:21 by Claireyyyy
[考研] 本科西工大 324求调剂 +4 wysyjs25 2026-04-10 4/200 2026-04-10 20:00 by 来看流星雨10
[考研] 265求调剂 +12 风说她早忘了 2026-04-10 13/650 2026-04-10 18:56 by chemisry
[考研] 0702物理学学硕299求调剂 +6 祁柒连 2026-04-06 6/300 2026-04-10 11:10 by Roomoo
[考研] 生物学求调剂 一志愿沪9,326分 +7 刘墨墨 2026-04-06 7/350 2026-04-10 08:11 by kangsm
[考研] 考研调剂 +13 冰冰,,, 2026-04-07 13/650 2026-04-09 17:01 by Lilly_Li
[考研] 265求调剂 +4 风说她早忘了 2026-04-07 4/200 2026-04-09 13:59 by only周
[考研] 一志愿华南师范大学0702物理学305调剂 +4 念常安 2026-04-07 6/300 2026-04-08 22:53 by bljnqdcc
[考研] 296求调剂 +3 汪!?! 2026-04-08 3/150 2026-04-08 22:00 by zhouyuwinner
[考研] 266调剂 +8 daya sun 2026-04-07 9/450 2026-04-08 20:27 by yutian743
[考研] 求调剂,现在还能填的 +3 上岸小莹加油 2026-04-08 3/150 2026-04-08 14:30 by zhq0425
[考研] 一志愿西南090202求调剂 +4 在线求有学上 2026-04-07 4/200 2026-04-07 19:47 by biomichael
[考研] 0854求调剂 +9 亨氏番茄沙司 2026-04-06 10/500 2026-04-07 14:37 by shdgaomin
信息提示
请填处理意见