Выпуск MATLAB для

Question

Рассчитайте периодограмму, используя переменные a и b, для последовательности частот f, между 1 и 20, с интервалами 0,2. Для каждого значения частоты, найти значение, через уравнение

T = (1/(4*pi*f))*atan(sum(sin(4*pi*f*t(i))))/(sum(cos(4*pi*f*t(i))) 

Затем вычислить значение periodgram для частоты f, P(f), через уравнение:

P(f) = ((sum(MasterYi(i)*sin(2*pi*f(t(i)-T))))^2)/(sum(sin^2(2*pi*f(t(i)-T))))+... ((sum(MasterYi(i)*cos(2*pi*f(t(i)-T))))^2)/(sum(cos^2(2*pi*f(t(i)-T)))) 

Где,

MasterYi = a * sin(2*pi*f(ti-T))+b*cos(2*pi*f(ti-T)) 

У меня есть этот код:

for f = 1:0.2:20 
%t is a value that we know 
t = 1:82; 
T(:,f) = (1./(4*pi.*f)).*tan(sum(sin(4*pi*f*t(:,f))))./(sum(cos(4*pi*f.*t(:,f)))) 

end 

На данный момент мы просто пытаемся получить значения для T, но он говорит:

индекс должен быть реальным положительные целые числа или логические выражения.

Answer 1

Я очень смущен относительно того, что происходит. Я делаю предположение, что t(:,f) в вопросе является ошибкой. С помощью этого предположения один из способов решить проблему стоимости T выглядит следующим образом:

% create a function handle to the T equation 
T = @(f,t) (1/(4*pi*f))*atan(sum(sin(4*pi*f*t)))/(sum(cos(4*pi*f*t))); 
% create the t vector 
t = 0:82; 
% create the f vector 
f = 1:.02:20; 
% pre-allocate the results Tresults 
Tresults = zeros(length(f),1); 
% loop over f finding T 
for x = 1:length(f) 
    Tresults(x) = T(f(x), t); 
end 

Проблема в том, нужно использовать целое число в качестве индекса для ваших результатов. f не может быть индексом, поскольку он содержит значения, такие как 1.2, 1.4, 1.6, которые не являются целыми числами. Вы можете следовать этому шаблону и распространить его на остальные ваши уравнения. Я бы продолжил, но в описании проблемы отсутствуют переменные, например ti, a, b. Есть крутые способы решить это, используя bsxfun, но это слишком продвинуто на этом этапе.