应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 求大神解释~~F=@(p,x)p(1)*x(:,1).^2+p(2)*x(:,1).*x(:,2)+p(3)*x(:,2).^2+p(4);
51/1返回列表
查看: 1686  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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)

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

» 猜你喜欢
1楼 2015-12-26 21:51:14
已阅   回复此楼   关注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的回帖
查看全部 7 个回答

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

木虫 (正式写手)

【答案】应助回帖

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

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

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

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

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

蓝色拟合的基础是数据根据到原点的具体把数据手动分几类，然后就能做了
求大神解释~~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的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 279求调剂 +9 张番茄不炒蛋 2026-04-11 9/450 2026-04-12 22:22 by paopaotu326
[考研] 326求调剂 +6 Shansyn 2026-04-10 6/300 2026-04-12 09:46 by hammer3
[考研] 307求调剂 +10 tzq94092 2026-04-10 10/500 2026-04-12 08:18 by wise999
[考研] 化学工程调剂289 +44 yang婷 2026-04-07 50/2500 2026-04-12 02:36 by 秋豆菜芽
[考研] 293求调剂 +8 勇远库爱314 2026-04-06 8/400 2026-04-11 20:25 by 蓝云思雨
[考研] 087100初试311求调剂 +4 任雅琴 2026-04-09 4/200 2026-04-11 10:33 by zhq0425
[考研] 生物学调剂 +8 小冉要努力 2026-04-10 9/450 2026-04-11 10:22 by wwj2530616
[考研] 中药学调剂 初试324 +4 洋甘菊、 2026-04-10 6/300 2026-04-11 09:41 by gong120082
[考研] 085402通信工程调剂，有4项学科竞赛国奖（电赛国二），硕士研究生调剂自荐信。 +5 m永o不v言o弃m 2026-04-09 5/250 2026-04-11 09:33 by zhq0425
[考研] 342电子信息专硕求调剂 +9 你让我怎么荔枝 2026-04-10 10/500 2026-04-11 08:33 by zhq0425
[考研] 085404 298分求调剂 +10 呼啦呼啦呼呼呼 2026-04-10 11/550 2026-04-10 16:44 by wangy0907
[考研] 机械还有还有名额吗？太难了 +6 笑笑袁 2026-04-10 6/300 2026-04-10 11:54 by 高维春
[考研] 已调剂 +18 柴郡猫_ 2026-04-09 19/950 2026-04-09 22:10 by 柴郡猫_
[考研] 材料专硕(0856) 339分求调剂 +9 哈哈哈鹅哈哈哈 2026-04-09 10/500 2026-04-09 20:01 by Orcid
[考研] 考研调剂 +13 冰冰，，， 2026-04-07 13/650 2026-04-09 17:01 by Lilly_Li
[考研] 材料307分求大佬组收留 +17 Hll胡 2026-04-07 17/850 2026-04-09 10:53 by liuhuiying09
[考研] 327求调剂 +10 Xxjc1107. 2026-04-06 11/550 2026-04-09 01:21 by lature00
[考研] 考研求调剂 +4 雯??? 2026-04-08 4/200 2026-04-08 21:44 by 土木硕士招生
[考研] 材料考研求调剂总分280 +30 mkjlz1 2026-04-06 35/1750 2026-04-08 21:25 by cyh—315
[考研] 313求调剂 +3 十六拾陆 2026-04-07 3/150 2026-04-07 23:20 by lbsjt
信息提示
关闭
请填处理意见
关闭