【求助】请大家帮我看看这个程序怎么只出一个波形
function Y1 = Lorenz_response(tspan);%%计算处于响应地位的Lorenz系统的数值解，并由此画出其相图yinit = ;% 初始化输入
y(1:3) = yinit;
tstart = 0; % 时间初始值
tstep = 1e-1; % 时间步长
wholetimes = 1e2; % 总的循环次数
steps = 1; % 每次演化的步数
iteratetimes = wholetimes/steps; % 演化的次数
S=output;
results=zeros(iteratetimes,3);
error=zeros(iteratetimes,3);
for i=1:50;
tspan = tstart:0.5*tstep:(tstart + tstep*steps);
= ode45('deng', tspan, y);
error(i,2:4)=Y1(3,:)-S(i+1,2:4);
error(i,1)=T(3)-S(i+1,1);
y = Y1(3,:);
y(1)=S(i+1,2);
tstart = tstart + tstep*steps; % 重新定义起始时刻
results(i,:)=Y1(3,:);
end
figure(1)
plot3(results(:,1),results(:,2),results(:,3));
error;
function s=output;
L=5;
n=0:(L-1);
a=1;
bit=5;
range=2^3;
tstart = 0; % 时间初始值
tstep = 1e-1; % 时间步长
wholetimes = 1e2; % 总的循环次数
tspan=tstart:tstep:wholetimes*tstep
= ode45('deng',tspan,);
s=;
interval=(a*2)/(range-1);
partition = [-a:interval:a];
interval2 = (a*2)/(range)
codebook=[-a:interval2:a];
= quantiz(s,partition,codebook);
%plot(n,s,'x',n,quantized,'.');
%legend('Original signal','Quantized signal');
s=(index-16)/4;
figure(2)
plot3(Y(:,1),Y(:,2),Y(:,3));
这个是经过AD与DA转换之后的程序！
错误：??? Index exceeds matrix dimensions.
Error in ==> Lorenz_response at 15
error(i,2:4)=Y1(3,:)-S(i+2,2:4); 错误信息不是很明显么，你的矩阵长度不相等，你看一下Y1是多少列的。
页:
[1]