Hi Debi
clear all
dx=.001;x=[-2*pi:dx:2*pi];
a0=5;
da=.1;
alow=a0-da;
ahigh=a0+da;
a=[alow:0.1*da:ahigh];
b0=3;
db=.2;
blow=b0-db;
bhigh=b0+db;
b=[blow:0.1*db:bhigh];
y2low=0;
y2high=0;
for k=1:1:numel(x)
y2=a*cos(x(k))+b*sin(x(k));
y2low = [y2low min(y2)];
y2high = [y2high max(y2)];
end
y=a0*cos(x)+b0*sin(x);
ny=[1:1:numel(y)];
plot(ny,y,'r');grid on
hold all;
plot(ny,y2low([2:end]),'b');
plot(ny,y2high([2:end]),'b');
.
the error function you want delY is
delY=y2high-y2low
if you find this answer useful would you please be so kind to mark my answer as Accepted Answer?
To any other reader, please if you find this answer
please click on the thumbs-up vote link
thanks in advance
John BG
