用matlab写一段程序：用bp整定pid的一个闭环控制系统，该闭环控制系统输入为阶跃响应，后经过bp-pid，再经过传递函数10/(s^3+2*s^2+s)后输出，并可以画出pid三个参数迭代变化曲线和输出响应曲线，并保证网络可以迭代多次，且网络输入向量到网络层的权值向量有3维

这是一个比较复杂的问题，需要一定的控制理论基础和MATLAB编程经验。以下是一个参考程序：

% 定义阶跃信号
t = 0:0.01:10;
u = ones(size(t));

% 定义传递函数
G = tf([10], [1, 2, 1, 0]);

% 定义PID控制器
Kp = 1;
Ki = 0.1;
Kd = 0.01;
pid = pid(Kp, Ki, Kd);

% 定义反馈系统
sys = feedback(pid*G, 1);

% 计算系统响应
[y, t] = step(sys, t);

% 绘制系统响应图像
plot(t, y);

% 定义训练数据
x = u;
d = y;

% 定义BP神经网络
net = newff(x, d, [3, 1], {'tansig', 'purelin'});

% 设置训练参数
net.trainParam.show = 50;
net.trainParam.lr = 0.05;
net.trainParam.epochs = 500;

% 开始训练
net = train(net, x, d);

% 获取网络输出
y_nn = sim(net, x);

% 绘制网络输出曲线
hold on;
plot(t, y_nn, 'r--');

% 绘制PID参数变化曲线
figure;
plot(1:net.trainParam.epochs, net.IW{1});
hold on;
plot(1:net.trainParam.epochs, net.LW{2, 1});
plot(1:net.trainParam.epochs, net.b{1});
plot(1:net.trainParam.epochs, net.b{2});

% 显示结果
disp(['Kp: ', num2str(Kp)]);
disp(['Ki: ', num2str(Ki)]);
disp(['Kd: ', num2str(Kd)]);

需要注意的是，这个程序只是一个参考，具体的实现会因为具体的控制问题而有所不同。