一个简单的单边带电台Matlab仿真

一个简单的单边带电台仿真

这个仿真以真实的音频信号作为输入，设计一个单边带发信机。将基带信号调制为SSB 信号后送入带通型高斯噪声信道，加入给定功率的噪声之后，再送入单边带接收机。单边带接收机将型号解调下来，通过计算机声卡将解调信号播放出来试听效果，从而对信道信噪比与解调音质之间的关系进行主观测试。

现设计一个单边带发信机、带通信道和相应的接收机，参数定位如下值。

(1)输入信号为一个话音信号，采样率为8000Hz。话音输入后首先进行预滤波，滤波

器是一个频率范围在[300, 3400]Hz的带通滤波器，其目的是将话音频频谱限制在

3400Hz以下。单边带调制的载波频率设计为10kHz，调制输出上边带。要求观测单边带调制前后的信号功率谱。

(2)信道是一个带限高斯噪声信道，其通带频率范围是[10000, 13500]Hz。能够根据信噪

比SNR 的要求加入高斯噪声。

(3)接收机采用相干解调方式。为了模拟载波频率误差对解调话音音质的影响，设本地

载波频率为，与发信机载波平率相差200Hz。解调滤波器设计为300～3400Hz的带通滤波器。

1.对音频输入与调制的仿真

根据设置的参数，系统中信号最高频率约为14kHz。为了较好地显示调制波形，系统仿真采样率设为50kHz，满足取样定理。由于话音信号的采样率为8000Hz，与系统仿真采样率不等，因此，在进行信号处理之前，必须将话音的采样率提高到50kHz，用插值函数来做这一任务。

先编写程序将基带音频信号读入，进行[300，3400]Hz的带通滤波，并将信号采样率提高到50kHz，进行单边带调制之后，将调制输出结果保存为wav 文件，文件名为。程序如下：

% FileName:

% 功能：音频录入与调制

clc; clear all;

% 录音，采样点数为1000，采样率为8000 % jilu = wavrecord(5*8000, 8000, 'double');

% wavwrite(jilu, '');

[wav, fs] = wavread(''); % 计算声音的时间长度

t_end = 1/fs * length(wav);

% 仿真系统采样率

Fs = 50000;

% 仿真系统采样时间点

t = 1/Fs:1/Fs:t_end;

% 设计300～3400Hz 的带通滤波器H(z) [fenzi, fenmu] = butter(3, [300 3400]/(fs/2));

% 对音频信号进行滤波

wav = filter(fenzi, fenmu, wav);

% 利用插值函数将音频信号的采样率提升为Fs=50kHz

wav = interp1([1/fs:1/fs:t_end], wav, t, 'spline');

% 音频信号的希尔伯特变换

wav_hilbert = imag(hilbert(wav));

% 载波频率

fc = 10000;

% 单边带调制

SSB_OUT = wav.*cos(2*pi*fc*t) - wav_hilbert.*sin(2*pi*fc*t);

figure(1);

subplot(2, 2, 1); plot(wav(53550:53750)); axis([0 200 ]);

subplot(2, 2, 2); psd(wav, 10000, Fs); axis([0 2500 -20 10]);

subplot(2, 2, 3); plot(SSB_OUT(53550:53750)); axis([0 200 ]);

subplot(2, 2, 4); psd(SSB_OUT, 10000, Fs);

% 将SSB 调制输出存盘备用

wavwrite*SSB_OUT, '');

程序仿真结果：

2.对指定信噪比信道的仿真

仿真指定信噪比信道，仿真函数如下。

% FileName:

% 功能：计算实际信噪比clear;

[in, Fs] = wavread(''); SNRdB = 20;

out = ChanelSimulink(in, SNRdB); wavwrite(out, Fs, '');

计算信噪比为20db 时的信道输出，将结果保存为文件。输出仿真的实际测量信噪比为db。

%

% 功能：仿真指定信噪比的信道

function out = ChanelSimulink(in, SNRdb)

% SNR_db 设定信噪比

% in 输入信号序列

% out 信道输出序列

% 系统采样率

Fs = 50000;

Power_of_in = var(in);

Power_of_noise = Power_of_in/(10.^(SNRdb/10));

% 信道带宽

bandwidth = 13500 - 10000;

% 噪声功率谱密度值W/Hz

NO = Power_of_noise/bandwidth;

Gause_noise = sqrt(NO*Fs/2) .* randn(size(in));

% 噪声通道10～

[num, den] = butter(4, [10000 13500]/(Fs/2));

signal_of_filter_out = filter(num, den, in);

noise_of_filter_out = filter(num, den, Gause_noise);

SNR_dB = 10*log10(var(signal_of_filter_out)/var(noise_of_filter _out));

% 测量得出信噪比

% 信道输出

out = signal_of_filter_out + noise_of_filter_out;

SNR_dB

clear;

Fs = 50000;

% 读入信道输出信号数据

[recvsignal, Fs] = wavread('');

t = (1/Fs:1/Fs:length(recvsignal)/Fs)';

% 本地载波频率

fc_local = 10000 -200;

% 本地载波

local_carrier = cos(2*pi*fc_local.*t);

% 相干解调

xianggan_out = recvsignal.*local_carrier;

% 设计300～3400Hz 的带通滤波器[fenzi, fenmu] = butter(3, [300 3400]/(Fs/2)); demod_out = filter(fenzi, fenmu, xianggan_out); sound(demod_out/max(demod_out), Fs);

wavwrite(demod_out, Fs, '');

subplot(1, 2, 1); plot(5*demod_out(53550:53750)); axis([0 200 ]);

subplot(1, 2, 2); psd(5*demod_out, 10000, Fs); axis([0 25000 -20 10]);

解调后图像，可见仿真结果还行！

