function varargout = untitled(varargin)
% UNTITLED M-file for untitled.fig
% UNTITLED, by itself, creates a new UNTITLED or raises the existing
% singleton*.
%
% H = UNTITLED returns the handle to a new UNTITLED or the handle to
% the existing singleton*.
%
% UNTITLED('CALLBACK',hObject,eventData,handles,...) calls the local
% function named CALLBACK in UNTITLED.M with the given input arguments.
%
% UNTITLED('Property','Value',...) creates a new UNTITLED or raises the
% existing singleton*. Starting from the left, property value pairs are
% applied to the GUI before untitled_OpeningFunction gets called. An
% unrecognized property name or invalid value makes property application
% stop. All inputs are passed to untitled_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 untitled
% Last Modified by GUIDE v2.5 09-Jun-2006 19:08:35
% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name', mfilename, ...
'gui_Singleton', gui_Singleton, ...
'gui_OpeningFcn', @untitled_OpeningFcn, ...
'gui_OutputFcn', @untitled_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 untitled is made visible.
function untitled_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 untitled (see VARARGIN)
% Choose default command line output for untitled
handles.output = hObject;
% Update handles structure
guidata(hObject, handles);
% UIWAIT makes untitled wait for user response (see UIRESUME)
% uiwait(handles.figure1);
% --- Outputs from this function are returned to the command line.
function varargout = untitled_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 Untitled_1_Callback(hObject, eventdata, handles)
% hObject handle to Untitled_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)
set(handles.text5, 'String', '点击Yes关闭当前窗口');
answer=questdlg('您要关闭窗口吗?','关闭','Yes','No','Cancel','Yes');%打开提问对话框,确认是否关闭当前窗口
switch answer %如果选择OK则关闭当前窗口
case 'Yes'
close
case 'No'
case 'Cancel'
end
% --------------------------------------------------------------------
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)
set(handles.text5, 'String', '在对话框中选择一种方式后点击OK');
a=handles.metricdata.density%把输入的值赋给a
if a>1 %判断输入是否合法
errordlg('输入必须在0-1之间','Error');
else if a errordlg('输入必须在0-1之间','Error');
else
str={'有噪声','无噪声'};
[Selection,ok]=listdlg('ListString',str,'Name','选择一种方式','PromptString','请选择','SelectionMode','Single');%打开菜单对话框,用以选择所要仿真的类型,对于升余弦滚降系统,分为"有噪声"和"无噪声"
alpha=a;%滚降系数
N=2^14; %采样点数
L=32; %每码元采样点
M=N/L; %码元数
Rb=2; %码速率
Ts=1/Rb; %码元间隔
dt=Ts/L; %时域采样间隔
df=1/(N*dt); %频域采样间隔
T=N*dt; %截短时宽
Bs=N*df/2; %频带宽度
Na=4; %画眼图一次四个码元
t=[-T/2+dt/2:dt:T/2]; %时域坐标
f=[-Bs+df/2:df:Bs]; %频域坐标
g1=sin(pi*t/Ts)./(pi*t/Ts);
g2=cos(alpha*pi*t/Ts)./(1-(2*alpha*t/Ts).^2);
g=g1.*g2;%升余弦脉冲波形
G=t2f(g,dt); %系统总的频谱
a=sign(randn(1,M));
imp=zeros(1,N);
imp(L/2:L:N)=a/dt;%产生冲击序列
S=t2f(imp,dt).*G;%升余弦信号的傅立叶变换
if ok==1
switch Selection%使用switch语句对菜单对话框进行控制
case 1%当选择"有噪声"时,执行以下的函数
s=f2t(S,dt);%升余弦信号的时域波形
s=real(s)+0.2*randn(1,N);
set(handles.text2, 'String', '噪声情况下的眼图');
axes(handles.axes2);%在第二个坐标系中绘制噪声情况下的眼图
cla;
hold on
grid
%画眼图
tt=[0:dt:Na*L*dt];
for jj=1:Na*L:N-Na*L
plot(tt,s(jj:jj+Na*L))
end
axes(handles.axes1);%在第一个坐标系中绘制系统总的频谱
cla;
plot(f,abs(G))
grid
axis([-2,2,min(abs(G)),max(abs(G))]);
case 2%当选择"无噪声"时,执行以下的函数
set(handles.text2, 'String', '无噪声情况下的眼图');
s=f2t(S,dt);%升余弦信号的时域波形
s=real(s);
axes(handles.axes2);%在第二个坐标系中绘制无噪声情况下的眼图
cla;
hold on
grid
%画眼图
tt=[0:dt:Na*L*dt];
for jj=1:Na*L:N-Na*L
plot(tt,s(jj:jj+Na*L))
end
axes(handles.axes1);%在第一个坐标系中绘制系统总的频谱
cla;
plot(f,abs(G))
grid
axis([-2,2,min(abs(G)),max(abs(G))]);
end;
else
end;
end
end
% --------------------------------------------------------------------
function Untitled_5_Callback(hObject, eventdata, handles)
% hObject handle to Untitled_5 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
a=handles.metricdata.density%把输入的值赋给a
set(handles.text5, 'String', '在对话框中选择一种方式后点击OK');
if a>1%判断输入是否合法
errordlg('输入必须在0-1之间','Error');
else if a errordlg('输入必须在0-1之间','Error');
else
str={'均衡前','均衡后'};
[Selection,ok]=listdlg('ListString',str,'Name','选择一种方式','PromptString','请选择','SelectionMode','Single');%打开菜单对话框,用以选择所要仿真的类型,对于多径传播,分为"均衡前"和"均衡后"
alpha=a;
N=2^14; %采样点数
L=32; %每码元采样点数
M=N/L; %码元数
Rb=2; %码速率
Ts=1/Rb; %码元间隔
dt=Ts/L; %时域采样间隔
df=1/(N*dt); %频域采样间隔
d=16; %半个码元采样点数
T=N*dt; %截短时宽
Bs=N*df/2; %频带宽度
Na=4; %画眼图一次四个码元
t=[-T/2+dt/2:dt:T/2]; %时域坐标
f=[-Bs+df/2:df:Bs]; %频域坐标
g1=sin(pi*t/Ts)./(pi*t/Ts);
g2=cos(alpha*pi*t/Ts)./(1-(2*alpha*t/Ts).^2);
g=g1.*g2;%升余弦脉冲波形
Cf=1-0.5*exp(-2*i*pi*f*d); %信道传输函数
a=sign(randn(1,M));
imp=zeros(1,N);
imp(L/2:L:N)=a/dt;%产生冲击序列
if ok==1
switch Selection%使用switch语句对菜单对话框进行控制
case 1%当选择"均衡前"时,执行以下的函数
set(handles.text2, 'String', '多径均衡前的眼图');
G=t2f(g,dt).*Cf; %系统总的频谱
S=t2f(imp,dt).*G; %序列在频域通过系统
s=f2t(S,dt); %转换到时域
s=real(s)+0*randn(1,N); %可加入随机噪声
%预设眼图的位置
axes(handles.axes2);%在第三个坐标系中绘制多径均衡前的眼图
cla;
hold on
grid
%画眼图
tt=[0:dt:Na*L*dt];
for jj=1:Na*L:N-Na*L
plot(tt,s(jj:jj+Na*L))
end
axes(handles.axes1);%在第一个坐标系中绘制系统总的频谱
cla;
plot(f,abs(G))
axis([-2,2,min(abs(G)),max(abs(G))]);
case 2%当选择"均衡后"时,执行以下的函数
set(handles.text2, 'String', '多径均衡后的眼图');
G=t2f(g,dt).*Cf.*1./Cf; %系统总的频谱
S=t2f(imp,dt).*G; %序列在频域通过系统
s=f2t(S,dt); %转换到时域
s=real(s);
axes(handles.axes2);%在第三个坐标系中绘制多径均衡后的眼图
cla;
hold on
grid
%画眼图
tt=[0:dt:Na*L*dt];
for jj=1:Na*L:N-Na*L
plot(tt,s(jj:jj+Na*L))
end
axes(handles.axes1);%在第一个坐标系中绘制系统总的频谱
cla;
plot(f,abs(G))
grid
axis([-2,2,min(abs(G)),max(abs(G))]);
end;
else
end
end
end
% --------------------------------------------------------------------
function run_Callback(hObject, eventdata, handles)
% hObject handle to run (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% --------------------------------------------------------------------
% --------------------------------------------------------------------
function close_Callback(hObject, eventdata, handles)
% hObject handle to close (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% --------------------------------------------------------------------
function Untitled_6_Callback(hObject, eventdata, handles)
% hObject handle to Untitled_6 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% --------------------------------------------------------------------
function Untitled_7_Callback(hObject, eventdata, handles)
% hObject handle to Untitled_7 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
set(handles.text5, 'String', '在对话框中选择一种方式后点击OK');
a=handles.metricdata.density%把输入的值赋给a
if a>1%判断输入是否合法
errordlg('输入必须在0-1之间','Error');
else if a errordlg('输入必须在0-1之间','Error');
else
str={'不加网孔','加网孔'};
[Selection,ok]=listdlg('ListString',str,'Name','选择一种方式','PromptString','请选择','SelectionMode','Single');%打开菜单对话框,用以选择所要仿真的类型,对于多径传播,分为"有均衡"和"无均衡"
alpha=a;
set(handles.text1, 'String', '网孔滤波特性');
N=2^14; %采样点数
L=32; %每码元采样点数
M=N/L; %码元数
Rb=2; %码速率
Ts=1/Rb; %码元间隔
dt=Ts/L; %时域采样间隔
df=1/(N*dt); %频域采样间隔
T=N*dt; %截短时宽
Bs=N*df/2; %频带宽度
Na=4; %画眼图一次四个码元
t=[-T/2+dt/2:dt:T/2]; %时域坐标
f=[-Bs+df/2:df:Bs]; %频域坐标
a=round(rand(1,M));
rz=zeros(1,N);
for tt=1:L
rz(tt+[0:M-1]*L)=a;
end
RZ=t2f(rz,dt);
WK=(pi*f*Ts)./sin(pi*f*Ts);%网孔均衡
g1=sin(pi*t/Ts)./(pi*t/Ts);
g2=cos(alpha*pi*t/Ts)./(1-(2*alpha*t/Ts).^2);
g=g1.*g2;%升余弦脉冲波形
H=t2f(g,dt); %系统总的频谱
G= RZ;%G为不归零脉冲序列
S2=G.*WK.*H;%加网孔,升余弦滚降
S3=G.*H;%不加网孔,直接滚降
s2=real(f2t(S2,dt));
s3=real(f2t(S3,dt));
tt=[0:dt:Na*L*dt];
if ok==1%判断对话框中OK键是否被按下
switch Selection
case 1
set(handles.text2, 'String', '不加网孔时的眼图');
axes(handles.axes1);%在第一个坐标系中绘制网孔滤波特性曲线
cla;
plot(f,abs(WK))
grid
axis([0,1/Ts,0,10])
%画眼图
axes(handles.axes2);%在第二个坐标系中绘制不加网孔时的眼图
cla;
for jj=1:Na*L:N/2
plot(tt,s3(jj:jj+Na*L));
hold on
end
case 2
set(handles.text2, 'String', '加网孔时的眼图');
axes(handles.axes1);%在第一个坐标系中绘制网孔滤波特性曲线
cla;
plot(f,abs(WK))
grid
axis([0,1/Ts,0,10])
%画眼图
axes(handles.axes2);%%在第二个坐标系中绘制加网孔时的眼图
cla;
for jj=1:Na*L:N/2
plot(tt,s2(jj:jj+Na*L));
hold on
end
end
else
end
end
end
function edit3_Callback(hObject, eventdata, handles)
% hObject handle to edit3 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hints: get(hObject,'String') returns contents of edit3 as text
% str2double(get(hObject,'String')) returns contents of edit3 as a double
density = str2double(get(hObject, 'String'));%获取输入的滚降系数的值,将其赋给density
handles.metricdata.density = density;
guidata(hObject,handles)
% --- Executes during object creation, after setting all properties.
function edit3_CreateFcn(hObject, eventdata, handles)
% hObject handle to edit3 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called
% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end
