2017-01-09 5 views
-3

L'erreurComment puis-je déboguer le code Matlab suivant?

fichier: parameter_estimation_1.m ligne: 8 Colonne: 10

expression inattendue Matlab.

cultures quand je lance le code suivant Matlab:

T = 0:0.25:5; % time vector (row) 
T = T'; % time vector (column) 
seed = [3;0.5]; % seed for noise generation 
randn('state',seed); % using the same seed each time 
uu = 0.5 1 0.25*randn(length(T),1); % mean of 0.5 with variance 
% of 0.25 
U = 2*round(uu)-1; % creates PRBS with -1 and 1 values % 
sim('est_vdv'); % runs simulation of linear van de vusse % diagram 
figure(1); % plot input-output data 
subplot(2,1,1),plot(tp,yp,'k',t,y,'ko'); 
xlabel('time, min'), ylabel('y') 
title('PRBS estimation example') 
subplot(2,1,2),plot(tp,up,'k'); xlabel('time, min'),ylabel('u') 
axis([0 5 -1.1 1.1]) 
% % generate phi matrix for estimation 
for j = 4:22; 
phi(j-3,:) = [y(j-2) y(j-3) u(j-2) u(j-3)]; 
end 
% 
theta = inv(phi'*phi)*phi'*y(3:21) % estimate parameters 
num = [theta(3) theta(4)]; % numerator of discrete transfer function 
den = [1 -theta(1) -theta(2)]; % denominator of discrete transfer function 
sysd = tf(num,den,0.25) % create discrete tf object 
tzero(sysd) % calculate zeros 
pole(sysd) % calculate poles 
syszpk = zpk(sysd) % zero-pole-k form 

Ce code est censé fonctionner en tandem avec un modèle Simulink intitulé « est_vdv » pour estimer les paramètres d'un modèle. Comment dois-je gérer cette erreur?

+2

vous voulez dire cette ligne? 'Subplot (2,1,1), plot (tp, yp, 'k', t, y, 'ko');' Il y a un coma au lieu d'un point virgule –

+0

ou this un: «uu = 0,5 1 0,25 * randn (longueur (T), 1); (les parenthèses manquantes) –

+0

En général, la plupart des utilisateurs de Stack Overflow n'aiment pas beaucoup les questions [mendier du poisson] (https://en.wiktionary.org/wiki/give_a_man_a_fish_and_you_feed_him_for_a_day;_teach_a_man_to_fish_and_you_feed_him_for_a_lifetime) - vous devriez apprendre à * attraper le vôtre *, dans ce cas en apprenant à déboguer votre propre code. –

Répondre

-1

Merci amis pour vos suggestions, j'ai été en mesure de comprendre ce qui a mal tourné avec la 5ème ligne, il aurait dû être

uu=0.5+0.25.*randn(length(T),1) 

Désolé, il a été à tort indiqué que l'erreur était dans la 8e ligne .