RESONANCE — Master the Frequency - MATLAB One-Button Game Mini Hack

Roatis
/
RESONANCE — Master the Frequency
about 16 hours ago
4
23
8
1
9305
game.m
function g=game()
g.init=@initState;g.update=@updateState;g.render=@renderFrame;
end

function s=initState(ax)
ax.YDir='reverse';ax.Color=[.91 .90 .84];
R(ax,[0 0 800 9],[.95 .41 .22]);
T(ax,40,54,'RESONANCE',[.12 .23 .22],32,'left');
T(ax,43,82,'INTERACTIVE PHYSICS LAB',[.40 .46 .42],11,'left');
T(ax,752,57,'SPACE = IMPULSE',[.77 .32 .19],12,'right');
R(ax,[34 116 732 272],[.10 .19 .18]);
R(ax,[34 116 732 5],[.95 .49 .24]);
s.lt=T(ax,59,151,'',[.82 .89 .80],12,'left');
s.sp=T(ax,740,151,'',[.54 .70 .63],10,'right');
s.lz=R(ax,[310 195 60 143],[.20 .33 .29]);
s.rz=R(ax,[570 195 60 143],[.20 .33 .29]);
line(ax,[205 735],[271 271],'Color',[.43 .55 .48],...
'LineWidth',2,'HitTest','off','PickableParts','none');
line(ax,[470 470],[194 339],'Color',[.37 .54 .46],...
'LineStyle',':','HitTest','off','PickableParts','none');
R(ax,[146 228 15 82],[.95 .49 .24]);
s.coil=line(ax,[161 446],[271 271],'Color',[1 .68 .34],...
'LineWidth',2.6,'HitTest','off','PickableParts','none');
s.ball=R(ax,[452 253 36 36],[.98 .85 .56]);
s.ball.Curvature=[1 1];
s.ball.EdgeColor=[.96 .50 .25];
s.ball.LineWidth=2;
s.dot=R(ax,[463 264 8 8],[.18 .28 .25]);
s.dot.Curvature=[1 1];
s.st=T(ax,400,177,'',[1 .83 .57],17,'center');
s.exp=T(ax,400,363,'',[.75 .87 .77],11,'center');
R(ax,[34 412 480 140],[.14 .24 .22]);
T(ax,56,440,'OSCILLOSCOPE / DISPLACEMENT',...
[.76 .85 .73],11,'left');
for y=468:24:533
line(ax,[58 494],[y y],'Color',[.23 .35 .31],...
'HitTest','off','PickableParts','none');
end
line(ax,[58 494],[497 497],'Color',[.40 .55 .45],...
'HitTest','off','PickableParts','none');
s.wave=line(ax,linspace(58,494,130),497*ones(1,130),...
'Color',[1 .66 .34],'LineWidth',2,...
'HitTest','off','PickableParts','none');
R(ax,[528 412 238 140],[.14 .24 .22]);
T(ax,547,440,'AMPLITUDE',[.76 .85 .73],11,'left');
R(ax,[546 475 204 15],[.29 .39 .33]);
s.z=R(ax,[600 475 80 15],[.31 .65 .42]);
s.met=R(ax,[546 474 3 18],[1 .66 .34]);
s.read=T(ax,549,461,'',[1 .88 .67],12,'left');
s.hb=R(ax,[546 521 1 8],[.49 .84 .58]);
T(ax,548,515,'STABILITY',[.67 .79 .70],9,'left');
s.fb=T(ax,400,399,'',[.72 .35 .22],10,'center');
s.ft=T(ax,400,585,'',[.22 .34 .30],11,'center');

s.names={'ENERGY TRANSFER','DAMPING',...
'PRECISION CONTROL','NATURAL FREQUENCY',...
'RESONANCE CHALLENGE'};

s.less1={...
'A PUSH ALONG THE MOTION ADDS KINETIC ENERGY.',...
'DAMPING REMOVES ENERGY FROM THE OSCILLATOR.',...
'AN OPPOSING IMPULSE CAN REDUCE KINETIC ENERGY.',...
'NATURAL FREQUENCY DEPENDS ON MASS AND STIFFNESS.',...
'PHASE-ALIGNED IMPULSES BUILD UP OSCILLATIONS.'};

s.less2={...
'BETTER TIMING GIVES A LARGER AMPLITUDE.',...
'WELL-TIMED IMPULSES REPLACE THE ENERGY LOST.',...
'BRAKING HELPS BRING THE AMPLITUDE INTO TARGET.',...
'HIGHER FREQUENCY MEANS A SHORTER PERIOD.',...
'THIS IS A DISCRETE ANALOGUE OF RESONANT DRIVING.'};

s.omega=[5 5 5 7.2 5.8];
s.damp=[.12 .22 .08 .12 .15];
s.push=[165 165 120 130 145];
s.low=[90 65 75 90 105];
s.high=[175 150 140 150 180];
s.duration=[.9 1 1.05 1.15 1.2];

s.level=1;s.total=0;s.tp=0;s.fails=0;
s=reset(s);s.mode=0;
end

function s=reset(s)
s.x=0;s.v=0;s.time=0;s.hold=0;s.np=0;
s.delta=0;s.cd=0;s.A=0;s.mode=1;
s.nb=0;s.sync=0;s.hist=zeros(1,130);
s.w=s.omega(s.level);s.beta=s.damp(s.level);
s.j=s.push(s.level);
end

function s=updateState(s,dt,pressed)
if s.mode==0
if pressed,s.mode=1;end
return
elseif s.mode==2
if pressed
if s.level==5
s.mode=4;
else
s.level=s.level+1;
s=reset(s);
end
end
return
elseif s.mode==3
if pressed,s=reset(s);end
return
elseif s.mode==4
if pressed
s.level=1;s.total=0;s.tp=0;s.fails=0;
s=reset(s);s.mode=0;
end
return
end

