作者：Sreetam Bhaduri

%% Wankel Engine Housing Profile

clear;

clc;

fprintf('It''s an no-parametric program. So please provide all the inputs in same system of measurement\n');

%% User Inputs

e=input('Please provide the eccentricity magnitude\n');

fprintf('K Factor=Radius of Rotor/Eccentricity\nIt''s value ranges in [0,4].\n');

k_factor=input('What''s the K Factor of the engine?\n');

%% Solver

theta=[0:pi/1000:2*pi/3];

R=e*k_factor;

if k_factor<=4

   for n=1:1:3

       x=(e*cos(3*theta))+(R*cos(theta+(2*n*pi/3)));

       y=(e*sin(3*theta))+(R*sin(theta+(2*n*pi/3)));

       %% Plotting parameters

       plot(x,y,'LineWidth',2);

       title('Wankel Engine Housing Profile');

       xlabel('X Axis');

       ylabel('Y Axis');

       hold on;

   end

else

   disp('Learn the basics for Wankel Engine, then try again.');

end

%% End of MATLAB Program

下载地址：

http://page2.dfpan.com/fs/elc5j2f21329216ec07/  

更多精彩文章请关注微信号：