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

登录注册

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

24小时热门版块排行榜

返回列表

13143747

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 697.4
散金: 300
帖子: 109
在线: 54.4小时
虫号: 2542177
注册: 2013-07-11
性别: GG
专业: 系统科学与系统工程

[求助] 信号分析中有两个for循环太浪费时间了求帮助怎么减少程序运行时间

我的想法：1：运用MATLAB GPU并行运算运算来把for循环放到gpu中计算；
2：将两个for循环改成矩阵运算。
但是没成功

close all;clear all;
fs=20000;%采样频率
f=fs/10;%调制频率
N=2000;%采样点数
t=0:1/fs:N/fs-1/fs;

fm=50;%模拟故障频率
b=fm*3000;%衰减率
Tm2=1/fm;%冲击周期
nn=Tm2*floor(t/Tm2);
x1=0.1*exp(-b*(t-nn).^2).*sin(2*pi*f*t);%模拟故障信号，冲击信号
x11=abs(hilbert(x1));%包络解调信号，参考故障诊断里的共振解调方法

%原始信号及频谱
figure;subplot(211);plot(t,x1);
subplot(212);hua_fft(x1,fs,2);

%输入噪声
%load 2000noisebearing.mat;

%% raw noise signal 含噪噪声及频谱

x33=x1+(0.15/0.6)*x2;
figure;subplot(211);plot(t,x33);
subplot(212);hua_fft(x33,fs,2);

Wn=[(1000)/(fs/2) (3000)/(fs/2)]; %对含噪信号先进行带通滤波，再解调
[b, a]=butter(3,Wn);
x333=filter(b,a,x33);
x3=abs(hilbert(x333));

figure;subplot(211);plot(t,x3);
subplot(212);hua_fft(x3,fs,2,0,1000);

S2=x3;%随机共振输入信号

%% Parameters tuning of 2bistable potential
%通过寻找信噪比最优来寻优
Max_SNR2=-50;
index2=1;

      for h=5:5:1000
         for r=0.001:0.001:0.06
         % twosr为随机共振算法
         s6=twosr(1,1,h/fs,r,S2);
         s6=s6-mean(s6);

         y=fft(s6,N);
         pyy=y.*conj(y)/N;

         P_S_1=pyy(round(fm*N/fs)+1);
         P_N_1=sum(pyy(1:N/2))-P_S_1;
         SNR_B2(int32(h/5),int32(r*1000))=10*log10(P_S_1/P_N_1);
         if 10*log10(P_S_1/P_N_1)>Max_SNR2
            Max_SNR2=10*log10(P_S_1/P_N_1);
            Maxh=h;
            Maxr=r;
         end

         end
      end

%得到SNR最优的随机共振输出信号
s8=twosr(1,1,Maxh/fs,Maxr,S2);
figure;subplot(211);plot(t,s8);
subplot(212);hua_fft(s8,fs,2,0,1500);

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : 2000noisebearing.mat

2015-10-08 21:43:45, 15.19 K

附件 2 : hua_fft.m

2015-10-08 21:43:45, 2.01 K

附件 3 : twosr.m

2015-10-08 21:43:45, 682 bytes

» 猜你喜欢

孩子确诊有中度注意力缺陷已经有12人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复
2026博士申请-功能高分子，水凝胶方向已经有6人回复
论文投稿，期刊推荐已经有4人回复
硕士和导师闹得不愉快已经有13人回复
请问2026国家基金面上项目会启动申2停1吗已经有5人回复
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

信号分析中有两个for循环太浪费时间了求帮助怎么减少程序运行时间已经有0人回复

1楼2015-10-08 21:41:32

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

Scorpio_0

银虫 (著名写手)

应助: 25 (小学生)
金币: 2558.7
散金: 110
红花: 16
沙发: 1
帖子: 1282
在线: 179.4小时
虫号: 3932734
注册: 2015-06-20
性别: GG
专业: 计算机科学

坐故障诊断的，哈哈

发自小木虫Android客户端

赞一下

回复此楼

Practice makes perfect

2楼2015-10-08 21:45:04

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

13143747

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 697.4
散金: 300
帖子: 109
在线: 54.4小时
虫号: 2542177
注册: 2013-07-11
性别: GG
专业: 系统科学与系统工程

引用回帖:

2楼: Originally posted by Scorpio_0 at 2015-10-08 21:45:04
坐故障诊断的，哈哈

对的

回复此楼

3楼2015-10-08 21:47:10

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

Scorpio_0

银虫 (著名写手)

应助: 25 (小学生)
金币: 2558.7
散金: 110
红花: 16
沙发: 1
帖子: 1282
在线: 179.4小时
虫号: 3932734
注册: 2015-06-20
性别: GG
专业: 计算机科学

引用回帖:

3楼: Originally posted by 13143747 at 2015-10-08 21:47:10
对的

告诉你，不要追求实效性，无意中，真想时效就和c++混合编程，matlab并行计算不行，数据集太少的时候还会更慢。还有就是并行计算很多要求，比如同层数据之间没有递归计算。

发自小木虫Android客户端

赞一下

回复此楼

Practice makes perfect

4楼2015-10-08 21:50:37

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

13143747

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 697.4
散金: 300
帖子: 109
在线: 54.4小时
虫号: 2542177
注册: 2013-07-11
性别: GG
专业: 系统科学与系统工程

引用回帖:

4楼: Originally posted by Scorpio_0 at 2015-10-08 21:50:37
告诉你，不要追求实效性，无意中，真想时效就和c++混合编程，matlab并行计算不行，数据集太少的时候还会更慢。还有就是并行计算很多要求，比如同层数据之间没有递归计算。
...

这个程序跑完20秒，我扩大了范围，跑了2000多秒……

发自小木虫Android客户端

赞一下

回复此楼

5楼2015-10-08 22:13:44

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