解调输出信号被保存为音频文件，并由sound 函数播放。聆听播放解调输出信号的声音可知，在20dB信道信噪比条件下，即使解调本地载波频率误差达到200Hz，声音仍然是清晰可懂的。

单边带电台频率

湛江海岸台：2182、6501、6200 舟山海翔船务公司：12408 厦门海岸台：8453、4314、12876 美军AFN夏威夷、关岛电台： 13362 & 10320（白天）、5765 & 6350（夜晚） 冰岛电台：13855 东亚太平洋航空天气台：8828 国际海事通讯波段：8195-8815、12330-13200 海岸电台 广州海岸电台 频道号 岸台发射频率(KHz) 岸台守听频率(KHz) 工作时间(北京时间) 天线方向 822 8782 8258 H24 全向 837 8716 8170 应邀而开 全向 1208 13098 12251 应邀而开 全向 1211 13107 12260 H24 全向 1225 13149 12302 H24 全向 1236 13182 12335 H24 全向 1616 17287 16405 1600-2300 全向 1653 17398 16516 H24 全向 1806 19770 18795 H24 全向 1809 19779 18804 H24 全向 1814 19794 18819 1600-2300 全向 2214 22735 22039 H24 全向 2220 22753 22057 H24 全向 每天0900-1600L每小时正点和1700-2300单小时正点在1616频道（发射频率17287 KHz）、1809频道（发射频率19779KHz）播发通话表、海上区域气象信息、新闻。 上海海岸电台 频道号 岸频率（kHz） 船频率（kHz） 工作时间（LT） 发射天线方向 405 4369 4077 1800-0600 全向 601 6501 6200 H24 全向 818 8770 8246 H24 全向 825 8791 8267 H24 全向 830 8806 8282 0600-1800 全向 1234 13176 12329 H24 全向 1238 13188 12341 H24 全向 1656 17407 16525 H24 全向 烟台海岸电台 台名 YANTAI RADIO 烟台话台 位置： 37°25′10〃N 121°30′12〃E RADIO 2182，2627 2182 J3E 全向 C409 4381 4089 J3E 全向 C816 8764 8240 J3E 全向 C1201 13077 12230 J3E 全向 C1631 17332 16450 J3E 全向 青岛海岸台 频道号 2182KHZ 2182KHZ J3E H3E 0000-2400 遇险 单边带无线电话SSB CH410 4384KHZ 4092KHZ J3E 常规 CH826 8794KHZ 8270KHZ CH1211 13107KHZ 12260KHZ 海南八所海岸电台 BASUO RADIO

基于MATLAB的FSK调制解调实现完整版

目录 一. FSK理论知识………………………………………………… 1.1FSK概念………………………………………………………………… 1.22FSK信号的波形及时间表示式………………………………………… 1.32FSK信号的产生方法…………………………………………………… 1.42FSK信号的功率谱密度………………………………………………… 1.52FSK信号的解调………………………………………………………… 1.6FSK的误码性能…………………………………………………………… 二.用MATLAB进行FSK原理及误码性能仿真……… 三、结论…………………………………………… 四、参考文献…………………………………………、 五、源程序……………………………………………

1、FSK理论知识 频率调制的最简单形式是二进制频率键控(FSK，frequency-shift keying)。FSK是调制解调器通过电话线路发送比特的方法。每个比特被转换为一个频率，0由较低的频率表示，1由较高的频率表示。 1.1、FSK概念 传“0”信号时，发送频率为f1的载波; 传“1”信号时，发送频率为f2的载波。可见，FSK是用不同频率的载波来传递数字消息的。 实现模型如下图： 1.2、2FSK信号的波形及时间表示式 根据上图模型的实现可以得到2FSK的信号波形如图：

2FSK信号的时间表达式为： 由以上表达式可见，2FSK信号由两个2ASK信号相加构成。 注意：2FSK有两种形式: （1）相位连续的2FSK； （2）相位不连续的2FSK。 在这里，我们只讨论相位不连续的频移键控信号，这样更具有普遍性。 1.3、2FSK信号的产生方法 2FSK信号的产生方法：2FSK信号可以两类方法来产生。 一是采用模拟调频的方法来产生（图1）；另一种方法是采用键控法（图2）； 图1.3-1 图1.3-2 1.4、2FSK信号的功率谱密度

中国业余电台常用频率

