| 查看: 2896 | 回复: 9 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
13143747银虫 (小有名气)
|
[求助]
信号分析中有两个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
- 附件 2 : hua_fft.m
- 附件 3 : twosr.m
2015-10-08 21:43:45, 15.19 K
2015-10-08 21:43:45, 2.01 K
2015-10-08 21:43:45, 682 bytes
» 猜你喜欢
咱们一起用铁证分析2026国家社科基金中标与否
已经有17人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有47人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有6人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有11人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有4人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有37人回复
-
好奇怪的filecode
已经有7人回复
-
filecode
已经有9人回复
-
请教兼职经验
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 信号分析中有两个for循环太浪费时间了求帮助怎么减少程序运行时间
已经有0人回复
13143747
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 697.4
- 散金: 300
- 帖子: 109
- 在线: 54.4小时
- 虫号: 2542177
- 注册: 2013-07-11
- 性别: GG
- 专业: 系统科学与系统工程
10楼2015-10-09 09:43:26
Scorpio_0
银虫 (著名写手)
- 应助: 25 (小学生)
- 金币: 2558.7
- 散金: 110
- 红花: 16
- 沙发: 1
- 帖子: 1282
- 在线: 179.4小时
- 虫号: 3932734
- 注册: 2015-06-20
- 性别: GG
- 专业: 计算机科学
2楼2015-10-08 21:45:04
13143747
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 697.4
- 散金: 300
- 帖子: 109
- 在线: 54.4小时
- 虫号: 2542177
- 注册: 2013-07-11
- 性别: GG
- 专业: 系统科学与系统工程
3楼2015-10-08 21:47:10
Scorpio_0
银虫 (著名写手)
- 应助: 25 (小学生)
- 金币: 2558.7
- 散金: 110
- 红花: 16
- 沙发: 1
- 帖子: 1282
- 在线: 179.4小时
- 虫号: 3932734
- 注册: 2015-06-20
- 性别: GG
- 专业: 计算机科学
|
告诉你，不要追求实效性，无意中，真想时效就和c++混合编程，matlab并行计算不行，数据集太少的时候还会更慢。还有就是并行计算很多要求，比如同层数据之间没有递归计算。 发自小木虫Android客户端 |
4楼2015-10-08 21:50:37
