MATLAB画李萨如图形

%% 画李萨如图

% 由于是按一定步长取散点作图
% x,y方向频率不宜取得过大，若实际频率为φx10000Hz φy15000Hz
% 则输入a=1,b=1.5或2，3亦或4，6等较小的等比例数值

a = input('请输入x方向信号频率:');
b = input('请输入y方向信号频率:');
phi = input('请输入相位差φy-φx（以°为单位）:');
t = 0:0.0001:2*pi; 
x = sin(2*pi*a.*t);
y = sin(2*pi*b.*t+phi*pi/180);
plot(x,y);

输入

>> parameter_equation
请输入x方向信号频率:1
请输入y方向信号频率:1.8
请输入相位差φy-φx（以°为单位）:0