版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

当前主题已经存档。

飞扬2282

荣誉版主 (著名写手)

应助: 0 (幼儿园)
贵宾: 2.315
金币: 10773.4
散金: 1185
红花: 17
帖子: 1897
在线: 705.2小时
虫号: 661639
注册: 2008-11-25
性别: GG
专业: 核技术及其应用
管辖: 学术会议

[交流] 【求助】【30金币求助】Matlab编程代码？

已知N个点的坐标（ui，vi）, 其中i=1，2，3...N;
由下列公式计算点p0,p1,p2,p3,p4的值。

p0*(ui)*(ui)-2*p1*ui-2*p2*vi+2*p3*ui*vi+p4=-(vi)*(vi)

matlab实现代码？

回复此楼

» 猜你喜欢

假如你的研究生提出不合理要求已经有12人回复
实验室接单子已经有7人回复
全日制（定向）博士已经有5人回复
萌生出自己或许不适合搞科研的想法，现在跑or等等看？已经有4人回复
Materials Today Chemistry审稿周期已经有4人回复
参与限项已经有3人回复
对氯苯硼酸纯化已经有3人回复
所感已经有4人回复
要不要辞职读博？已经有7人回复
北核录用已经有3人回复

1楼 2009-04-10 16:24:13

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

dnp

荣誉版主 (知名作家)

小木虫浪子

应助: 74 (初中生)
贵宾: 18.476
金币: 18799.9
散金: 2040
红花: 77
沙发: 6
帖子: 7304
在线: 593.8小时
虫号: 437184
注册: 2007-08-26
性别: GG
专业: 生物信息学
管辖: 有机交流

★ ★ ★ ★ ★
飞扬2282(金币+5,VIP+0):具体问题？ 4-10 17:16

和前面那个帖子变化一下即可，仍然是采用最小二乘法来计算：

x1 = u(1:N).*u(1:N);
x2 = -2*u(1:N);
x3 = -2*v(1:N);
x4 = 2*u(1:N)*v(1:N);
x5 = ones(N,1);
y = -v(1:N).*v(1:N);

X = [x1 x2 x3 x4 x5];

b = inv(X'*X)*X'*y; % or b = pinv(X)*y;

注意：x1,x2,x3,x4,x5必须都为列向量，y也必须得为列向量；

得到的b值中的各个元素就是p0,p1,p2,p3,p4的值。

赞一下(6人)

回复此楼

What would Jesus do?

2楼2009-04-10 16:59:20

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

飞扬2282

荣誉版主 (著名写手)

应助: 0 (幼儿园)
贵宾: 2.315
金币: 10773.4
散金: 1185
红花: 17
帖子: 1897
在线: 705.2小时
虫号: 661639
注册: 2008-11-25
性别: GG
专业: 核技术及其应用
管辖: 学术会议

引用回帖:

Originally posted by dnp at 2009-4-10 16:59:
和前面那个帖子变化一下即可，仍然是采用最小二乘法来计算：

x1 = u(1:N).*u(1:N);
x2 = -2*u(1:N);
x3 = -2*v(1:N);
x4 = 2*u(1:N)*v(1:N);
x5 = ones(N,1);
y[/ ...

写成如下，执行不了啊，有什么问题吗？

N=6;
u=[1 3 3 4 5 6]';
v=[3 4 5 7 8 10]';
x1 = u(1:N).*u(1:N);
x2 = -2*u(1:N);
x3 = -2*v(1:N);
x4 = 2*u(1:N).*v(1:N);
x5 = ones(N,1);
y = -v(1:N).*v(1:N);
X = [x1 x2 x3 x4 x5]';
b = inv(X'*X)*X'*y

赞一下

回复此楼

3楼2009-04-10 17:16:15

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

dnp

荣誉版主 (知名作家)

小木虫浪子

应助: 74 (初中生)
贵宾: 18.476
金币: 18799.9
散金: 2040
红花: 77
沙发: 6
帖子: 7304
在线: 593.8小时
虫号: 437184
注册: 2007-08-26
性别: GG
专业: 生物信息学
管辖: 有机交流

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
飞扬2282(金币+25,VIP+0):thanks~~ 4-10 18:53

引用回帖:

Originally posted by 飞扬2282 at 2009-4-10 17:16:

写成如下，执行不了啊，有什么问题吗？

N=6;
u=[1 3 3 4 5 6]';
v=[3 4 5 7 8 10]';
x1 = u(1:N).*u(1:N);
x2 = -2*u(1:N);
x3 = -2*v(1:N);
x4 = 2*u(1:N).*v(1:N);
x5 = ones(N,1);
y = -v(1: ...

X = [x1 x2 x3 x4 x5]';

把这句代码改成：

X = [x1 x2 x3 x4 x5];

即可~~

前面的u, v都是列向量，所以X应该是6行5列的矩阵。
注意：
X的行数和y的行数必须一致~~X不能转置。

引用回帖:

>> N=6;
u=[1 3 3 4 5 6]';
v=[3 4 5 7 8 10]';
x1 = u(1:N).*u(1:N);
x2 = -2*u(1:N);
x3 = -2*v(1:N);
x4 = 2*u(1:N).*v(1:N);
x5 = ones(N,1);
y = -v(1:N).*v(1:N);
X = [x1 x2 x3 x4 x5];
b = inv(X'*X)*X'*y

b =

2.3946
-0.6497
0.0201
-1.5886
-2.9967

[ Last edited by dnp on 2009-4-10 at 17:28 ]

赞一下(5人)

回复此楼

What would Jesus do?

4楼2009-04-10 17:22:00

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

embrace5758

应助: (幼儿园)
在线:
虫号: 705680

学了~~~

回复此楼

5楼2009-04-13 10:59:58

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

相关版块跳转我要订阅楼主飞扬2282 的主题更新

返回列表