中国业余电台常用频率(短波) 80m波段:3.840MHz 3.843MHz 3.850MHz 3.855MHz 40m波段:7.030MHz(CW) 7.050MHz 7.O53MHz 7.055MHz 7.060MHz 7.068MHz 20m波段:14.180MHz 14.225MHz 14.270MHz 14.330MHz 15m波段:21.400MHz 10m波段:29.600MHz(FM) 14.180MHz是"B NET"工作频率.每天8:00开始.经常有老业余家出现14.330MHz是"CRSA NET"的工作频率.每周二10:00开始. 国内短波电台热闹的几个频点和守听！ 一、国内比较热闹的短波频点有： 1、7.050MHZ 模式：LSB 单边带 主要是早上和晚上很热闹，传播好的时候，比本地的中继台信号还要好，国内HAM默认的公用呼叫频率，大家常在7.050呼叫到对方后，都喜欢转到7.055、7.060、7.065、7.070、7.080做通联QSO。 2、14.270MHZ 模式：USB 单边带 主要是白天很热闹，传播好的时候，比本地的中继台信号还要好，国内HAM默认的公用呼叫频率，大家常在14.270呼叫到对方后，都喜欢转到14.270、14.275、14.280、14.265做通联QSO。 3、14.180MHZ 模式：USB 单边带 主要是早上9点前热闹，传播好的时候，比本地的中继台信号还要好，国内1级BA大哥和2级BD大哥HAM默认的公用呼叫频率，大家可以在早上收听他们的通联，可以学到更多的知识。 4、14.330MHZ 模式：USB 单边带 中国无线电运动协会台网专用频点，每周二10点有BY1PK主持，发布总部的通知和点名，大家可以收听哦。 5、21.400MHZ 模式：USB 单边带 传播特点：白天有传播，偶尔晚上的传播非常好你的日文好，级别是2级的话，可以到21.200---21.300MHZ之间，可以和日本的友台通联，日本的HAM有很多使用21MHZ的5W手持机，你随时可以呼叫日本友台，毫不费力！ 6、29.600MHZ 模式：FM 白天有传播，特别是下午2点到旁晚的6点传播，比打电话还清楚，4级火腿可以在这频段上合法使用，很有挑战性哦 二、短波传播判别心得 1、10米波29.600MHZ 模式：FM 常常是下午和旁晚才有传播，当然偶尔上午和晚上12点前的传播厉害到5W可以呼叫全球任何地区的电台，通话质量比本地台还清楚。我常常在呼叫前，收听29.600附近的鱼船电台，若能清楚收听到他们的通话时，说明这时29.600有传播了，你就可以呼叫全球电台了，成功率很高哦！哈哈，但是10米波的传播时间持续时间可能很短哦，大家要抓紧时间通联，重要的是把对方的呼号抄下。功率嘛，不需要大功率，我常常使用5W，都很好质量了。 2、20米波14.270MHZ 模式：USB 白天传播好，有无传播，可以收听广州的13.149MHZ，海事中转台，它使用率很高，能清楚收听的话，这时你的通话质量一定好，野外白天设台通联，主要就是使用这个波段和21MHZ 3、40米波7.050MHZ 模式；LSB 早晚上传播好，早上5：30-8：30前是最热闹的，通话质量很好，象听本地FM调频广播电台无线电短波电台是呼叫远方电台的主要工具，常常用来呼叫全球电台，晚上传播好的时候，频率上是密密麻麻的信号，只要你的级别允许你使用的频率范围，都可以用英文呼叫国外电台，是无线电通联技术训练和英语学习的好伙伴。晚上主要在14.150---14.250MHZ之间，很多国外电台都在此呼叫，是练习英语口语和听力的好机会，凌晨主要在7MHZ 上，很是热闹的。 晚上我常在14.195、14.200、14.202、14.210等频点收听国外电台每天晚上很多不同国家的无线电爱好者活跃在这里也听到了很有趣的不同国家版本的英语口语！ 无线电短波通联我为它着迷、为他发狂！有趣、刺激、是激发自己不断创新、提高的源头每每

基于MATLAB SIMULINK的FM调制解调

摘要 在模拟通信系统中，由模拟信源产生的携带有信息的消息经过传感器转换成电信号。模拟基带信号在经过调制将低通频谱搬移到载波频率上适应信道，最终解调还原成电信号。本文应用了频率调制法产生调制解调信号。本论文中主要通过对SIMULINK工具箱的学习和使用，利用其丰富的模板以及本科对通信原理知识的掌握，完成了FM信号的调制与解调，以及用SIMULINK进行设计和仿真。首先利用简单的正玄波信号发生器作为信源，对模拟信号进行FM调制解调原理的仿真。 关键词：调制解调；FM ；MATLAB；SIMULINK仿真

Abstract In the simulation of communication systems, generated by the analog source carrying a message through the sensor into electrical signals. Analog baseband signal after the modul- -ation of the low pass spectrum to carrier frequency to adapt to the channel, the final reducti- -on into electrical signal demodulation. This paper applied the frequency modulation method to generate the signal modulation and demodulation. Mainly through the study and use of SIMULINK toolbox in this thesis, with its rich template and undergraduate course on comm--unication theory knowledge,the modulation and demodulation of FM signal, as well as the design and simulation with SIMULINK. Firstly, sine wave signal generator is simple as the source, simulation FM modulation anddemodulation principle of analogue signals. Then, using the song as the source. Keywords: modulation and demodulation；FM; MATLAB; SIMULINK simulation

matlab矩阵操作汇总

matlab矩阵操作大全 1.1.1数值矩阵的生成 1.实数值矩阵输入 MATLAB的强大功能之一体现在能直接处理向量或矩阵。当然首要任务是输入待处理的向量或矩阵。 不管是任何矩阵（向量），我们可以直接按行方式输入每个元素：同一行中的元素用逗号（，）或者用空格符来分隔，且空格个数不限；不同的行用分号（；）分隔。所有元素处于一方括号（［］）内；当矩阵是多维（三维以上），且方括号内的元素是维数较低的矩阵时，会有多 重的方括号。如： >> Time = [11 12 1 2 3 4 5 6 7 8 9 10] Time = 11 12 1 2 3 4 5 6 7 8 9 10 >> X_Data = [2.32 3.43 ； 4.37 5.98] X_Data = 2.43 3.43 4.37 5.98 >> vect_a = [1 2 3 4 5] vect_a = 12~ 3 4 5 >> Matrix_B = [1 2 3 ； >> 2 3 4 ； 3 4 5] Matrix_B = 1 2 3 2 3 4 3 4 5 >> Null_M = [ ] %生成一个空矩阵

2 .复数矩阵输入 复数矩阵有两种生成方式： 第一种方式 1/ 13

例1-1 >> a="2".7。b=13/25 >> C 二[1,2*a+i*b,b*sqrt(a) 5.4000 + 0.5200i 0.8544 5.3000 4.5000 第2种方式 例1-2 >> R=[1 2 3 。 4 5 6], M=[11 12 13 。 14 15 16] R = I 2 3 4 5 6 M = II 12 13 14 15 16 >> CN="R"+i*M CN = 1.0000 +11.0000i 2.0000 +12.0000i 3.0000 +13.0000i 4.0000 +14.0000i 5.0000 +15.0000i 6.0000 +16.0000i si n(pi/4),a+5*b,3.5+1] C= 1.0000 0.7071

