Issue in finding Power Spectral Density of fitted curve

1 次查看（过去 30 天）
Rohan Gupta
Rohan Gupta 2017-7-4
编辑： dpb 2017-7-7
I have a respiration signal which is sinusoidal in shape. I have fitted the curve using Fourier series. Since the fitted curve is in cfit format, I am not able to find it's Power Spectral Density.
Can anyone help me out with this.

请先登录，再回答此问题。

采纳的回答

dpb
dpb 2017-7-4
编辑：dpb 2017-7-5
Evaluate the fit at the desired frequency for however long you have data and pass that result to fft
ADDENDUM: See the comment below for specific, but just use
yhat=resp_f(t);
to get predicted y at times t
  2 个评论
Rohan Gupta
Rohan Gupta 2017-7-6
Thanks, this has really solved my problem.
Just one more query, you have mentioned below that, I know the frequency content of the fitted signal from the series fit. I would like to know as how can I know the frequency of fitted signal without further calculation.
dpb
dpb 2017-7-6
编辑：dpb 2017-7-7
If you did a Fourier fit, there's no content in the fitted signal other than that in the specific terms in the fitted equation and they're all known identically.
ADDENDUM
If you subsequently do actually calculate a PSD by either method, you'll see the peaks at the terms in the fitted model.
The spectrum will, however, not be single spikes at only those frequencies but will be smeared into adjacent frequency bins. These are simply artifacts, however, not actual "real" energy and are due to the windowing and symmetry effects inherent in the computation of a sampled signal.

请先登录，再进行评论。

更多回答（1 个）

Image Analyst
Image Analyst 2017-7-4
Did you pass the signal into pwelch()? (Requires Signal Processing Toolbox.)
  2 个评论
Rohan Gupta
Rohan Gupta 2017-7-5
编辑：dpb 2017-7-5
Yes, I have tried that. When we give curve fitted signal (which is in cfit format) as an input to pwelch(), it returns an error. It would be really greatful if you can help me out with this.
Following is the code and also the respiratory signal which is used is attached.
load '0009_8min.mat';
t = [1:1:1800]./300;
Fs = 300;
resp = signal.co2.y;
resp_signal = resp(1:1800,:);
[resp_f,gof,output] = fit(t',resp_signal,'fourier8'); % Curve fitted respiration signal
plot(resp_signal);
figure
plot(resp_f);
[pxx2,w2] = pwelch(resp_f,[],[],[],Fs);
The returned error is:
Error in data_read (line 50)
[pxx2,w2] = pwelch(resp_f,[],[],[],Fs);
Reference to non-existent field 'internal'.
dpb
dpb 2017-7-5
编辑：dpb 2017-7-5
But you now already know what the frequency content of the fitted signal is from the series fit anyway????
Still, pwelch isn't enabled to process a cfit object, evaluate the fit and pass that...
t = [1:1:1800]./300;
[resp_f,gof,output] = fit(t',resp_signal,'fourier8');
[pxx2,w2] = pwelch(resp_f(t),[],[],[],Fs);
...

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Spectral Estimation 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by