Question

matlab 哪里错了 matlab大神

>> N=1000;T=500;h=0.001;alpha=4;
>> v=2;vt=0.05;
>> x=rand(1,N);
>> xs=sort(x);
>> mox1=mean(xs);
>> y=x;
>> for k=1:T
z(k)=1/cos(y(k)).^2;
end;
>> for b=1 T
Fz(b)=tan(b);
end;

>> zs=sort(z);
>> moz=mean(zs);
>> mo=[num2str(moz),''];
>> mo=mo(1:6);
>> Dz=sum((zs-moz).^2)./T;
>> dsp=[num2str(Dz),''];
>> dsp=dsp(1:6);
>> k=0;T=-1;
>> L=sqrt(T).*max(abs(zs-k./T));
>> kolm=[num2str(L),''];
>> kolm=kolm(1:6);
>> disp('С В П М О Ф.');
С В П М О Ф.
>> disp('');
>> disp(['М О',mo,'']);
М О1.5889
>> disp(['Д',dsp,'']);
Д0.3939
>> disp(['В К К',kolm,'']);
В К К0+3.39
>> plot(zs,[0:T-1]/T,'-r',[0:tan(1):tan(1)*length(Fz)-1]/(length(Fz)),Fz,'-g');
错误使用 plot
矢量长度必须相同。

Answer 1

因为  >> k=0;T=-1;  使得T不再是500了，导致最后 [0:T-1]/T 生成矩阵时出错。

在最后画图前加T=500 或把plot中的T改成500;  例如

T=500; 

plot(zs,[0:T-1]/T,'-r',[0:tan(1):tan(1)*length(Fz)-1]/(length(Fz)),Fz,'-g');

结果：