Index exceeds matrix dimensions.

clc;
clear;

%Исходные данные
F=100; %частота следования импульсов, Гц
alpha=0.5; %Коэффициент скругления
U0=1; % Амплитуда несущего колебания, В
f0=100000; %частота несущего колебания, Гц

Fs=10*f0; %частота дискретизации, Гц
T=1/F; %Период следования импульсов, с
Ts=1/Fs; %Период дискретизации, с

f=140:0.1:160; % массив частот
M=length(f);
j=complex(0, 1);

%Добавление шума
SNR=-15; % добавление шума, дБ

% 1. Расчёт импульса приподнятого косинуса
delay=3; %задержка
g=rcosine(F, 200*F, ‘sqrt’, alpha, delay); %импульс приподнятого косинуса
g=g/max(g);
N1=length(g);
t1=0:(1/(200*F)):(N1-1)*(1/(200*F));

% 2. Синтез случайного кода
Nsym=100; %число символов
p=ceil(2.*rand(1,Nsym));
p=p-1;

%———————————-Модуляция DQPSK————
% 3. Расчёт символов a
a=zeros(1, Nsym);
for n=1:Nsym %Вычисление A и B
if p(n)==0
a(n)=-1;
end
if p(n)==1
a(n)=1;
end
end
n=1:Nsym;

figure(1)
stem(n, a);
%xlim([0, 20]);
ylim([-2, 2]);
xlabel(‘n’);
ylabel(‘a’);
grid on

% 4. Синтез модулирующего сигнала
N=Nsym*N1;
t=0:Ts:(N-1)*Ts;
Sa=zeros(1,N);
for n=1:Nsym
for k=1:N1
Sa(k+(n-1)*200*F)=a(n).*g(k)+Sa(k+(n-1)*200*F);
end
end

figure(2)
plot(t, Sa);
title(‘График модулирующего сигнала символов a’);
xlim([0, 10]);
ylim([-2, 2]);
xlabel(‘t, с’);
ylabel(‘Sa(t), В’);
grid on

% 5. Балансная модуляция
Ug=U0*cos(2*pi*f0*t); %несущее колебание
Uw=Sa.*Ug; %радиосигнал
P=mean(Uw.*Uw); % мощность в Вт
PF=10*log10((max(Uw.*Uw))/P); %пик-фактор
Ps=10*log10(P); % мощность в дБ
Noise=wgn(1, N, Ps-SNR);
Uw=Uw+Noise; % добавление шума

figure(3)
plot(t, Uw);
xlim([0, 20]);
ylim([-2, 2]);
xlabel(‘t, с’);
ylabel(‘Uw(t), В’);
grid on

%———————————-Демодуляция BPSK————
% 6. Балансная демодуляция
t=0:Ts:(N-1)*Ts;
Uda=Uw.*cos(2*pi*f0*t);
Sfa=filter(g, 1, Uda);
Sfa=Sfa/max(Sfa)*max(Sa);
t=0:(1/(200*F)):(N-1)*(1/(200*F));

figure(4)
subplot(2, 1, 1);
plot(t, Sfa);
title(‘График демодулированного сигнала a’);
xlim([0, 0.23]);
ylim([-2, 2]);
xlabel(‘t, с’);
ylabel(‘Sa(t), В’);
grid onclc;
clear;

Fs=10*f0; %частота дискретизации, Гц
T=1/F; %Период следования импульсов, с
Ts=1/Fs; %Период дискретизации, с

f=140:0.1:160; % массив частот
M=length(f);
j=complex(0, 1);

%Добавление шума
SNR=-15; % добавление шума, дБ

% 2. Синтез случайного кода
Nsym=100; %число символов
p=ceil(2.*rand(1,Nsym));
p=p-1;

figure(1)
stem(n, a);
%xlim([0, 20]);
ylim([-2, 2]);
xlabel(‘n’);
ylabel(‘a’);
grid on

% 4. Синтез модулирующего сигнала
N=Nsym*N1;
t=0:Ts:(N-1)*Ts;
Sa=zeros(1,N);
for n=1:Nsym
for k=1:N1
Sa(k+(n-1)*200*F)=a(n).*g(k)+Sa(k+(n-1)*200*F);
end
end

figure(2)
plot(t, Sa);
title(‘График модулирующего сигнала символов a’);
xlim([0, 10]);
ylim([-2, 2]);
xlabel(‘t, с’);
ylabel(‘Sa(t), В’);
grid on

figure(3)
plot(t, Uw);
xlim([0, 20]);
ylim([-2, 2]);
xlabel(‘t, с’);
ylabel(‘Uw(t), В’);
grid on

Fs=10*f0; %частота дискретизации, Гц
T=1/F; %Период следования импульсов, с
Ts=1/Fs; %Период дискретизации, с

f=140:0.1:160; % массив частот
M=length(f);
j=complex(0, 1);

%Добавление шума
SNR=-15; % добавление шума, дБ

% 2. Синтез случайного кода
Nsym=100; %число символов
p=ceil(2.*rand(1,Nsym));
p=p-1;

figure(1)
stem(n, a);
%xlim([0, 20]);
ylim([-2, 2]);
xlabel(‘n’);
ylabel(‘a’);
grid on

% 4. Синтез модулирующего сигнала
N=Nsym*N1;
t=0:Ts:(N-1)*Ts;
Sa=zeros(1,N);
for n=1:Nsym
for k=1:N1
Sa(k+(n-1)*200*F)=a(n).*g(k)+Sa(k+(n-1)*200*F);
end
end

figure(2)
plot(t, Sa);
title(‘График модулирующего сигнала символов a’);
xlim([0, 10]);
ylim([-2, 2]);
xlabel(‘t, с’);
ylabel(‘Sa(t), В’);
grid on

figure(3)
plot(t, Uw);
xlim([0, 20]);
ylim([-2, 2]);
xlabel(‘t, с’);
ylabel(‘Uw(t), В’);
grid on

figure(4)
subplot(2, 1, 1);
plot(t, Sfa);
title(‘График демодулированного сигнала a’);
xlim([0, 0.23]);
ylim([-2, 2]);
xlabel(‘t, с’);
ylabel(‘Sa(t), В’);
grid on matlab MATLAB Answers — New Questions

Cart

Cart