模拟连续波CW信号生成block的work函数

东方耀

模拟连续波CW信号生成block的work函数

1. #include <iostream>
   
2. #include <complex>
   
3. #include <vector>
   
4. #include <unistd.h>
   
5. #include <string.h>
   
6. 
   
7. #define GR_M_PI 3.14159265358979323846            /* pi */
   
8. 
   
9. 
   
10. //complex 头文件定义了用于处理复数的功能。complex<T> 模板类型的实例表示的是复数，这里定义了 3 个特化类型：
    
11. //complex<float>、complex<double>、complex<long double>
    
12. 
    
13. typedef std::vector<void*> gr_vector_void_star;
    
14. //typedef std::vector<const void*> gr_vector_const_void_star;
    
15. typedef std::complex<float> gr_complex;
    
16. 
    
17. 
    
18. int work(int noutput_items, gr_vector_void_star& output_items)
    
19. {
    
20.     // 谁调这个函数，谁就来决定output_items的数目
    
21.     //std::cout << "python里的tb.start()来调work()吗？noutput_items=" << noutput_items << std::endl;
    
22.     //std::cout << "到底谁来反复调这个work函数开始？这个值noutput_items输出数目谁定的？=" << noutput_items << std::endl;
    
23.     //std::cout << "1个包的长度d_packet_len界面参数定的=" << d_packet_len << std::endl;
    
24.     //output_items[0]指向模块的输出缓冲区
    
25. 
    
26.     gr_complex* out = (gr_complex*)output_items[0];  // out就是一个复数数组了
    
27.     //out = new gr_complex[4096];  
    
28.     std::cout << "输出缓冲区数组长度(这里的out还是首元素地址)=" << sizeof(out) / sizeof(out[0]) << std::endl;
    
29. 
    
30.     //我加一句 需要重新编译并安装 make  sudo make install 起作用了
    
31.     //sizeof(gr_complex)=8
    
32.     //std::cout << "sizeof(gr_complex)=" << sizeof(gr_complex) << std::endl;
    
33.     //noutput_items=4096   noutput_items=3072  noutput_items=3073  noutput_items=4095
    
34.     // std::cout << "输出样点数noutput_items=" << noutput_items << std::endl;
    
35. 
    
36.     // Set output to zero 每次都清空一下
    
37.     memset(out, 0, noutput_items * sizeof(gr_complex));
    
38.    
    
39.     int d_samp_rate=32000;   //采样率  32k
    
40.     float d_amplitude=1.0;   //振幅
    
41.     std::vector<float> d_frequency;
    
42.     d_frequency.push_back(1000);  //频率1k
    
43.     int d_num_freq = d_frequency.size();// 频点个数
    
44.     std::vector<gr_complex> d_phase;   //相位
    
45.     d_phase.resize(d_num_freq);
    
46.    
    
47. 
    
48.     // Integrate phase for iq signal iq信号的集成相位
    
49.     for (int i = 0; i < noutput_items; ++i) {
    
50.         
    
51.         // Write sample  这里才是真正的写样本
    
52.         for (int k = 0; k < d_num_freq; ++k) { // Go through frequencies
    
53.             //std::cout << i << "----" << k << "----" << d_phase[k] << std::endl;
    
54.             // 遍历 所有的频率
    
55.             out[i] += d_amplitude / (float)d_num_freq * exp(d_phase[k]);
    
56.             // output amplitude is normed on d_amplitude
    
57.             d_phase[k] = gr_complex(0, std::fmod(imag(d_phase[k]) + 2 * GR_M_PI * d_frequency[k] /
    
58.                                                                (float)d_samp_rate,
    
59.                                         2 * GR_M_PI));
    
60.             //std::cout << "相位(只改变虚部)：" << d_phase[k] << std::endl;
    
61.         }
    
62.         //std::cout << "输出的值out(复数形式)=" << out[i] << std::endl;
    
63.     }
    
64. 
    
65.     // Tell runtime system how many output items we produced.
    
66.     // 告诉运行时系统我们生成了多少输出项
    
67.     return noutput_items;
    
68. }
    
69. 
    
70. 
    
71. 
    
72. int main(){
    
73.     //模拟连续波CW信号生成block的work函数
    
74.     int noutput_items = 4096;   //这个值在框架里会波动 4095 4096
    
75.     gr_complex out[noutput_items];  //复数数组 栈  int a[10];
    
76.     std::cout << "sizeof(gr_complex)=" << sizeof(gr_complex) << std::endl;
    
77. 
    
78.     //模拟一下 定时调用work
    
79.     for(int i=0;i<5;i++){
    
80.         //work(int noutput_items, gr_vector_const_void_star& input_items, gr_vector_void_star& output_items)
    
81.         //gr_vector_const_void_star input;
    
82.         gr_vector_void_star output;  
    
83.         output.push_back(out);//代表1个输出端口
    
84.         work(noutput_items, output);
    
85.         sleep(1); //停留n秒
    
86.     }
    
87. 
    
88.     return 0;
    
89. }
    

复制代码

zouqiqi

让天下人人学会人工智能！人工智能的前景一片大好！