sim.rar_信号的时频图_时频图_画信号时频图

%%generate fec_data figure, plot constellations clc ; clear all ; close all ; modetype = 1 ; %1 16qam 2 qpsk 调制方式标志 fs = 903225 ; % 信号的采样频率 oversample = 4 ; % o&m算法需要4倍过采样 即每个符号4个采样 param_files = load('.\file_name.txt') ; %导入文本文件，读取起始、结束、中心频率、符号速率、信噪比信息 hd_srrc_rx = sfc_srrc ; %匹配滤波器 % [ uw_bin_q ] = fun_uw_bit( '.\data_out\uw\p20Q_uw.txt' ,64,15); %已找到的信号UW字库 % uw_pre_post = [uw_bin_q] ; %前后导 % uw_pre = uw_pre_post(1:end,1:40) ; %前导 % uw_post = uw_pre_post(1:end,41:60) ; %后导 % len_pre = size(uw_pre,2) ; %矩阵列数，即uw_pre长度 % len_post = size(uw_post,2) ; %矩阵列数，即uw_post长度 % res_stat_file= '.\test\uw_test_bin.csv' %res_idx = 1; res_table = {}; cnt = 0 ; plot_en = 1 ; %for file_id=1:4000 %4000条数据 for file_id = 1:13267 [sym_st,sym_et,sym_fc,sym_rate,sym_snr,str_sym_st,str_sym_et,str_sym_fc,str_sym_rate,str_snr ] = param_load( param_files, file_id); file_name = ['st(',str_sym_st,')_et(',str_sym_et,')_fc(',str_sym_fc,')_symrate(',str_sym_rate,')_snr(',str_snr,')'] ; fig_name = ['.\picture\fig', num2str(file_id), '.jpg'] ; data_obj = load(['.\data\',file_name,'.mat']) ; data_din = (data_obj.data_fram)/2^15 ; [ afc_o ,dds_o ] = z_afc( data_din,0,sym_fc/fs) ; %dds， afc_o 混频后的输出，dds_o即为混频的相干信号 hd_flt = fir_flt(sym_rate/fs) ; %低通滤波器 hd_flt为低通滤波器传输函数 sym_rate/fs为归一化后的信号带通截止频率 flt_o = step(hd_flt,afc_o) ; %阶跃响应 flt_o为低通滤波器输出 %oversample = 4; 即sps，每个符号4采样 flt_resample = resample(flt_o,sym_rate*oversample,fs) ; %resample重采样 flt_resample为重采样后的输出 过采样率：length(flt_o)/length(flt_resample) = fs/sym_rate*oversample; srrc_o = step(hd_srrc_rx,flt_resample) ; %匹配滤波 hd_srrc_rx， srrc_o为匹配滤波器后的输出 din_sof= [1;zeros(length(srrc_o)-1,1)] ; % ;英文状态下分号 :英文状态下冒号 din_sof为length(srrc_o)*1，第一位为1，剩下全为0 [ sym_rec,sym_sof,epsi_om,epsi_correct] = sfr_om(srrc_o,srrc_o,din_sof,length(din_sof)/4) ; %符号同步 o&m算法（也叫平方律非线性算法） sym_rec为符号同步后的输出 plot(sym_rec(85:length(sym_rec)-50), '.'); saveas(gcf, fig_name); end