clc
clear all
close all
%% 参数
f=[0.0076,0.012,0.016]; %频率
FFi=[0 0.5 1 1.5]; %相位
A= 255 ; %背景光强；
B= 255 ; %振幅
[~,p]=size(f);%x频
[~,b]=size(FFi);%x步
%% 生成图片大小
X_s=460; %像素行数
Y_s=819; %像素列数
%% Function
for xx=1:p
for yy=1:b
BW(1:X_s,1:Y_s)=uint8(255);
for x=1:Y_s
I(x)=A+B*cos(2*pi*f(xx)*x+FFi(yy)*pi);
end
I=I/(A+B);
I_x=1:819;
gcf=figure (1);
plot(I_x,I);
for i=1:X_s
for j=1:Y_s
BW(i,j)=I(j)*BW(i,j);
end
end
% figure (2);
% imshow(BW);
% figure (3);
% plot(I_x,BW(2,:));
saveas(gcf,[ '频率' num2str(f(xx)) '相位' num2str(FFi(yy)) '.jpg'])
imwrite(BW,[ '频率' num2str(f(xx)) '相位' num2str(FFi(yy)) '.bmp']);
end
end
close 1;
