帮忙检查一行代码，谢谢

zanjero · 发表于 2006-12-22 15:54

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

Qo,Qc为１列n行的序列,欲求DC为（１x１）的矩阵
DC=1-E(qc(i)-qo(i))^2/E(qo(i)-mean(qo))^2
其中E为求和-_-
dc=1-sum((qc-qo).^2)/sum(qo-mean(qo))^2
上面这样写不对，请帮改一下

ps:qo(i)-mean(qo)就是距平吧，matlab中有没有这样的函数呢

dy123 · 发表于 2006-12-23 09:50

>> qc=[1 2 3]';
>> qo=[4 5 6]';
>> d=1-sum((qc-qo).^2)/sum((qo-mean(qo)).^2)

d =

-12.5000

lxq · 发表于 2006-12-23 10:51

不太明白是用于哪方面的?

zanjero · 发表于 2006-12-23 16:24

谢谢dy123

原帖由 lxq 于 2006-12-23 10:51 发表
不太明白是用于哪方面的?

呵呵，专业比较冷门，水文上算效率系数的

zanjero · 发表于 2006-12-23 16:45

假设qc,qo都有5列，下面的代码怎么简化啊？
for i=1:5
mean_qo=mean(qo(:,i))
d(i)=1-sum((qc-qo).^2)/sum((qo-mean_qo).^2)
next

账号		自动登录	找回密码
密码			我要加入

[编程技巧] 帮忙检查一行代码，谢谢

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

看看我做的

评分