I get warning: Matrix is close to singular or badly scaled. results may be inaccurate. RCOND=6.564544e-61.
2 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
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
0 Kommentare
Antworten (0)
Siehe auch
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!