使用MATLAB将波形数据存入TXT

cyy13623649348

这里使用MATLAB产生两个不同频率的正弦信号s1与s2，将两个信号相加得到s，即s=s1+s2，最后将信号s的数据存入TXT文件。

1. f1=1000;     %信号1频率为1kHz
   
2. f2=3000;    %信号2频率为3kHz
   
3. Fs=20000;  %采样频率为20KHz
   
4. N=12;         %量化位数
   
5. %产生信号
   
6. t=0:1/Fs:0.1;
   
7. c1=2*pi*f1*t;
   
8. c2=2*pi*f2*t;
   
9. s1=sin(c1);%产生正弦波
   
10. s2=sin(c2);%产生正弦波
    
11. s=s1+s2;   %产生两个单载波合成后的信号
    
12. %产生随机序列信号
    
13. 
    
14. %归一化处理
    
15. s=s/max(abs(s));
    
16. 
    
17. %12比特量化
    
18. Q_s=round(s*(2^(N-1)-1));
    
19. 
    
20. %求信号的幅频响应
    
21. m_s=20*log(abs(fft(Q_s,1024)))/log(10); m_s=m_s-max(m_s);
    
22. 
    
23. %设置幅频响应的横坐标单位为Hz
    
24. x_f=[0:(Fs/length(m_s)):Fs/2];
    
25. %只显示正频率部分的幅频响应
    
26. mf_s=m_s(1:length(x_f));
    
27. 
    
28. %绘制时域图形
    
29. figure(1)
    
30. plot(t,s,'-.');
    
31. set(gca,'XLim',[0 0.01]);%x轴的数据显示范围
    
32. xlabel('时间');ylabel('归一化幅值');title('Matlab仿真合成信号时域图');
    
33. grid;
    
34. %绘制幅频响应曲线
    
35. figure(2)
    
36. plot(x_f,mf_s,'-.');
    
37. xlabel('频率(Hz)');ylabel('幅度(dB)');title('Matlab仿真合成单频信号的频谱');
    
38. legend('输入信号频谱');
    
39. grid;
    
40. 
    
41. %将生成的数据以十进制数据格式写入txt文件中
    
42. fid=fopen('J:\E4_7_Int_s.txt','w');   %自己改一下存放目录，这里我存在J盘
    
43. fprintf(fid,'%8d\r\n',Q_s);
    
44. fprintf(fid,';');
    
45. fclose(fid);
    
46. 
    
47. %将生成的数据以二进制数据格式写入txt文件中
    
48. fid=fopen('J:\E4_7_Bin_s.txt','w');
    
49. for i=1:length(Q_s)
    
50.     B_s=dec2bin(Q_s(i)+(Q_s(i)<0)*2^N,N)
    
51.     for j=1:N
    
52.        if B_s(j)=='1'
    
53.            tb=1;
    
54.        else
    
55.            tb=0;
    
56.        end
    
57.        fprintf(fid,'%d',tb);  
    
58.     end
    
59.     fprintf(fid,'\r\n');
    
60. end
    
61. fprintf(fid,';');
    
62. fclose(fid);

复制代码

   下面是运行后产生的波形和频谱：