|
原帖由 wwbeyondww 于 2009-1-4 13:40 发表
想了一下,我是找不到好办法了,搞了个循环,
y=[-0.5572 -3.3084 -1.6488 -1.3941 -0.0379 -0.2872 -0.7380......];%大小为256 的向量(在附件中)
x=linspace(0,20,256); ...
怎麽没人发现异样! 楼上提供的答案是对的, 但所给的程序是错的!
我想可能复制错了吧!
不好意思! 帮忙改了下!
load y.txt; x=linspace(0,20,256); x=x'; zero_1=(zeros(1,length(x))); plot(x,y,'b',x,zero_1,'r');
a=y(1:end-1).*y(2:end); b=find(a<=0); xx=zeros(size(b));
c=[y(b),y(b+1)]; d=[x(b),x(b+1)];
for k=1:length(b), xx(k)=interp1([c(k,1), c(k,2)],[d(k,1),d(k,2)],0); end
[ 本帖最后由 ChaChing 于 2009-1-5 08:48 编辑 ] |
|