智能控制（神经网络）作业

智控制作业
学生姓名： 学 号： 专业班级：
72 采BP网路RBF网路DRNN网路逼线性象分进行matlab仿真
（） 采BP网络仿真
网络结构261采样时间1ms输入信号权值初值机取值
仿真m文件程序：
BP simulation
clear all
clear all

xite05
alfa05
w1rands(26) value of w1initially by random
w1_1w1w1_2w1

w2rands(61) value of w2initially by random
w2_1w2w2_2w2_1


dw10*w1

x[00]＇

u_10
y_10

I[000000]＇ input of yinhanceng cell
Iout[000000]＇ output of yinhanceng cell
FI[000000]＇

ts0001
for k111000

time(k)k*ts
u(k)05*sin(3*2*pi*k*ts)
y(k)(u_109*y_1)(1+y_1^2)

for j116
I(j)x＇*w1(j)
Iout(j)1(1+exp(I(j)))
end

yn(k)w2＇*Iout output of network
e(k)y(k)yn(k) error calculation
w2w2_1+(xite*e(k))*Iout+alfa*(w2_1w2_2) rectify of w2

for j116
FI(j)exp(I(j))(1+exp(I(j))^2)
end

for i112
for j116
dw1(ij)e(k)*xite*FI(j)*w2(j)*x(i) dw1 calculation
end
end
w1w1_1+dw1+alfa*(w1_1w1_2) rectify of w1

jacobian information
yu0
for j116
yuyu+w2(j)*w1(1j)*FI(j)
end
dyu(k)yu

x(1)u(k)
x(2)y(k)

w1_2w1_1w1_1w1
w2_2w2_1w2_1w2
u_1u(k)
y_1y(k)
end
figure(1)
plot(timey＇r＇timeyn＇b＇)
xlabel(＇times＇)ylabel(＇y and yn＇)
figure(2)
plot(timeyyn＇r＇)
xlabel(＇times＇)ylabel(＇error＇)
figure(3)
plot(timedyu)
xlabel(＇times＇)ylabel(＇dyu＇)

运行结果：




（二） 采RBF网络仿真
网路结构241采样时间1ms输入信号权值初值机取值高斯基函数初值
仿真m文件程序：
RBF simulation
clear all
clear all

alfa005
xite05
x[00]＇

b15*ones(41)
c05*ones(24)
wrands(41)

w_1ww_2w_1
c_1cc_2c_1
b_1bb_2b_1
d_w0*w
d_b0*b
y_10

ts0001
for k112000
time(k)k*ts
u(k)05*sin(1*2*pi*k*ts)
y(k)(u(k)09*y_1)(1+y_1^2)

x(1)u(k)
x(2)y_1

for j114
h(j)exp(norm(xc(j))^2(2*b(j)*b(j)))
end
ym(k)w＇*h＇
em(k)y(k)ym(k)

for j114
d_w(j)xite*em(k)*h(j)
d_b(j)xite*em(k)*w(j)*h(j)*(b(j)^3)*norm(xc(j))^2
for i112
d_c(ij)xite*em(k)*w(j)*h(j)*(x(i)c(ij))*(b(j)^2)
end
end

ww_1+d_w+alfa*(w_1w_1)
bb_1+d_b+alfa*(b_1b_2)
cc_1+d_c+alfa*(c_1c_2)

Jacobian information
yu0
for j114
yuyu+w(j)*h(j)*(c(1j)x(1))b(j)^2
end
dyu(k)yu

y_1y(k)

w_2w_1
w_1w
c_2c_1
c_1c
b_2b_1
b_1b

end
figure(1)
plot(timey＇r＇timeym＇b＇)
xlabel(＇time(s)＇)ylabel(＇y and ym＇)
figure(2)
plot(timeyym＇r＇)
xlabel(＇time(s)＇)ylabel(＇error＇)
figure(3)
plot(timedyu＇b＇)
xlabel(＇time(s)＇)ylabel(＇dyu＇)
仿真结果：

























































（三） DRNN网路仿真
采样时间1ms输入信号权值初值机取值取
仿真m文件程序：

DRNN simulation
clear all
clear all

wdrands(71)
worands(71)
wirands(37)

xjzeros(71)
xj_1xj

u_10y_10
xitei035
xited035
xiteo035

ts0001
for k111000
time(k)k*ts
u(k)05*sin(4*pi*k*ts)
y(k)(u_109*y_1)(1+y_1^2)

Ini[u_1y_11]＇
for j117
sj(j)Ini＇*wi(j)+wd(j)*xj(j)
end
for j117
xj(j)(1exp(sj(j)))(1+exp(sj(j)))
end

Pj0*xj
for j117
Pjwo(j)*(1+xj(j))*(1xj(j))*xj_1(j)
end

Qij0*wi
for j117
for i113
Qij(ij)wo(j)*(1+xj(j))*(1xj(j))*Ini(i)
end
end

ymk0
for j117
ymkymk+xj(j)*wo(j)
end

ym(k)ymk
e(k)y(k)ym(k)

wowo+xiteo*e(k)*xj
wdwd+xited*e(k)*Pj
wiwi+xitei*e(k)*Qij
xj_1xj
u_1u(k)
y_1y(k)
end
figure(1)
plot(timey＇r＇timeym＇b＇)
xlabel(＇times＇)ylabel(＇y and ym＇)
figure(2)
plot(timeyym＇r＇)
xlabel(＇times＇)ylabel(＇error＇)
仿真结果：


文档香网(httpswwwxiangdangnet)户传

《香当网》用户分享的内容，不代表《香当网》观点或立场，请自行判断内容的真实性和可靠性！
该内容是文档的文本内容，更好的格式请下载文档