| 查看: 2910 | 回复: 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
» 猜你喜欢
93BebMhtakh前后11位开头都是大写
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有86人回复
-
感觉是下周放榜了
已经有11人回复
-
filecode=后面第一个是大写字母
已经有10人回复
-
时间戳变了，能看出什么问题？
已经有4人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 信号分析中有两个for循环太浪费时间了求帮助怎么减少程序运行时间
已经有0人回复
dk1013
木虫 (正式写手)
- 应助: 61 (初中生)
- 金币: 4285.5
- 散金: 246
- 红花: 13
- 帖子: 988
- 在线: 3204.6小时
- 虫号: 980198
- 注册: 2010-03-23
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
9楼2015-10-09 09:10:03
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
