test4_2_1.m

来自「两点间的头镇试验源码」· M 代码 · 共 13 行

clc,clear
t=linspace(0,pi,20000);
y1=sin(t)'*[1,-1];  % 包络线
y=sin(t).*sin(9*t);

% 求过零点
t0=t(find(abs(y)<1e-3));  % 用搜索法找零点，当然数值法找零点可用函数fzero，但需要分别给自变量初始点
% t0=pi*(0:9)/9;  % 用数学知识求得过零点的准确值，不实用
y0=sin(t0).*sin(9*t0);
plot(t,y1,'r:',t,y,'b',t0,y0,'bo','markersize',10)
axis([0,pi,-1,1])