TCR-230 单边带短波通信电台试验器使用说明

机载电子设备维修技术人员培训系列教材 机载电子设备试验器使用操作说明 TCR-230单边带短波电台试验器 作 者： 贺 军 笔 名： STDHJ ( 作者：贺军 中信海直通用航空维修工程有限公司 电子工程师 ) TCR-230单边带短波通讯电台校验技术条件和校验规程 一、TCR-230单边带短波电台校验技术条件： 1.1、概述： HF-230是一部轻型机载远程通信的单边带电台，电台工作频率为2.0000MHZ到29.9999MHZ范围内，道波道频率间隔100Hz、预编程波道16个。可提供100W的峰包功率，在电台频率范围内，可直接调谐频率使用，也可以预先编程16个波道；工作方式有USB（上边带）、AM（兼容调幅），为适应沿海公用通信和海陆通信的需要，电台还可以工作在TEL PLT CAR或TEL SUP CAR （电话方式、前者为导频的上边带，后者为抑制载波的上边带）；电台广泛地使用了数字信号处理技术和集成电路，频率合成器为收发机提供了高稳定和高准确度的频率源；电台的器件、电路板的结构组件设计为便于调整和拆卸、修理的模块化结构。 1.2、适用范围： HF-230单边带电台由控制盒、收发机部件、功率放大器部件、自动天线调谐器及电缆、馈线、天线组成，本技术条件适用于TCR-230单边带收发机在试验室进行检验和修理的技术条件依据。 1.3、引用文件： Collins HF-230 High-Frequency Communications System 第一册 TCR-230收发机部分（1984年2月修改版） 版本：523-0772873-00211A 23-19-10 1.4、高频电台技术指标: 1.4.1、电源电压：直流22V到30.5V；额定27.5V±2%DC。 1.4.2、接收电流：在直流27.5V时1.9A。 1.4.3、发射电流：在直流27.5V时 2.2A。 1.4.4、频率范围： 2.000MHZ-29.9999MHZ。 1.4.5、频率稳定性：±20HZ。

基于MATLAB的FSK调制解调1

基于MATLAB的FSK调制解调 学生姓名：段斐指导老师：吴志敏 摘要本课程设计利用MATLAB集成环境下的M文件，编写程序来实现FSK 的调制解调，并绘制出解调前后的时域和频域波形及叠加噪声时解调前后的时频波形，并观察解调前后频谱有何变化以加深对F SK信号解调原理的理解。对信号叠加噪声，并迚行解调，绘制出解调前后信号的时频波形，改变噪声功率迚行解调，根据运行结果和波形来分析该解调过程的正确性及信道对信号传输的影响。完成整个FSK的调制解调过程。程序开发平台为MATLAB7.1，使用其自带的M文件实现。运行平台为Windows 2000。 关键词：程序设计；FSK ；调制解调；MATLAB7.1；M文件 1引言 本课程设计是利用MATLAB集成环境下的M文件，编写程序来实现FSK 的调制解调，并绘制出解调前后的时域和频域波形及叠加噪声时解调前后的时频波形，根据运行结果和波形来分析该解调过程的正确性及信道对信号传输的影响。 1.1课程设计目的 此次课程设计的目的是熟悉MATLAB中M文件的使用方法，编写M文件实现FSK的调制和解调，绘制出FSK信号解调前后在时域和频域中的波形，观察调解前后频谱的变化，再对信号迚行噪声叠加后解调同样绘制解调前后的

信号时频波形，最后改变噪声功率迚行调解，分析噪声对信号传输造成的影响，加深对FSK信号解调原理的理解。 1.2课程设计要求 熟悉MATLAB中M文件的使用方法，并在掌握FSK调制解调原理的基础上，编写出F SK调制解调程序。在M文件环境下运行程序绘制出F SK信号解调前后在时域和频域中的波形，观察波形在解调前后的变化，对其作出解释，同时对信号加入噪声后解调，得到解调后的时频波形，分析噪声对信号传输造成的影响。解释所得到的结果。 1.3课程设计步骤 本课程设计采用M文件编写的方法实现二迚制的FSK的调制与解调，然后在信号中叠加高斯白噪声。一，调用dmode函数实现FSK的解调，并绘制出F SK信号调制前后在时域和频域中的波形，两者比较。二，调用ddemod函数解调，绘制出F SK信号解调前后在时域和频域中的波形，两者比较。三，调用awgn函数在新海中叠加不同信噪比的噪声，绘制在各种噪声下的时域频域图。最后分析结果。 1.4设计平台简介 Matlab是美国MathWorks公司开发的用于概念设计,算法开发,建模仿真,实时实现的理想的集成环境。是目前最好的科学计算类软件。 作为和Mathematica、Maple并列的三大数学软件。其强项就是其强大的矩阵计算以及仿真能力。Matlab的由来就是Matrix + Laboratory = Matlab，这个软件在国内也被称作《矩阵实验室》。Matlab提供了自己的编译器：全面兼容C++以及Fortran两大语言。Matlab 7.1于2005.9最新发布-完整版,提供了

矩阵在matlab中的基本命令

一、矩阵的表示 在MATLAB中创建矩阵有以下规则： a、矩阵元素必须在”[ ]”内； b、矩阵的同行元素之间用空格（或”,”）隔开； c、矩阵的行与行之间用”;”（或回车符）隔开； d、矩阵的元素可以是数值、变量、表达式或函数； e、矩阵的尺寸不必预先定义。 二，矩阵的创建： 1、直接输入法 最简单的建立矩阵的方法是从键盘直接输入矩阵的元素，输入的方法按照上面的规则。建立向量的时候可以利用冒号表达式，冒号表达式可以产生一个行向量，一般格式是：e1:e2:e3，其中e1为初始值，e2为步长，e3为终止值。还可以用linspace函数产生行向量，其调用格式为：linspace(a,b,n) ，其中a和b是生成向量的第一个和最后一个元素，n是元素总数。 2、利用MATLAB函数创建矩阵 基本矩阵函数如下： (1) ones()函数：产生全为1的矩阵，ones(n)：产生n*n维的全1矩阵，ones(m,n)：产生m*n 维的全1矩阵； (2) zeros()函数：产生全为0的矩阵；

