Error using fft on code
% Not sure what I am doing wrong to be getting an error with the data type
% Load data from Signal.mat file
load(‘Spec.mat’);
% Number of rows (process realizations)
num_rows = size(Spec, 1);
% Number of samples per row
num_samples = size(Spec, 2);
% Calculate the power spectral density (PSD) for each row
psd_rows = zeros(num_rows, num_samples);
for i = 1:num_rows
psd_rows(i, 🙂 = (abs(fftshift(fft(Spec(i, :)))).^2) / num_samples; <————- error here
end
% Average the PSDs of all rows
psd_average = mean(psd_rows, 1);
% Sampling frequency (assuming 40 kHz)
Fs = 40e3;
% Frequency resolution
df = Fs / num_samples;
% Frequency axis
f_axis = (-num_samples/2 : num_samples/2 – 1) * df;
% Plot individual and composite PSDs
figure;
hold on;
for i = 1:num_rows
plot(f_axis, fftshift(psd_rows(i, :)));
end
plot(f_axis, fftshift(psd_average), ‘k’, ‘LineWidth’, 2);
xlabel(‘Frequency (Hz)’);
ylabel(‘Power/Frequency’);
title(‘Power Spectral Density (PSD)’);
legend(‘Individual PSDs’, ‘Composite PSD’);% Not sure what I am doing wrong to be getting an error with the data type
% Load data from Signal.mat file
load(‘Spec.mat’);
% Number of rows (process realizations)
num_rows = size(Spec, 1);
% Number of samples per row
num_samples = size(Spec, 2);
% Calculate the power spectral density (PSD) for each row
psd_rows = zeros(num_rows, num_samples);
for i = 1:num_rows
psd_rows(i, 🙂 = (abs(fftshift(fft(Spec(i, :)))).^2) / num_samples; <————- error here
end
% Average the PSDs of all rows
psd_average = mean(psd_rows, 1);
% Sampling frequency (assuming 40 kHz)
Fs = 40e3;
% Frequency resolution
df = Fs / num_samples;
% Frequency axis
f_axis = (-num_samples/2 : num_samples/2 – 1) * df;
% Plot individual and composite PSDs
figure;
hold on;
for i = 1:num_rows
plot(f_axis, fftshift(psd_rows(i, :)));
end
plot(f_axis, fftshift(psd_average), ‘k’, ‘LineWidth’, 2);
xlabel(‘Frequency (Hz)’);
ylabel(‘Power/Frequency’);
title(‘Power Spectral Density (PSD)’);
legend(‘Individual PSDs’, ‘Composite PSD’);Â % Not sure what I am doing wrong to be getting an error with the data type
% Load data from Signal.mat file
load(‘Spec.mat’);
% Number of rows (process realizations)
num_rows = size(Spec, 1);
% Number of samples per row
num_samples = size(Spec, 2);
% Calculate the power spectral density (PSD) for each row
psd_rows = zeros(num_rows, num_samples);
for i = 1:num_rows
psd_rows(i, 🙂 = (abs(fftshift(fft(Spec(i, :)))).^2) / num_samples; <————- error here
end
% Average the PSDs of all rows
psd_average = mean(psd_rows, 1);
% Sampling frequency (assuming 40 kHz)
Fs = 40e3;
% Frequency resolution
df = Fs / num_samples;
% Frequency axis
f_axis = (-num_samples/2 : num_samples/2 – 1) * df;
% Plot individual and composite PSDs
figure;
hold on;
for i = 1:num_rows
plot(f_axis, fftshift(psd_rows(i, :)));
end
plot(f_axis, fftshift(psd_average), ‘k’, ‘LineWidth’, 2);
xlabel(‘Frequency (Hz)’);
ylabel(‘Power/Frequency’);
title(‘Power Spectral Density (PSD)’);
legend(‘Individual PSDs’, ‘Composite PSD’); fft MATLAB Answers — New Questions
​