s.time=s.time+dt;s.total=s.total+dt;
s.cd=max(0,s.cd-dt);

if pressed && s.cd<=0
s.delta=s.v*s.j+.5*s.j^2;

if s.delta<0
s.nb=s.nb+1;
end

if s.v>0 && abs(s.x)<40
s.sync=s.sync+1;
end

s.v=s.v+s.j;
s.np=s.np+1;
s.tp=s.tp+1;
s.cd=.5;
end

s.v=s.v+(-s.w^2*s.x-2*s.beta*s.v)*dt;
s.x=s.x+s.v*dt;
s.A=sqrt(s.x^2+(s.v/s.w)^2);
s.hist=[s.hist(2:end) s.x];

ok=s.A>=s.low(s.level) && s.A<=s.high(s.level);

if s.level==3
ok=ok && s.nb>=1;
elseif s.level>=4
ok=ok && s.sync>=s.level-2;
end

if ok
s.hold=s.hold+dt;
else
s.hold=0;
end

if s.hold>=s.duration(s.level)
s.mode=2;
elseif s.A>245
s.mode=3;
s.fails=s.fails+1;
end
end

function renderFrame(~,s)
p=470+s.x;

s.ball.Position=[p-18 253 36 36];
s.dot.Position=[p-7 264 8 8];
s.coil.XData=linspace(161,p-20,32);
s.coil.YData=271+[0 repmat([10 -10],1,15) 0];
s.wave.YData=497-.19*s.hist;

lo=s.low(s.level);
hi=s.high(s.level);

s.lz.Position=[470-hi 195 hi-lo 143];
s.rz.Position=[470+lo 195 hi-lo 143];
s.z.Position=[546+204*lo/245 475 204*(hi-lo)/245 15];
s.met.Position=[546+min(s.A,245)/245*204 474 3 18];
s.hb.Position=[546 521 max(1,204*s.hold/s.duration(s.level)) 8];

s.lt.String=sprintf('EXPERIMENT %d / 5 - %s',...
s.level,s.names{s.level});

s.sp.String=sprintf('f0 = %.2f Hz | beta = %.2f /s',...
s.w/(2*pi),s.beta);

s.read.String=sprintf('A = %.0f / %.0f - %.0f',...
s.A,lo,hi);

if s.mode==0
s.fb.String='';
s.st.String='PRESS SPACE TO START';
s.exp.String='CONTROL A SIMULATED MASS-SPRING SYSTEM';
s.ft.String='FIVE EXPERIMENTS / ONE BUTTON';

elseif s.mode==1
s.st.String='CONTROL THE OSCILLATION';

if s.np==0
s.fb.String='WATCH THE MOTION BEFORE YOUR FIRST PUSH';
elseif s.A<lo
s.fb.String='AMPLITUDE LOW: ADD ENERGY WITH A TIMED PUSH';
elseif s.A>hi
s.fb.String='AMPLITUDE HIGH: TIME A BRAKING IMPULSE';
elseif s.level==3 && s.nb==0
s.fb.String='IN TARGET: ONE BRAKING IMPULSE REQUIRED';
elseif s.level>=4 && s.sync<s.level-2
s.fb.String='IN TARGET: SYNCHRONIZE NEAR THE CENTER';
elseif s.delta>=0
s.fb.String='IN TARGET: LAST IMPULSE ADDED ENERGY';
else
s.fb.String='IN TARGET: LAST IMPULSE REMOVED ENERGY';
end

if s.level==1
msg='PUSH WITH THE MOTION TO ADD ENERGY';
elseif s.level==2
msg='COMPENSATE FOR STRONGER DAMPING';
elseif s.level==3
msg='BRAKE ONCE, THEN HOLD THE TARGET';
elseif s.level==4
msg='FASTER OSCILLATION: MAKE 2 SYNCHRONIZED PULSES';
else
msg='MAKE 3 PHASE-ALIGNED PULSES AND STABILIZE';
end

s.exp.String=msg;

s.ft.String=sprintf(...
'TIME %.1fs  PULSES %d  BRAKES %d  SYNC %d',...
s.time,s.np,s.nb,s.sync);

elseif s.mode==2
s.st.String=sprintf('LESSON %d / 5 - COMPLETE!',s.level);
s.exp.String=s.less1{s.level};
s.fb.String=s.less2{s.level};
s.ft.String='PRESS SPACE TO CONTINUE';

elseif s.mode==3
s.fb.String='TIP: PUSH WHILE MOVING LEFT FAST TO BRAKE';
s.st.String='EXCESSIVE OSCILLATION';
s.exp.String='AMPLITUDE EXCEEDED THE SAFETY LIMIT';
s.ft.String='SPACE TO RETRY THIS EXPERIMENT';

else
s.fb.String='';
s.st.String='ALL FIVE EXPERIMENTS COMPLETE!';
s.exp.String='YOU MASTERED ENERGY, DAMPING AND IMPULSE TIMING';
s.ft.String=sprintf(...
'TOTAL %.1fs / %d PULSES / %d FAILS - SPACE TO RESTART',...
s.total,s.tp,s.fails);
end
end

function h=R(ax,pos,c)
h=rectangle(ax,'Position',pos,'FaceColor',c,...
'EdgeColor','none','HitTest','off','PickableParts','none');
end

function h=T(ax,x,y,z,c,sz,al)
h=text(ax,x,y,z,'Color',c,'FontSize',sz,...
'FontName','Consolas','HorizontalAlignment',al,...
'HitTest','off','PickableParts','none');
end
Play Game
Preview

Play Game opens this in MATLAB Drive. Click Run this app there to play.