(3) rand()函数：产生在（0，1）区间均匀分布的随机阵； (4) eye()函数：产生单位阵； (5) randn()函数：产生均值为0，方差为1的标准正态分布随机矩阵。 3、利用文件建立矩阵 当矩阵尺寸较大或为经常使用的数据矩阵，则可以将此矩阵保存为文件，在需要时直接将文件利用load命令调入工作环境中使用即可。同时可以利用命令reshape对调入的矩阵进行重排。reshape(A,m,n)，它在矩阵总元素保持不变的前提下，将矩阵A重新排成m*n的二维矩阵。 二、矩阵的简单操作 1．获取矩阵元素 可以通过下标（行列索引）引用矩阵的元素，如Matrix(m,n)。 也可以采用矩阵元素的序号来引用矩阵元素。 矩阵元素的序号就是相应元素在内存中的排列顺序。 在MATLAB中，矩阵元素按列存储。 序号(Index)与下标(Subscript )是一一对应的，以m*n矩阵A为例，矩阵元素A(i,j)的序号为(j-1)*m+i。 其相互转换关系也可利用sub2ind和ind2sub函数求得。 2．矩阵拆分

全国电台频率表

全国电台频率表

0.5-1.8M 广播频段 MF(AM) 525KHz-1605KHz 9KHz 中波调幅广播 北京电台首都生活广播603 中央电台中国之声639 中央电台经济之声720 中央电台文艺之声(第9套) 747 北京电台新闻频道828/846 国际电台HIT-FM 900 北京电台体育频道927 国际电台ROUND THE CLOCK 846/1008 中央电台朝鲜语广播1206 北京电台经济频道1026 国际电台英语频道1251 1.8-30M 短波波段 2182kHz 为无线电话的国际遇险频率 4125 kHz频率用作2182 kHz频率的补充 HF(SW) 3.5-29.7MHz 9KHz 短波调幅广播及单边带通讯 玩具遥控器 通信设备带宽(KHz)发射功率(mW)频率(MHz) 遥控设备带宽(KHz) 发射功率(W) 频率 (MHz) 1 <1 2 ≤100 26.965 1 ≤8 ≤1 26.975 2 <12 ≤100 26.985 2 ≤8 ≤1 26.995 3 <12 ≤100 27.005 3 ≤8 ≤1 27.015 4 <12 ≤100 27.02 5 4 ≤8 ≤1 27.045 5 <12 ≤100 27.055 5 ≤8 ≤1 27.065 6 <12 ≤100 27.075 6 ≤8 ≤1 27.095 7 <12 ≤100 27.105 7 ≤8 ≤1 27.115 8 <12 ≤100 27.125 8 ≤8 ≤1 27.145 9 <12 ≤100 27.165 9 ≤8 ≤1 27.195 10 <12 ≤100 27.185 10 ≤8 ≤1 27.255 30-76（59）M 50M业余波段 无绳电话使用频率划分表一 座机频率(MHz) 手机发射频率(MHz)

基于MATLAB的FSK调制解调 (1)

基于MATLAB的FSK的实验报告 姓 1.1

实现对FSK的MATLAB仿真. 重点研究问题: (1) 对FSK的概念、组成以及性能分析方法有深入的研究； (2) FSK调制与解调的原理及应用MATLAB软件实现仿真的方案. 1.2 FSK信号的调制方法 移频键控(FSK)：用数字调制信号的正负控制载波的频率。当数字信号的振幅为正时载波频率为f1，当数字信号的振幅为负时载波频率为 f2。有时也把代表两个以上符号的多进制频率调制称为移频键控。移频键控能区分通路，但抗干扰能力不如移相键控和差分移相键控。他的主要调制方法有以下两种: 方法一： 用一个矩形脉冲序列对一个载波进行调频。 图2-3 2FSK信号的产生(一) 方法二:键控法 图2-4 2FSK信号的产生(二) 键控法是利用矩形脉冲()t b来控制开关电路对两个不同的独立频率源进行选通。

1.3 FSK解调的方法 常见的FSK解调方法有两种:相干解调法与非相干解调法.现在我将对这两种解法。 1.4 设计总思路 如下图所示,我将FSK的调制与FSK的解调独立开作为两个子函数,其中FSK调制的输出即可作为FSK解调的输入信号.最后设计一主函数main将两个子函数同时调用完成整个仿真过程。 图3-1 设计总思路图 2.1 FSK调制的仿真设计 本文主要是对2FSK进行调制,而2FSK可看做是基带信号与载波频率的结合就可.FSK的产生思路参考的是键控法,如图4

图3-2 2FSK信号的产生(二) 2.2 FSK解调的仿真设计 如上图所示的FSK信号的相干检测原理图,FSK信号可以采用两个乘法检测器进行相干检测. 上图中输入信号为2FSK信号加上噪声组成 带通滤波器2的设计类似滤波器1,只是更改频率为fc2就可.

Matlab界面、矩阵操作

科学计算与MATLAB语言

教学内容 第1讲介绍MATLAB概述与运算基础第2讲介绍MATLAB程序设计 第3讲MATLAB文件操作 第4讲绘图形功能 第5讲线形代数中的数值计算问题 第6讲数据处理方法与多项式 第7讲MATLAB符号计算 第8讲MATLAB的图形用户界面设计

