@KolyaIce

Нужно построить спектр функции .
Проблема в том, что график суммы двух синусоид в скобках строится, а вот при умножении на скосинусоиду после скобки выдает ошибку:
Error using *
Inner matrix dimensions must agree.
Error in Spectr2 (line 18)
Signal=(Ak+A1*cosd((F1*360).*T+Phi1)+A2*cosd((F2*360).*T+Phi2))*cosd((F0*360).*T +Phi0);% Массив сигнала
Помогите разобраться!
Вот исходный код (за основу брал тот, что нашел в интернете http://habrahabr.ru/post/112068/ - изначально там сумма двух косинусоид)

clear all% Очистка памяти
%% Параметры
Tm=2;% Длина сигнала (с)
Fd=512;% Частота дискретизации (Гц)
Ak=1;% Постоянная составляющая
A1=0.7;% Амплитуда сигнала 1
A2=0.5;% Амплитуда сигнала 2
A0=1;% Амплитуда сигнала 0
F0=6000000; % Частота сигнала 0 (Гц)
F1=1000;% Частота сигнала 1 (Гц)
F2=2000;% Частота сигнала 2 (Гц)
Phi0=0; % Начальная фаза сигнала 0 (Градусов)
Phi1=0;% Начальная фаза сигнала 1 (Градусов)
Phi2=0;% Начальная фаза сигнала 2 (Градусов)
FftL=1024;% Количество линий Фурье спектра
%% Генерация рабочих массивов
T=0:1/Fd:Tm;% Массив отсчетов времени
Signal=(Ak+A1*cosd((F1*360).*T+Phi1)+A2*cosd((F2*360).*T+Phi2))*cosd((F0*360).*T +Phi0);% Массив сигнала !!!! ЗДЕСЬ ГДЕ ТО ОШИБКА !!!!!!
%% Спектральное представление сигнала
FftS=abs(fft(Signal,FftL));% Амплитуды преобразования Фурье сигнала
FftS=2*FftS./FftL;% Нормировка спектра по амплитуде
FftS(1)=FftS(1)/2;% Нормировка постоянной составляющей в спектре
FftSh=abs(fft(Signal,FftL));% Амплитуды преобразования Фурье смеси сигнал+шум
FftSh=2*FftSh./FftL;% Нормировка спектра по амплитуде
FftSh(1)=FftSh(1)/2;% Нормировка постоянной составляющей в спектре
%% Построение графиков
subplot(2,1,1);% Выбор области окна для построения
plot(T,Signal);% Построение сигнала
title('Сигнал');% Подпись графика
xlabel('Время (с)');% Подпись оси х графика
ylabel('Амплитуда');% Подпись оси у графика
subplot(2,1,2);% Выбор области окна для построения
plot(T,Signal);% Построение смеси сигнал+шум
title('Сигнал+шум');% Подпись графика
xlabel('Время (с)');% Подпись оси х графика
ylabel('Амплитуда');% Подпись оси у графика

F=0:Fd/FftL:Fd/2-1/FftL;% Массив частот вычисляемого спектра Фурье
figure% Создаем новое окно
subplot(2,1,1);% Выбор области окна для построения
plot(F,FftS(1:length(F)));% Построение спектра Фурье сигнала
title('Спектр сигнала');% Подпись графика
xlabel('Частота (Гц)');% Подпись оси х графика
ylabel('Амплитуда');% Подпись оси у графика
subplot(2,1,2);% Выбор области окна для построения
plot(F,FftSh(1:length(F)));% Построение спектра Фурье сигнала
title('Спектр сигнала');% Подпись графика
xlabel('Частота (Гц)');% Подпись оси х графика
ylabel('Амплитуда');% Подпись оси у графика

0

@SSC

Попробуй так

3

@KolyaIce

Помогло! Спасибо! Не думал, что все так просто! Я изначально ставил точку между скобками, а вот перед периодом Т не убрал!

0

Зосима

для ознакомления
СЕКРЕТНО! Зловещая точка претыкания

Как оформлять код в сообщениях

0

@KolyaIce

Теперь у меня другая проблема в частотной области. Центральная частота должна быть порядка 10^5 Гц, а на графике сотни Гц. Что то наверное в спектральном представлении сигнала напутал.
Учитывая Signal=5*10^-4*(Ak+A1*cosd((2*pi*1000)*T)+A2*cosd((4*pi*1000)*T)).*cosd((2.56*pi*10^5)*T)
Я новичок, не бейте сильно))

0

@SSC

1 Определись в градусах или радианах ты проводишь вычисления (cos-радианы cosd-градусы)
2 Частота дискретизации (Fd у Вас в программе) должна быть не менее чем в 3 раза больше (а желательно больше), чем несущая

3

@KolyaIce

Спасибо! Подправил
Fd=4194304 - Частота дискретизации (Гц)
FftL=524288 - Количество линий Фурье спектра, иначе спектр очень размазанный
Signal=5*10^-4*(Ak+A1*cosd((F1*360)*T)+A2*cosd((F2*360)*T)).*cosd((F0*360)*T) - Теперь здесь частота в Герцах (модулирующие -1 кГц и 2 кГц, несущая - 465 кГц)

Теперь все как надо!

0