iddlZddlmZeddj ddZedddzj eZed dj ddZ ed ddzj eZ e dd Z e dd Z Gd d Z e ddddddZdZdZdZdZ ej&j)dddZeeZeeZej1eej3eej4deej7ee\ZZZeredz Zej?eeZej@ejBeejDeejFejHzz dzZ%edzdk(r-edkDrdndZ&erdndZ'e(deddeddedd ed!d"e&d#e' ed$kDred%krese(d&edd'edyedz Zed(kDr e(d)y@)*N)readz../X_train.waviz../y_train.wav z ../X_test.wavz ../y_test.wavicdeZdZ ddZdZdZdZdZdZdZ dZ d Z d Z dd Z d Zd Zy)CC_MLPClassifierctjj||dz|_tjd|f|_tjj||dz|_tjd|f|_tjj||dz|_tjd|f|_ ||_ ||_ ||_ g|_ g|_y)N{Gz?r)nprandomrandnW1zerosb1W2b2Mbcollapse_threshold energy_budget cycle_check M_history error_history)self input_size hidden_size output_sizerrrs c01.py__init__zCC_MLPClassifier.__init__s))//*k:TA((A{+,))//+{;dB((A{+,k:TA1k*+#5*&c.tjd|S)Nr)r maximumrxs rreluzCC_MLPClassifier.relu$szz!Qr c6tj|dkDddS)Nrr)r wherer#s rrelu_derivativez CC_MLPClassifier.relu_derivative'sxxAq!$$r ctj|tj|ddz }|tj|ddz S)NrTaxiskeepdims)r expmaxsum)rr$exp_xs rsoftmaxzCC_MLPClassifier.softmax*s:q266!!d;;<rvve!d;;;r cvtj||j|jz|_|j |j|_tj|j |j|jz|_ |j|j|_ |jSN) r dotrrz1r%a1rrz2r1output)rXs rforwardzCC_MLPClassifier.forward.sv&&DGG$tww.))DGG$&&$''*TWW4ll477+ {{r c|jd}tj|tj|dzz |z }|S)Nrg& .>)shaper r/log)ry_truey_predmlosss r compute_losszCC_MLPClassifier.compute_loss5s< LLOvv} 5566: r c|jd}||z }tj|jj||z }tj |dd|z }tj||j j}||j|jz} tj|j| |z } tj | dd|z } | | ||fS)NrTr*) r<r r4r6Tr/rr(r5) rr9r>r?r@dz2dW2db2da1dz1dW1db1s rbackwardzCC_MLPClassifier.backward:s LLOvoffTWWYY$q(ffSq4014ffS$''))$D((11ffQSS#"ffSq4014Cc!!r c|j|}|j|||\}}}}|xjd|zzc_|xjd|zzc_|xjd|zzc_|xj d|zzc_y)Nr )r:rLrrrr)rr9r>r?rJrKrFrGs rupdatezCC_MLPClassifier.updateEsra!]]1ff=S#s 4#: 4#: 4#: 4#:r cR|j|}tj|dS)Nrr+)r:r argmax)rr9 probabilitiess rpredictzCC_MLPClassifier.predictMs  Q yyQ//r cT|j|}tj||k(Sr3)rSr mean)rr9r>ys rscorezCC_MLPClassifier.scoreQs" LLOwwqF{##r crtj||j|jz}|j||z }tj tj j|d}|jj|t|j|jkDr|jjd|j}||jkDrG|sE|j|xx||zzcc<|xj||j dzz c_d}nd}|||fS)z Conditional Collapse of M based on entropy/error. CCT Principles: - Only update when error > collapse_threshold - Detect periodicity to avoid oscillation - Pay with energy only when necessary rrPrTFr r4rrrrUlinalgnormrappendlenrpop_detect_oscillationr) rr9yt energy_weight activationerrentropy is_periodiccollapse_occurreds rcollapse_manifestation_Mz)CC_MLPClassifier.collapse_manifestation_MWsVVAtww'$&&0 ffRj:%''"))..1.56 !!'* t!! "T%5%5 5    " "1 %..0  T,, ,[ FF2J-#- -J FFmchhAh&66 6F $ !& );66r ct|j|jkrytj|j}tj |}tj tj |}|dkr|dkDryy)z/Detect if M is oscillating (periodic behavior).F皙?r T)r]rrr arrayvardiff)rhistoryvariance diff_variances rr_z$CC_MLPClassifier._detect_oscillationysp t!! "T%5%5 5((4--.66'?rwww/0  c>md2r ctj||j|jz}|j||z }tj tj j|d}|jj|t|j|jkDr|jjd|j}|rd}n|dkrd}n||jkrd}nd}|dkDrB|j|xx||zzcc<|xj||j dzz c_|||fS) a Adaptive M update using CCT energy efficiency. Strategy: - If system is stable (low entropy): minimal update - If system is unstable (high entropy): full update - If system is oscillating: freeze M, let MLP adapt rrPrgrigMbP?r g?rY)rr9r`rbrcrdrelrs rupdate_M_adaptivez"CC_MLPClassifier.update_M_adaptivesVVAtww'$&&0 ffRj:%''"))..1.56 !!'* t!! "T%5%5 5    " "1 %..0  B s]B t.. .BB 8 FF2J"s( "J FFb3888++ +FK''r N)g?? )r )__name__ __module__ __qualname__rr%r(r1r:rBrLrNrSrWrgr_rrr rrr sM$'" , %< "0$ 7D"&(r rdrtg?rs)rrrrrri`r ACTIVE CONSERVED OSCILLATINGSTABLEzStep 4dz | Score: z.4fz | Entropy: z | LR: z.3fz | Energy: z | State: g\(\?riu ✓ COLLAPSE ACHIEVED: Score z , Entropy izMax iterations reached))numpyr scipy.io.wavfilerreshapeX_trainastypeinty_trainX_testy_testX20yt20rftotal_collapse_opstotal_energy_savedoscillation_freezesir randintidxr9r`r:rNeyerrrdlr_usedrerWrUrr4rr total_error energy_status osc_statusprintrxr rrsf!   # + +B 4  !! $q ( 0 0 5 o q ! ) )"c 2   "Q & . .s 3 Udm et}c(c(N  ))  Auc *C A BIIaLHHQr 2%&$7$72$>!GWkq  GGArNE"''133r7fbffQo&;