第一讲MATLAB概述与运算基础MATLAB语言是一种广泛应用于工程计算及数值分析领域的新型高级语言，自1984年由美国MathWorks公司推向市场以来，历经十多年的发展与竞争，现已成为国际公认的最优秀的工程应用开发环境。MATLAB功能强大、简单易学、编程效率高，深受广大科技工作者的欢迎。

在欧美各高等院校，MATLAB已经成为线性代数、自动控制理论、数字信号处理、时间序列分析、动态系统仿真、图像处理等课程的基本教学工具，成为大学生、硕士生以及博士生必须掌握的基本技能。 MATLAB: 是英文MATrix LABorotory(矩阵实验室) 的缩写。

一. MATLAB特点:

1. 数值计算和符号计算功能 MATLAB的数值计算功能包括：矩阵运算、多项式和有 理分式运算、数据统计分析、数值积分、优化处理等。符 号计算将得到问题的解析解。 2.MATLAB语言 MATLAB除了命令行的交互式操作以外，还可以程序方式 工作。使用MATLAB可以很容易地实现C或FORTRAN语言的几乎全部功能，包括Windows图形用户界面的设计。 3.图形功能 MATLAB提供了两个层次的图形命令：一种是对图形句柄进行的低级图形命令，另一种是建立在低级图形命令之上的高级图形命令。利用MATLAB的高级图形命令可以轻而易举地绘制二维、三维乃至四维图形，并可进行图形和坐标的标识、视角和光照设计、色彩精细控制等等。

电台频率表

8967kHz：全军聊天频率。0700至2200 日本气象：13597、9970 台湾气象：13900 日本报纸：12745 东南亚语种、英语海岸电台强信号经常出没波段：16590-16860（不详）！未知电台（13288）专门播送怪语言 天津海岸台：13092 上海海岸台：8773 烟台海岸台：13077、17332、8764、4381 湛江海岸台：2182、6501、6200 舟山海翔船务公司：12408 厦门海岸台：8453、4314、12876 美军AFN夏威夷、关岛电台： 13362 & 10320（白天）、5765 & 6350（夜晚） 冰岛电台：13855 东亚太平洋航空天气台：8828 国际海事通讯波段：8195-8815、12330-13200 海岸电台 广州海岸电台 频道号岸台发射频率(KHz) 岸台守听频率(KHz) 工作时间(北京时间) 天线方向 822 8782 8258 H24 全向 837 8716 8170 应邀而开全向 1208 13098 12251 应邀而开全向 1211 13107 12260 H24 全向 1225 13149 12302 H24 全向 1236 13182 12335 H24 全向 1616 17287 16405 1600-2300 全向 1653 17398 16516 H24 全向 1806 19770 18795 H24 全向 1809 19779 18804 H24 全向 1814 19794 18819 1600-2300 全向 2214 22735 22039 H24 全向 第1 页单边带SSB电台频率表(参考表).txt 2220 22753 22057 H24 全向

基于MATLAB的ASK调制解调实现

基于MATLAB的ＡSK调制解调实现

————————————————————————————————作者：————————————————————————————————日期： ?

长沙理工大学 《通信原理》课程设计报告 学院专业 班级学号 学生姓名指导教师 课程成绩完成日期2016年1月8日

课程设计成绩评定 学院专业 班级学号 学生姓名指导教师 课程成绩完成日期20１6年1月8日 指导教师对学生在课程设计中的评价 评分项目优良中及格不及格课程设计中的创造性成果 学生掌握课程内容的程度 课程设计完成情况 课程设计动手能力 文字表达 学习态度 规范要求 课程设计论文的质量 指导教师对课程设计的评定意见 综合成绩指导教师签字２01６年1月8日

课程设计任务书 城南学院通信工程专业 课程名称通信原理课程设计时间２０15/2016学年第一学期１7~19 周 学生姓名指导老师 题目基于MＡTLAＢ的ＡＳK调制解调实现 主要内容: 利用MATＬAＢ集成环境下的M文件，编写程序来实现AＳＫ的调制解调， 要求采样频率为３６0ＨZ,并绘制出解调前后的时域和频域波形及叠加噪声时解 调前后的时频波形,根据运行结果和波形来分析该解调过程的正确性及信道对信 号传输的影响。 要求： 1)熟悉MATLAB中M文件的使用方法,并在掌握ASK调制解调原理 的基础上，编写出ASＫ调制解调程序。 ２)绘制出ASK信号解调前后在时域和频域中的波形,并观察解调前后频谱有何变化以加深对ＡSK信号解调原理的理解。 3）对信号叠加噪声，并进行解调，绘制出解调前后信号的时频波形,改变噪声功率进行解调,分析噪声对信号传输造成的影响。 4)在老师的指导下，要求独立完成课程设计的全部内容,并按要求编写课 程设计学年论文,能正确阐述和分析设计和实验结果。 应当提交的文件: （１）课程设计学年论文。 (2)课程设计附件。

matlab中的矩阵的基本运算命令

