| 查看: 515 | 回复: 2 | ||
[求助]
噪声分布的程序
|
|
1、如果知道一组噪声的时序，怎么确定它服从什么分布，好像是有一个这样的程序。金币不够多，本来写了给100个，发了以后才发现一共就60多个了，还改不了了，只能给50个了，如果以后攒够了再补上，谢谢啦。 [ Last edited by 李三石 on 2013-5-28 at 12:17 ] |
回复此楼» 猜你喜欢
三无产品还有机会吗
已经有6人回复
-
垃圾破二本职称评审标准
已经有7人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有7人回复
-
博士申请都是内定的吗？
已经有14人回复
-
谈谈两天一夜的“延安行”
已经有13人回复
-
氨基封端PDMS和HDI反应快速固化
已经有11人回复
-
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有11人回复
-
论文投稿求助
已经有4人回复
-
Applied Surface Science 这个期刊。有哪位虫友投过的能把word模板发给我参考一下嘛
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- MATLAB中产生高斯白噪声，涉及到awgn和wgn函数
已经有3人回复
- Matlab中添加乘性噪声的问题
已经有5人回复
- 求助Matlab图像加入噪声的问题
已经有1人回复
- 重金求助Matlab图像加噪声问题
已经有8人回复
- 【求助】如何用Matlab产生高斯白噪声
已经有6人回复
2楼2013-06-07 10:07:42
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
李三石: 金币+100, ★★★★★最佳答案 2013-10-23 18:43:10
ben_ladeng: 请认真阅读此帖，http://emuch.net/bbs/viewthread.php?tid=3711903。 并与两日内解释您和ID 李三石 之间的关系。 2013-10-24 09:51:17
nono2009: 金币-50, 马甲存档, 马甲转移金币行为：http://emuch.net/bbs/viewthread.php?tid=1996653 2013-10-28 14:44:22
nono2009: 金币-80, 扣除转移金币 2013-10-28 14:44:40
李三石: 金币+100, ★★★★★最佳答案 2013-10-23 18:43:10
ben_ladeng: 请认真阅读此帖，http://emuch.net/bbs/viewthread.php?tid=3711903。 并与两日内解释您和ID 李三石 之间的关系。 2013-10-24 09:51:17
nono2009: 金币-50, 马甲存档, 马甲转移金币行为：http://emuch.net/bbs/viewthread.php?tid=1996653 2013-10-28 14:44:22
nono2009: 金币-80, 扣除转移金币 2013-10-28 14:44:40
|
function [shang]=jss(xdate) m=2; n=length(xdate); r=0.2*std(xdate); cr=[]; gn=1; gnmax=m; while gn<=gnmax x2m=zeros(n-m+1,m);%存放变换后的向量 d=zeros(n-m+1,n-m);% 存放距离结果的矩阵 cr1=zeros(1,n-m+1);%存放 k=1; for i=1:n-m+1 for j=1:m x2m(i,j)=xdate(i+j-1); end end x2m; for i=1:n-m+1 for j=1:n-m+1 if i~=j d(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离 k=k+1; end end k=1; end d; for i=1:n-m+1 [k,l]=size(find(d(i,<r));%将比R小的个数传送给L cr1(1,i)=l; end cr1; cr1=(1/(n-m))*cr1; sum1=0; for i=1:n-m+1 if cr1(i)~=0 sum1=sum1+log(cr1(i)); end end cr1=1/(n-m+1)*sum1; cr(1,gn)=cr1; gn=gn+1; m=m+1; end cr; shang=cr(1,1)-cr(1,2); function [shang]=ybs(xdate) m=2; n=length(xdate); r=0.2*std(xdate); cr=[]; gn=1; gnmax=m; while gn<=gnmax x2m=zeros(n-m+1,m);%存放变换后的向量 d=zeros(n-m+1,n-m);% 存放距离结果的矩阵 cr1=zeros(1,n-m+1);%存放 k=1; for i=1:n-m+1 for j=1:m x2m(i,j)=xdate(i+j-1); end end x2m; for i=1:n-m+1 for j=1:n-m+1 if i~=j d(i,k)=max(abs(x2m(i,-x2m(j,));%计算各个元素和响应元素的距离 k=k+1; end end k=1; end d; for i=1:n-m+1 [k,l]=size(find(d(i,<r));%将比R小的个数传送给L cr1(1,i)=l; end cr1; cr1=(1/(n-m))*cr1; sum1=0; for i=1:n-m+1 sum1=sum1+cr1(i); end end cr1=1/(n-m+1)*sum1; cr(1,gn)=cr1; gn=gn+1; m=m+1; end cr; shang=-log(cr(1,1)/cr(1,2)); |
3楼2013-10-23 18:34:56