【数字信号调制】基于matlab GUI AM+DSB+SSB+ASK+FSK+PSK调制解调【含源码 058期】.zip

function varargout = communication(varargin) % COMMUNICATION MATLAB code for communication.fig % COMMUNICATION, by itself, creates a new COMMUNICATION or raises the existing % singleton*. % % H = COMMUNICATION returns the handle to a new COMMUNICATION or the handle to % the existing singleton*. % % COMMUNICATION('CALLBACK',hObject,eventData,handles,...) calls the local % function named CALLBACK in COMMUNICATION.M with the given input arguments. % % COMMUNICATION('Property','Value',...) creates a new COMMUNICATION or raises the % existing singleton*. Starting from the left, property value pairs are % applied to the GUI before communication_OpeningFcn gets called. An % unrecognized property name or invalid value makes property application % stop. All inputs are passed to communication_OpeningFcn via varargin. % % *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one % instance to run (singleton)". % % See also: GUIDE, GUIDATA, GUIHANDLES % Edit the above text to modify the response to help communication % Last Modified by GUIDE v2.5 02-Sep-2019 08:36:43 % Begin initialization code - DO NOT EDIT gui_Singleton = 1; gui_State = struct('gui_Name', mfilename, ... 'gui_Singleton', gui_Singleton, ... 'gui_OpeningFcn', @communication_OpeningFcn, ... 'gui_OutputFcn', @communication_OutputFcn, ... 'gui_LayoutFcn', [] , ... 'gui_Callback', []); if nargin && ischar(varargin{1}) gui_State.gui_Callback = str2func(varargin{1}); end if nargout [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:}); else gui_mainfcn(gui_State, varargin{:}); end % End initialization code - DO NOT EDIT % --- Executes just before communication is made visible. function communication_OpeningFcn(hObject, eventdata, handles, varargin) % This function has no output args, see OutputFcn. % hObject handle to figure % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % varargin command line arguments to communication (see VARARGIN) % Choose default command line output for communication handles.output = hObject; % Update handles structure guidata(hObject, handles); % UIWAIT makes communication wait for user response (see UIRESUME) % uiwait(handles.figure1); % --- Outputs from this function are returned to the command line. function varargout = communication_OutputFcn(hObject, eventdata, handles) % varargout cell array for returning output args (see VARARGOUT); % hObject handle to figure % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Get default command line output from handles structure varargout{1} = handles.output; % -------------------------------------------------------------------- function menu_1_Callback(hObject, eventdata, handles) % hObject handle to menu_1 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % -------------------------------------------------------------------- function Untitled_2_Callback(hObject, eventdata, handles) % hObject handle to Untitled_2 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) %=========程序代码 ======== clc ;clear all; hold on; Fs=960; %采样频率 N=960; %采样点 n=0:N-1;t=n/Fs; %时间序列 A0=10; %载波信号振幅 A1=1; %调制信号振幅 fc=120; %载波信号频率 fm=30; %调制信号频率 f=n*Fs/N; %频率 w0=2*fc*pi; w1=2*fm*pi; Uc=A0*cos(w0*t); %载波信号 C1=fft(Uc); %对载波信号进行傅里叶变换 cxf=abs(C1); %进行傅里叶变换 % figure(1); hold on; subplot(6,2,1); plot(t,Uc); title('载波信号波形'); axis([0 0.1 -10 10]); subplot(6,2,2); plot(f(1:N/2),cxf(1:N/2));title('载波信号频谱'); % subplot(6,2,2) % f = fftshift(C1); % w = linspace(-Fs/2, Fs/2, N);%频率坐标，单位Hz % plot(w,abs(f)); % title('信号的频谱'); % xlabel('频率(Hz)'); mes=1+A1*cos(w1*t); %调制信号 C2=fft(mes); % 对调制信号进行傅里叶变换 zxc=abs(C2); subplot(6,2,3); plot(t,mes); title('调制信号');axis([0 0.5 0 2]); subplot(6,2,4); plot(f(1:N/2),zxc(1:N/2)); title('调制信号频谱'); axis([0 50 0 500]) Uam=modulate(mes,fc,Fs,'am');%AM 已调信号 C3=fft(Uam); % 对AM已调信号进行傅里叶变换 asd=abs(C3); % figure(3) subplot(6,2,5);plot(t,Uam); grid on; axis([0 0.5 -2 2]);title('AM已调信号波形'); subplot(6,2,6);plot(f(1:N/2),asd(1:N/2)),grid; title('AM已调信号频谱'); Dam=demod(Uam,fc,Fs,'am'); %对AM调制信号进行解调 C4=fft(Dam); % 对AM解调信号进行傅里叶变换 wqe=abs(C4); % figure(4) subplot(6,2,7); plot(t,Dam); grid on; title('AM解调信号波形'); axis([0 0.5 0 2]); subplot(6,2,8); plot(f(1:N/2),wqe(1:N/2)),grid; title('AM解调信号频谱'); k=awgn(Uam,10); %加大噪声，信噪比为10 pp=fft(k,960); zs=abs(pp); % figure(6) subplot(6,2,9);plot(t,k);axis([0 0.5 -2 2]) title('加噪声后得到AM信号时域波形');%加大噪声后得到AM信号时域波形 subplot(6,2,10);plot(f(1:N/2),zs(1:N/2)); title('加噪声后得到AM信号频域波形'); %加大噪声后得到AM信号频域波形 grid on; qaz=demod(k,fc,Fs,'am'); %加大噪声后解调 % figure(7) subplot(6,2,11);plot(t,qaz);title('加大噪声后解调得到信号时域波形'); %加大噪声后解调得到信号时域波形 axis([0 0.5 -1 2]) wsx=fft(qaz,960); edc=abs(wsx); subplot(6,2,12); plot(f(1:N/2),edc(1:N/2)); hold off; title('加大噪声后解调得到信号频域波形'); %加大噪声后解调得到信号频域波形 grid on; % -------------------------------------------------------------------- function Untitled_3_Callback(hObject, eventdata, handles) % hObject handle to Untitled_3 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) fm=10;fc=40; am=sqrt(2); %振幅 Fs=300; %采样频率Fs，载波频率fc，信号频率fm wc=2*pi*fc; wm=fm*2*pi; N=300; n=0:N-1; t=n/Fs; %时间序列 f=n*Fs/N; %频率 %基带信号时域 sm=am*cos(wm*t); subplot(321); plot(t,sm); title('基带信号'); xlabel('t'); axis([0 1 -2 2]); grid on %基带信号频域 S=fft(sm,300);%300点的fft SG=abs(S); subplot(322); plot(f(1:N/2),SG(1:N/2)); %SSB信号频域波形 xlabel('Frequency(HZ)'); title('基带信号频域波形 '); axis([0 100 0 300]); set(gca,'XTick',0:20:100); grid on; %SSB调制信号时域 %%s=modulate(sm,fc,Fs,'amssb'); s=modulate(sm,fc,Fs,'amssb'); %对调制信号进行调制 S=fft(s,300); SG=abs(S); subplot(323); plot(t,s); %SSB信号时域波形 title('SSB信号时域波形 '); xlabel('t'); grid on; subplot(324); plot(f(1:N/2),SG(1:N/2)); %SSB信号频域波形 xlabel('Frequency(HZ)'); title('SSB信号频域波形 '); axis([0 100 0 300]); set(gca,'XTick',0:20:100); grid on; %------------------------------------------------------------------------- %解调 fm=10;%信号频率fm fc=40;%载波频率fc am=sqrt(2); Fs=300; %采样频率Fs wc=2*pi*fc; wm=fm*2*pi; N=300; n=0:N-1; t=n/Fs; %时间序列 f=n*Fs/N; sm=am*cos(wm*t); s=modulate(sm,fc,Fs,'amssb'); sd=demod(s,fc,Fs,'amssb'); %对SSB信号进行解调 SD=fft(sd,300); SDG=abs(SD); subplot(325); plot(t,sd); %解调后的时域波形 title('解调后的时域波形'); xlabel('t'); axis([0 1 -2 2]); grid on; subplot(326); plot(f(1:N/2),SDG(1:N/2)); %解调后的频域波形 title('解调后的频域波形'); xlabel('Frequency(HZ)'); axis([0 100 0 300]); set(gca,'XTick',0:20:100); grid on; % -------------------------------------------------------------------- function Possion_Callback(hObject, eventdata, handles) % hObject handle to Possion (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) clear all; clc; U = rand; lamda = 2; % 这里给定值 T = 10; % 这里给