1.1 矩阵的表示 1.2 矩阵运算 1.2.14 特殊运算 1．矩阵对角线元素的抽取 函数diag 格式X = diag(v,k) %以向量v的元素作为矩阵X的第k条对角线元素，当k=0时，v为X的主对角线；当k>0时，v为上方第k条对角线；当k<0时，v为下方第k条对角线。 X = diag(v) %以v为主对角线元素，其余元素为0构成X。 v = diag(X,k) %抽取X的第k条对角线元素构成向量v。k=0：抽取主对角线元素；k>0：抽取上方第k条对角线元素；k<0抽取下方第k条对角线元素。 v = diag(X) %抽取主对角线元素构成向量v。 2．上三角阵和下三角阵的抽取 函数tril %取下三角部分 格式L = tril(X) %抽取X的主对角线的下三角部分构成矩阵L L = tril(X,k) %抽取X的第k条对角线的下三角部分；k=0为主对角线；k>0为主对角线以上；k<0为主对角线以下。函数triu %取上三角部分 格式U = triu(X) %抽取X的主对角线的上三角部分构成矩阵U U = triu(X,k) %抽取X的第k条对角线的上三角部分；k=0为主对角线；k>0为主对角线以上；k<0为主对角线以下。3．矩阵的变维 矩阵的变维有两种方法，即用“：”和函数“reshape”，前者主要针对2个已知维数矩阵之间的变维操作；而后者是对于一个矩阵的操作。 （1）“：”变维 （2）Reshape函数变维 格式 B = reshape(A,m,n) %返回以矩阵A的元素构成的m×n矩阵B B = reshape(A,m,n,p,…) %将矩阵A变维为m×n×p×… B = reshape(A,[m n p…]) %同上 B = reshape(A,siz) %由siz决定变维的大小，元素个数与A中元素个数 相同。 （5）复制和平铺矩阵 函数repmat 格式 B = repmat(A,m,n) %将矩阵A复制m×n块，即B由m×n块A平铺而成。 B = repmat(A,[m n]) %与上面一致 B = repmat(A,[m n p…]) %B由m×n×p×…个A块平铺而成 repmat(A,m,n) %当A是一个数a时，该命令产生一个全由a组成的m×n矩阵。 1.3 矩阵分解 1.3.1 Cholesky分解 函数chol 格式R = chol(X) %如果X为n阶对称正定矩阵，则存在一个实的非奇异上三角阵R，满足R'*R = X；若X非正定，则产生错误信息。 [R,p] = chol(X) %不产生任何错误信息，若X为正定阵，则p=0，R与上相同；若X非正定，则p为正整数，R是有序的上三角阵。 1.3.2 LU分解

AM调制解调及matlab仿真程序和图

（1）所用滤波器函数：巴特沃斯滤波器 % 注: wp（或Wp）为通带截止频率 ws(或Ws)为阻带截止频率 Rp为通带衰减 As为阻带衰减 %butterworth低通滤波器原型设计函数要求Ws>Wp>0 As>Rp>0 function [b,a]=afd_butt(Wp,Ws,Rp,As) N=ceil((log10((10^(Rp/10)-1)/(10^(As/10)-1)))/(2*log10(Wp/Ws))); %上条语句为求滤波器阶数 N为整数 %ceil 朝正无穷大方向取整 fprintf('\n Butterworth Filter Order=%2.0f\n',N) OmegaC=Wp/((10^(Rp/10)-1)^(1/(2*N))) %求对应于N的3db截止频率 [b,a]=u_buttap(N,OmegaC); （2）傅里叶变换函数 function [Xk]=dft(xn,N) n=[0:1:N-1]; k=[0:1:N-1]; WN=exp(-j*2*pi/N); nk=n'*k; WNnk=WN.^(nk); Xk=xn*WNnk; 设计部分： 1.普通AM调制与解调 %单音普通调幅波调制y=amod(x,t,fs,t0,fc,Vm0,ma)要求fs>2fc %x调制信号，t调制信号自变量,t0采样区间，fs采样频率， %fc载波频率，Vm0输出载波电压振幅，ma调幅度 t0=0.1;fs=12000; fc=1000;Vm0=2.5;ma=0.25; n=-t0/2:1/fs:t0/2; x=4*cos(150*pi*n); %调制信号 y2=Vm0*cos(2*pi*fc*n); %载波信号figure(1) subplot(2,1,1);plot(n,y2); axis([-0.01,0.01,-5,5]); title('载波信号'); N=length(x); Y2=fft(y2); subplot(2,1,2); plot(n,Y2); title('载波信号频谱'); %画出频谱波形y=Vm0*(1+ma*x/Vm0).*cos(2*pi*fc*n); figure(2) subplot(2,1,1);plot(n,x) title('调制信号'); subplot(2,1,2) plot(n,y) title('已调波信号'); X=fft(x);Y=fft(y);

电机与拖动基础及MATLAB仿真习题答案(第四章)

4-14 一台直流电动机技术数据如下：额定功率PN=40kW ，额定电压UN=220V ，额定转速nN=1500r/min ，额定效率η=87.5%，求电动机的额定电流和额定负载时的输入功率？ 解： （1）额定电流 （2）输入功率 4-15 一台直流发电机技术数据如下：额定功率PN=82kW ，额定电压UN=230V ，额定转速nN=970r/min ，额定效率η=90%，求发电机的额定电流和额定负载时的输入功率？ 解： （1）额定电流 （2）输入功率 4-16 已知一台直流电机极对数p=2，槽数Z 和换向片数K 均等于22，采用单叠绕组。试求：（1）绕组各节距；（2）并联支路数。 解：（1）第一节距54 24222y 1=-=±=εp z ，为短距绕组。 单叠绕组的合成节距及换向器节距均为1，即1y ==k y 第二节距415y 12=-=-=y y （2）并联支路数等于磁极数，为4。 4-17 已知直流电机极数2p=6，电枢绕组总导体数N=400，电枢电流Ia=10A ，气隙每极磁通Φ=2.1×10-2Wb ，试求：（1）采用单叠绕组时电枢所受电磁转矩；（2）绕组改为单波保持支路电流ia 不变时的电磁转矩。 解: 电枢绕组为单叠绕组时,并联支路对数a=p=3, 电磁转矩 m N I a pN T a ?=?????=Φ=38.1310021.03 14.3240032π 如果把电枢绕组改为单波绕组, 保持支路电流a i 的数值不变,则电磁转矩也不变,仍为 13.369m N ?,因为无论是叠绕组还是波绕组,所有导体产生的电磁转矩的方向是一致 的, 保持支路电流a i 不变,就保持了导体电流不变,也就保持了电磁转矩不变。 也可以用计算的方法: 单叠绕组时并联支路数为6, 支路电流 A I i a a 6106== 改为单波绕组, 保持支路电流a i 的数值不变,仍为 A 610,而并联支路数为2 (a=1), 电枢电流A i I a a 3 102== 电磁转矩 m N T ?=?????=38.133 10021.0114.324003。 A U P I N N N N 79.207875.022010403 =??==ηkW I U P N N 71.4579.2072201=?=?=A U P I N N N 5.35623010823 =?==KW P P N 11.911==η

