| 查看: 600 | 回复: 2 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
噪声分布的程序
|
|||
|
1、如果知道一组噪声的时序，怎么确定它服从什么分布，好像是有一个这样的程序。金币不够多，本来写了给100个，发了以后才发现一共就60多个了，还改不了了，只能给50个了，如果以后攒够了再补上，谢谢啦。 [ Last edited by 李三石 on 2013-5-28 at 12:17 ] |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
时间戳又变了8-15
已经有10人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有57人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有9人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
关于Filecode分析方法
已经有13人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有8人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- MATLAB中产生高斯白噪声，涉及到awgn和wgn函数
已经有3人回复
- Matlab中添加乘性噪声的问题
已经有5人回复
- 求助Matlab图像加入噪声的问题
已经有1人回复
- 重金求助Matlab图像加噪声问题
已经有8人回复
- 【求助】如何用Matlab产生高斯白噪声
已经有6人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
李三石: 金币+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
2楼2013-06-07 10:07:42
