I get warning: Matrix is close to singular or badly scaled. results may be inaccurate. RCOND=6.564544e-61.

Question

wenlong Ma 2018-8-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/413052-i-get-warning-matrix-is-close-to-singular-or-badly-scaled-results-may-be-inaccurate-rcond-6-56454

I don't know what's wrong with my code, I attached my code! thanks a lot for your help!!!

load('MessdatenMadrid.mat'); nfft = 2^12; % 10 kHz Samplerate der Daten [TF_data, TF_w] = tfestimate(X_Achse,Y_Achse, hann(nfft), nfft/2, nfft, 1e4); % frequency-response (frd) model object FRF = frd(TF_data, TF_w);

% Der Rest ist deshalb eher didaktisch und f?r die Anschauung % Erzeuge kontinuierliche Zeit ?bertragungsfunktion (gebrochen rationales Polynom), % die einem bestimmten komplexen Frequenzgang entspricht rng = 1:length(TF_w) ; % hier evtl. obere Frequenzen abschneiden n = 9; % Anzahl Nullstellen, sollte um 1 kleiner sein als m = 10; % Anzahl Pole ? Diese hier variieren, so dass ?bererinstimmung zwischen TF_Data und Fit bis ~ 1 kHz gut [Num,Den] = invfreqs(TF_data(rng), TF_w(rng), n, m);

Fit = idpoly([], Num, [], [], Den, [], 0); % Polynom zum Plotten und Vergleichen mit TF_Data Root_Fit_Den = roots(Den); % komplexe Pole Root_Fit_Nom = roots(Num); % komplexe zero Position