完整word版,msk的调制解调MATLAB源代码

msk的调制解调MATLAB源代码 function out = delay(data,n,sample_number) %data:延迟的数据 %n:延迟码元个数 %sample_number:码元采样个数 out = zeros(1,length(data)); out(n*sample_number+1:length(data)) = data(1:length(data)-n*sample_number); function [data_diff] = difference(data) %差分编码 %************************************************************************* * %data 输入信号 %data_diff 差分编码后信号 %************************************************************************* *

%-------------------------------------------------------------------------- data_diff = zeros(1,length(data)); data_diff(1) = 1 * data(1); %1为差分编码的初始参考值 for i = 2:length(data) data_diff(i) = data_diff(i-1) * data(i); end %************************************************************************* * function [signal_out,I_out,Q_out] = mod_msk(data,data_len,sample_number,Rb) %MSK基带调制 %************************************************************************* * % data 调制信号 % data_len 码元个数 % sample_number 每个码元采样点数

16QAM调制解调(MATLAB)

题目： 基于MATLAB 的16QAM 及32QAM 系统的仿真 原理： QAM 是一种矢量调制，将输入比特映射到一个复平面，形成复数调制信号，然后将I 信号和Q 信号（实部虚部）分量采用幅度调制，分 别对应调制在相互正交的两个载波（cos t ω，sin t ω）上。下图为MQAM 的调制原理图。 MQAM 的信号表达式： ()()( )cos sin 1,2,...,, 0C S C S i i T C i T C S i i s t a g t t a g t t i M t T a a ωω=-=≤≤与 上述表达式可以看出，QAM 为两个正交载波振幅相位调制的结合。波形矢量可以表示为： ()()()11221,2,...,, 0i i i S s t s f t s f t i M t T =+=≤≤

( )()( )()()()()()12110 220 cos ,0sin ,01,2,...,1,2,...,S S T C S T C S T i i T i i f t t t t T f t t t t T s s t f t dt i M s s t f t dt i M ωω=≤≤= ≤≤====?? MQAM 信号最佳接收： 实验仿真条件： 码元数量设定为10000个，基带信号频率1HZ ，抽样频率32HZ ，载波频率4HZ 。 实验结果分析：

对于QAM ，可以看成是由两个相互正交且独立的多电平ASK 信号叠加而成。因此，利用多电平误码率的分析方法，可得到M 进制QAM 的误码率为： ])(1l o g 3[)1 1(0 22n E L L e r f c L P b e -- = 式中，M L =，Eb 为每码元能量，n 0为噪声单边功率谱密度。 通过调整高斯白噪声信道的信噪比SNR （Eb/No ），可以得到如图所示的误码率图： -1-0.50 0.51 1.52 2.5 10 -3 10 -2 10 -1 10 QAM 信号误码率分析 信噪比 误码率

matlab矩阵的表示和简单操作

matlab矩阵的表示和简单操作 一、矩阵的表示 在MATLAB中创建矩阵有以下规则： a、矩阵元素必须在”[ ]”内； b、矩阵的同行元素之间用空格（或”,”）隔开； c、矩阵的行与行之间用”;”（或回车符）隔开； d、矩阵的元素可以是数值、变量、表达式或函数； e、矩阵的尺寸不必预先定义。 二，矩阵的创建： 1、直接输入法 最简单的建立矩阵的方法是从键盘直接输入矩阵的元素，输入的方法按照上面的规则。建立向量的时候可以利用冒号表达式，冒号表达式可以产生一个行向量，一般格式是：e1:e2:e3，其中e1为初始值，e2为步长，e3为终止值。还可以用linspace函数产生行向量，其调用格式为：linspace(a,b,n) ，其中a和b是生成向量的第一个和最后一个元素，n是元素总数。 2、利用MATLAB函数创建矩阵 基本矩阵函数如下： (1) ones()函数：产生全为1的矩阵，ones(n)：产生n*n维的全1矩阵，ones(m,n)：产生m*n维的全1矩阵； (2) zeros()函数：产生全为0的矩阵； (3) rand()函数：产生在（0，1）区间均匀分布的随机阵； (4) eye()函数：产生单位阵； (5) randn()函数：产生均值为0，方差为1的标准正态分布随机矩阵。 3、利用文件建立矩阵 当矩阵尺寸较大或为经常使用的数据矩阵，则可以将此矩阵保存为文件，在需要时直接将文件利用load命令调入工作环境中使用即可。同时可以利用命令reshape对调入的矩阵进行重排。reshape(A,m,n)，它在矩阵总元素保持不变的前提下，将矩阵A重新排成m*n 的二维矩阵。 二、矩阵的简单操作 1．获取矩阵元素 可以通过下标（行列索引）引用矩阵的元素，如Matrix(m,n)。 也可以采用矩阵元素的序号来引用矩阵元素。 矩阵元素的序号就是相应元素在内存中的排列顺序。 在MATLAB中，矩阵元素按列存储。 序号(Index)与下标(Subscript )是一一对应的，以m*n矩阵A为例，矩阵元素A(i,j)的序号为(j-1)*m+i。 其相互转换关系也可利用sub2ind和ind2sub函数求得。