[問題] MATLAB響應圖問題已回收
clear;
clf;
W=[0.6858 -0.0549;0.0562 0.9414]; %K矩陣
N=20; %取樣時間
D=[-1 0;0 -1]; %E矩陣
A=[2 0;0 3]; %A矩陣
C=[0.1 0;0 0.1]; %C矩陣
W1=[-0.34 0;-0.1 0.1]; %W1矩陣
W2=[0.1 0.2;-0.15 -0.18]; %W2矩陣
W3=[0.41 -0.5;0.69 0.31]; %W3矩陣
Q=A+(W*D);
e0=[-1;2]; %初始狀態值
for k=0:N
if k==0
e_temp=e0;
else
t=pi;
d_=fix((1+3*sin(k*pi/2))/1);
tau_=fix((3+5*sin(k*pi/2))/1);
h_=fix((2+4*sin(k*pi/2))/1);
% if e_temp < k
% syms e_temp
% e_temp =
% -Q*e_temp+C*E(tau_)+W1*0.5*(abs(e_temp+[1;1])+abs(e_temp-[1;1]
% ))+W2*0.5*(abs(E(h_)+[1;1])+abs(E(h_)-[1;1]))+W3*0.5*(int(abs(
% e_temp+[1;1])+abs(e_temp-[1;1]),e_temp,t-d_,t))-W*e_temp ;
% else
%syms e_temp
e_temp=
-Q*e0+C*E(t-tau_)+W1*0.5*(abs(e0+[1;1])+abs(e0-[1;1]))+W2*0.5*(abs((t-h_)+[1;1])+abs((t-h_)-[1;1]))+W3*0.5*(int(abs(e_temp+[1;1])+abs(e_temp-[1;1]),e_temp,t-d_,t))-W*e0;
%e_temp
=W1*0.5*(abs(e_temp+[1;1])+abs(e_temp-[1;1]))+W2*0.5*(abs((t-h_)+[1;1])+abs((t-h_)-[1;1]))+W3*0.5*(int(abs(e_temp+[1;1])+abs(e_temp-[1;1]),e_temp,t-d_,t))
;
%e_temp =
-Q*e0+W1*0.5*(abs(e_temp+[1;1])+abs(e_temp-[1;1]))+W2*0.5*(abs((t-h_)+[1;1])+abs((t-h_)-[1;1]))+W3*0.5*(int(abs(e_temp+[1;1])+abs(e_temp-[1;1]),e_temp,t-d_,t))-W*e0
;
% e_temp =
-Q*e_temp+C*E(tau_)+W1*0.5*(abs(e_temp+[1;1])+abs(e_temp-[1;1]))+W2*0.5*(abs(E(h_)+[1;1])+abs(E(h_)-[1;1]))+W3*0.5*(int(abs(e_temp+[1;1])+abs(e_temp-[1;1]),e_temp,t-d_,t))-W*e_temp
;
% end
e0=e_temp;
end
if k == 0
E = e_temp';
else
E = [E;e_temp'];
end
end
figure(1)
plot(0:1:N,E(:,1) ,0:1:N,E(:,2))
% title('Time response of discrete-time recurrent neural network')
xlabel('time(sec)')
ylabel('amplitude')
這是有關類神經網路的響應圖程式，其中C*E(t-tau_)加上後，程式就跑不出來了，有哪位大大可以幫我解答一下!!
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 125.224.176.161
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章
