From e15de78c005e16dfdcd26a0e85a3630a85de50ff Mon Sep 17 00:00:00 2001 From: Nicolas D Mower Date: Thu, 13 Apr 2017 18:41:26 -0400 Subject: [PATCH] contents --- boussinesq_lookup.m | 50 ++++++++++++++++++++++++++++++++++++++++++++ cost_logistic.m | 21 +++++++++++++++++++ initializeData.m | 26 +++++++++++++++++++++++ least_squares.m | 8 +++++++ plot1.png | Bin 0 -> 20567 bytes problem_1_data.m | 23 ++++++++++++++++++++ setdefaults.m | 3 +++ 7 files changed, 131 insertions(+) create mode 100644 boussinesq_lookup.m create mode 100644 cost_logistic.m create mode 100644 initializeData.m create mode 100644 least_squares.m create mode 100644 plot1.png create mode 100644 problem_1_data.m create mode 100644 setdefaults.m diff --git a/boussinesq_lookup.m b/boussinesq_lookup.m new file mode 100644 index 0000000..5fe7f64 --- /dev/null +++ b/boussinesq_lookup.m @@ -0,0 +1,50 @@ +function [sigma_z]=boussinesq_lookup(q,a,b,z) +fmn= [0.1,0.02926,0.03007,0.03058 + 0.2,0.05733,0.05894,0.05994 + 0.3,0.08323,0.08561,0.08709 + 0.4,0.10631,0.10941,0.11135 + 0.5,0.12626,0.13003,0.13241 + 0.6,0.14309,0.14749,0.15027 + 0.7,0.15703,0.16199,0.16515 + 0.8,0.16843,0.17389,0.17739]; + +m = a/z; +n = b/z; + +if (n>1.5) + n=1.6; +elseif (1.3<=n) + if (n<=1.5) + n=1.4; + end +elseif 1.3<=n + n=1.2; +end + +c = ones(1,4); +d = ones(1,4); + +for i=1:4 +[~,p]=min(abs(m-fmn(:,1))); +M =fmn(p,1); +fmn(p,1)=0; +if n==1.2 + t=fmn(p,2); +elseif n==1.4 + t=fmn(p,3); +elseif n==1.6 + t=fmn(p,4); +end + +d(i) = M; +c(i) = t; +end + +b1=c(1); +b2=(c(2)-c(1))/(d(2)-d(1)); +b3=(((c(3)-c(2))/(d(3)-d(2)))-((c(2)-c(1))/(d(2)-d(1)))/(d(3)-d(1))); +b4=(((c(4)-c(3))/(d(4)-d(3)))-((c(3)-c(2))/(d(3)-d(2)))-((c(2)-c(1))/(d(2)-d(1))))/(d(4)-d(1)); +f3=b1+(b2*(m-d(1)))+(b3*(m-d(1))*(m-d(2)))+(b4*(m-d(1))*(m-d(2))*(m-d(3))); +sigma_z=q*f3; + +end \ No newline at end of file diff --git a/cost_logistic.m b/cost_logistic.m new file mode 100644 index 0000000..0932489 --- /dev/null +++ b/cost_logistic.m @@ -0,0 +1,21 @@ +function [cost, gradient] = cost_logistic(a, x, y) +cost = 0; +gradient = 0; +t = a(1)+a(2).*x; +sigma = 1./(1+exp(-t)); +cost = sum(-y.*log(sigma)- (1-y).*log(1-sigma)); +costFun = @ (a) sum(-y.*log((1./(1+exp(-(a(1)+a(2).*x)))))-(1-y).*log(1-(1./(1+exp(-(a(1)+a(2).*x)))))); +gradient = (1/length(x))*sum((sigma-y).*t); +ai = [0 0]; +% Set options for fminunc +options = optimset('GradObj','on','MaxIter',400); +% Run fminunc to obtain the optimal theta +% This function will return theta and the cost +[theta, cost] = fminunc(costFun, ai); +t = theta(1)+theta(2).*x; +sigma = 1./(1+exp(-t)); +plot(x,y,'xb', x, sigma); +title('Regression') +xlabel('Temp (Degrees F)') +ylabel('Pass or Fail (1 or 0)') +end \ No newline at end of file diff --git a/initializeData.m b/initializeData.m new file mode 100644 index 0000000..01b2f4f --- /dev/null +++ b/initializeData.m @@ -0,0 +1,26 @@ +c = [1.0000 53.0000 1.0000; + 2.0000 57.0000 1.0000; + 3.0000 58.0000 1.0000; + 4.0000 63.0000 1.0000; + 5.0000 66.0000 0; + 6.0000 66.8000 0; + 7.0000 67.0000 0; + 8.0000 67.2000 0; + 9.0000 68.0000 0; + 10.0000 69.0000 0; + 11.0000 69.8000 1.0000; + 12.0000 69.8000 0; + 13.0000 70.2000 1.0000; + 14.0000 70.2000 0; + 15.0000 72.0000 0; + 16.0000 73.0000 0; + 17.0000 75.0000 0; + 18.0000 75.0000 1.0000; + 19.0000 75.8000 0; + 20.0000 76.2000 0; + 21.0000 78.0000 0; + 22.0000 79.0000 0; + 23.0000 81.0000 0]; +a = c(:,1); +x = c(:,2); +y = c(:,3); \ No newline at end of file diff --git a/least_squares.m b/least_squares.m new file mode 100644 index 0000000..57e99ab --- /dev/null +++ b/least_squares.m @@ -0,0 +1,8 @@ +function [a,fx,r2] = least_squares(Z,y) +a = Z\y; +Sr = sum((y-Z*a).^2); +r2 = 1-Sr/sum((y-mean(y)).^2); +x = Z(:,2); +fx = a(1)+(a(2)*x)+(a(3)*x.^2); +end + diff --git a/plot1.png b/plot1.png new file mode 100644 index 0000000000000000000000000000000000000000..153bba4cbff763b1cff02cf68d38c48f3f907c37 GIT binary patch literal 20567 zcma&O2RxSl`!;?V5oOO3vUik~E$-}( z&*%I6p8xOld!E;Gze?e{uKT**<2;Y^IF9p*)KHTrz`u--AP9k?0!kA>F!>M!qu@Ly z{3ItjVhjGkb5XeKh9HEU=)V}rT!b|6LtJ-7)!Vr9*aVcc*!$7;#_%I%_uKmJTF%z) zUgoY=i0UnEm0L30SG8@euBw{5+v+*H+PiUIRdl*{O^%r+r+0Bem*yxur&U6OD8L_X*OB>OTKh#eoG<)ZdqSSN=hQ( zPIo_|<}l!2Vg?2bBuM%y`s>yI``2eHVhGauBr2-yWZ!Vnm5hwc&CSiwaO%Og?zD`I zRY^Yo`O5|B*^*?eDuZ4(5y#LFY#b&rG2ge3&vA0<-@WTN`^I~s+Di0-Vp}85C1Iy4 zq@>gsM%7r~7$~FW{10Z{g%hoc5hJz)yK(JTGbMaee0dUkD~+4$JyvzA^i@#t@$t@Y zw!=dWz8sMJ6CLq;_hMu!;$(xR1$4 zywk$CsQ&HSx7X_QQUTm&L)}kl*4NixtII#VoSc&K;?+@;$7*>+#obVwEbkxVRrek| zcrg0+Wo~X>uBM`Lkt;GX@@|baL79WIv#EhWU!Xi35{yHTa2B`)lQ_}jVb4|C_5OX$ zLsr#vX49;U3^A)da(;gP`dkI(or43t2-NAa(@fLeLKLs782o^jon2!)9I<6nP%z`u zW-BOEPOXk@a{tL_DtZ5Nfpy>$7?_wce*fR@i}>H($J+2xa_hI3R;}Jc2{=vPtW3|y zpgE(Z74s~g5JlBmbm6O(Js`LNzd4e*a3RXNVW!50NRLNID23wp;o26a!SEXOL1%LhkCJ?C2e@$m4B27Z(jWm(X2 z9Mzm9uqcy~l9Kzd_2ETEz)g~plMC1l7eA#D-kH9mrKNRtdJ@_5;lmn?9875}z0~RV zM0I_^c!k|Rn{D>msv-;&6Y5UVgWx zQPj%HYMSDWHJPT23?^5SQlYrl@a|G$YFgS|A0JV@_377c+uPeash7NUwbT@%D0wVc z1|$xQR81!(bDlLH)0Tcf6nY8H#hgJ&uT?d zQWC$-r)TM>UK`8#^07O=eod_AFDZ#v&H9Hl{@0>2&0fwwNw1LX ztxfz|Of|KBsj8y#a3dl@e>OpZSyEE+)2B}i3=FW`>4Y4{%U-tH-MjbBv8N`zy0P)* z{=&(@Tx#}cRlteQ{sIi1v#o91T<7fU?95E2ub9hM>HUTM^(jtI>_IvhpQ|TJQfD0O z?Bfm2^fYJ1FJBJy_Qp$9kq#wu={1hQ3aUT;VtV#lPxs-$hH;vTuBqw3tmNLX)=maj zY`M#umKMPo65pzO?+N1-I!&l3C}>69Z8xT04~YL#k(LfRPTeYhkyBMwRa=|XvN85D z`&j;67{ND_v%_z(JpL%uw06Y%So(%qrXbJ95Kc<$$0}jtuNgYGP|(v4e46?B^Ct!d z#>mKs5p9~v=vqB1D=X#SSKHZXvtMR+x7?8(kE0q2zoz(TEg~5k2glmV>SLx@s@Nt! z2L}fX^QYe4q2hFWeEuPwucAqaPk0D+*m@|HiGH&uk8fs1#-EMmO;`E{E!go2+J=UP zMc-#zd|$kHF|ALMD)h{-!KtoRD^&cgjl7W+5m9U~F45@72n!3#Z8%#nx z0|mX|1S7W1Pnlvn^j0=D{SErc@6Y9!!nWKRtuz+1>>=L&T{9r5tfFH15+aiC?(!{^ zX31_pZ`=2kz~4WMU@pZx*UVDe8msbn4C`HJ%bi<~!e|cfL@mX@Sk^d9RNw#dLN!Z# zqO;R-y72)C2?;u&mg^@ah)l^>7AFcav zr12}WTN)a~2adzj+|6IVnuQ{VYsJOI5q~Ci%5^MEHJ9)Gn;Xnshill8mj7UctqMkIS^$GE(={4moCEc=T4o|K1ve(%0%^ z<#*rVcUf#lkT9j6a0prTC0RFE`RuQiEOhHHhIG2g`|hvZuO}q&_V#|@on%vI?HsLB zv8yo*Tz}eyn6)e?JqWvU2}s>z|Douen7hDi%h2`&EzC(M%CH zxqvCSyBO;09UnfF!-xL0G|=X=X9p2$qR5Z+{MPi)LRS>T`E1_}kt4^7^74+34oyu> zZXO=N($Vs+C@Rw?_edHyi!cJ3l%%BEhuao)xR^Qe`bHF@*l8*$iHU9y${+{GnTrgI zd2NXQ`oQ1E#P#dP$Ee6ihPexv1g*XlAdrhW}b&5e&gOG)XrRy$Z8`nECC{45P82MZ@H4&GKvEAd#I zJ*fWSx2M7}S48rRX5n_e6sdpv`4wq0kHfn^-w{ecKrl2k^qUYnTHe)l2X;A~q~9Tg zl&Q%{bPRwPIQ{tCg~k#BWN`K(FN6Q?ki5ffe`aO`xo4s9rM0ZAY|eX$R%lzX&B7w$ zc6@@251zau;u9N36nwtggPc&cN1pThc+R__ZxR}pkAAPd^!N9N4KnAuDWRt~HaMu^ zf1C8u1Zw#E_uV-vl2fp+s$&NYkt@@aqT)cgSr-v_*upTqg&yZfaV`%@tY(w_= zd7Ys*+v5A#mx;4(JFG$95I<^ET*a5qxp-TD&%o(3)6FuiwT;in$dD7p!9;>EuD?+v zKna+&Ba(@iB@e%19gPg2!AEk7>1$G59YLMe&h8&;ugpJ6k>p4e6Z~&C!9n~f{_?gd z#?9$3>T=Q44-U7?PjK?XTAOVv-u_C+gdrI9%6vZkfrXHbZN_*}gC@w6$stzAfmE8m zKHHi|^_je#eqgBQ@1pVC^KAqV&Z+R@5(#PPYooflWbpp;K58}LFQl)^+jRz>s!`GXn;ZCMc!2TFsY_|j1%YMG^W+yy+lv(fAFM}wa~Z+y&}kq?r$ zlMW2s|5J3n!~8?yWt@xDE~}$FGa=Z>RhJGa7K@8`2!f4)_g`P%US$(Je|&=-L25C5 z6afC7SNMPW>Zqx0QnX2oAko;Qni(IWK3*8lDBOh0L}QR@u2!hCt8f&gAp{8Hj3mREv!aQYt7^ELqxfoTU^N2DRL!*2hF_&%tr7)io+6 z;sT=Ms=3dN-x_9X6sHy=$VDo6uMnRDPNQ|2l^Nw+@f%gddQdqCJPsMwRr5?4dh678 zKaxgYiKBQ>&jCSju<;bd^bM3u))p;`ZoVRaB}{`pd{{>2GixB#)c4!656%%GITUem zwJPfr;sq>7kTnBC-GysB9HwvYX}#N{phkks;3rIMV*MQ>t?|5KTuJvH;39~Xg2MVW zi#Pn0Uonv4(2#-)(-h)?sZKkKMg|=Vgb1Rb(Ddc=0^#jN5B+HX^mvM6uNtus>6*B> zT_g$sJG-pRXBzGe7;hpx#R2s1kl>JlK=i>uad9EC>kEX>E)yc%1bB-2+T0t&#L1zR zqzLtOW@ZNT^#?H+81Tp`y^x?796ZHA^!e^r6cljK*OrzK3bCZtRws~G!$g?i{|W{M z(- zM2NKERfbDK-P_x8xtq{)3EGc7rCAAl6VD~Cn+!kQeliQ1^@xTj!x9QI5U$rn->o-L{!wT z-QDLmtomAf4`S#fCWnXJA<~?l{3#pfWM*cTkeKQ3w_h19A)=QQ_gvElq+667aG98R z@G~nG&c#NLRm!t`ekJ@1!t5{C5QH{lOFhz$Ow-Tr_&t4~k%2*Jo?0xOgzifd3G)tY z3XkfVn&9AIe?LD$I`N+54_OjZ@7|$veqjZ~k2yf4ck<46{u&{c1$R+$E_oW!RFR9! zRu*1d8o)%BQ%ggbla+!WKX!F+9WidSGJ zjkJVI_t^@mp?(C7iiH;|tmXE`(9qBavdkR=ilsgZ75#Eo7LS}AFhgYy zBDms}3VHmKwK*kVO{%G>0bGJtcHpTuvzJ{Gr!C= zZ^7?)`@MWDo&S2>7%3T<*T(cE3JQ<;4j7esv$v1)^71^NjIz&vVn&F$Tk%S&t9O>N z1GBE*>xT7$iG_X@sif_k{@GkwS{fT2-8m+`aNz1B&E4)`kT)!XvTB2|gv4}!mzZR5g|73kYK}pG(RMXkbiVBh=-%x=J z6y3}IHd7JFF)bxcKrxz6EVn;EH5v=F%=PEjI%cVnx7;sEscOiXbmV3_}IK{$2auCD<=7T zc5Ugx@t(NYkRs=wmKFxz;OB=@`ege(y?WY>+v`-SX*U*se9X|RH2OJNSGJN4bWQhS ze=1=8E7a8bCMK=--p73Y{Q2xp+u6O;uMrB2@fYdDpM7U_gt;4dc3qY+{zJCZ+?Ov} zjPcKJK49eGQ4PRjwCXTrlWlwTfmS-45dutUh=9$f=h4y8*#W;J2TDtz#DLOASxpUM zO~XpMj*bqEkV8altc^3UFF-eNNp%^CFjS2}O*Anvu~t!1^7FZfTMDmfrNvjJ!-|() z9Oo!7d~04eoL-yZdxmXR9ulae?a{kYH)hh`(^FjbO_1ksZK9^r8TG!q+Ztj8{IM^2 z3zQLO$s-M3Az1Hshe|!f(uwJK#@~ZQ}#RpgLuM?l~Al;MMO^hun4usB?xuGJglVC zskzfppbL)9wL#?uL6DY~R@iy=7>djDJT#Bw#?43%gl=zKV?mIa!8MW}(UXNs*fNJw zxY(h_<==4P#tX&arp*LkSR7O}uV24T5+H=BicW017rf7%zVH2X!+ZN|wBM|M`sZ~z zWLYCF?kgU(3MS%zO9t~cddKf#;^C41Jpulc{^%_;GAKbnswpkqF1GpROF(W=4m!cQ zM?1@umXA`;whVdxWCTJd^e8a_O`4&@NgW9yfR&*P)F zuZ@7d=rW8M#9PsT<3`x2U(}^$Ub`)+Ws3Id@m;@uXZu)phU`@IQkaFw8Uds&00ZLu zc-KForF+KB5`MibleCU&1Eenf<@4bnU zwoqs`I*u-*gQu<3OF!E^yo!QZ4%*MTCBv9UW%4k$5P!4f zfR02*PtTr(oA~BMR*az0ZB>OIKYpOx_f|$yIQ>rc#vmR?X0V81|NZU#S1Fmv$vSFk zG1{Dt^)l}6wE!(LGc)%sux`0k*Ce$G*B8HuX?CLo;4JXW3RbN+0qub1Ybz1WsA$w;{XU3aTA~Pt)jPK#fu3g0-#N%#YN!|#s%x^4}E;HZa$D%-L4DW zDoUcpT~bZ`jvI%s_s1el48b5h<+xKrbB~Blyz%}QR^WLUVJ%ur^h{3rEcPYK^vmV& z?*uENLSPNrX6ja2=@}b;)IF_#^~%hc1x_Bo$j{{qD&Ttq0pKbGsc!p;>%ba&xS(Ie z4nmm59cH9RtuYrFY0zUd3eC+2XpD}dw%B-hwKy>3=mI52@bARJHmSu!|Hux0WQ88( zQP@lC@6eaPBqg<_hXs!=C31|>8%F?m_VwG|d>n;2GO^o8P&)&|Eip}Mm2@g!S<`qV zXe+Js&;j!K-Mg1B%6qL0a{#*g$na>}ELd?GA*E62xZ#Bia_iugz^dh6Y(y<@f1(YeZ3zchLcOL=Fq&b}#9(J-WK)aU)9B`st* z6L!c7o1jX%01%DEm6cB4x}S=kKOapBn*ynO__2Db@j?FU*GB*up`#89DhdAX>muKbDqo&z&o2jAMr_MT-U)ySHIW8M76v zt*rr3Mo*s&HH0e*B&7 z54XQ*mbe40%1uKO9~#;nxQBa2fCWl=E3HG5iDN?MY06Z^<9f1E2wzi%D3IX7wt@0#Mff@1{{>6*CTU&?6 zzXs5dVnolBV>rIo5;i1=GvSfFWPhw0rt6c?#FV-D{MP)t$YtQQ4PLu0PC2)=0Qy*)_LrfmWE zixAn?zF6Ngqgn*+53=yHgoMw9%2a78kVe$A#9#X!IE^;jp+U5#B~S*vQ7xGJ)-?@K zsl$lqw_->)19RBOKsvv~EQr=TJlB1ow5+bKUY}|(ZFI@~n~0(#Dve$ejVlU%fi{r6H8kDg(SmuxKakhR6x40pp}%B)j>MT$WTUQWo5l=JM}r( zn9=#ppjM=ona*RQhnY-9@y;p0&cR{2+Ui2T97@2vBLs>WkAWTGumNzUgw61oXo{ zza_VSr9I_bTv}>iBqcsC^eraC>bK9PY~mFP-YA5WlyKDfN0tJpv7? zjTGUfshN3!O;43*6_F0HW&3p`_gGdBNS!^8U{%j54}O(vTs1*3Z8*(# zVetBW@vvbzkSdtuqHh8=ym#-OnVFfl_fZV3Sf|DI%1GI3QPRi|7B;8P`@28Cv_Z`T zsJr^b3+u^R;C#7(Ef5nEgWy;gKtWCU!3jYmxGN2Gbn49Cg&iI~+#j8t4S>Y~((6}F z1!lCA`6bbdWMpj2d9L*Y-~m zrqxoCiz{&K?(HS;QQfakgOrmdj_r7!lUo%DymL`9RIfOmPPgBiayd*70KQb(;^l$r< zktISMXDcGTTug2)+-aEv((t15fL8$CyE4vbg*$3ZG=BWwcKf()N5nFkoFzS3EzeNUF~4$==f2|!O^kV z{bz@cew9V2XE2soUhvH8F$d3Jx?qeT&%9eQu&tT$Kv?|n;e%{AA<0l0zYPSKN|P2* zyJ2=le0+SnBKF*1twtBB#5<4Kv=C}~sIlL?c^)x6Gvm9ravKGkgDnrOJ;o;|k37P{ zU0Yw5m64g`QCnVy@orxRgm)t?I^u;|X5Z}O!~{wFv9VHrD(^slf1*<1WDIf8c^C;% z1YrB;&j}F;iHWaCX=pN@JUPoGfjkufXj`Rh9^;q@z2&+TNQ2_%GvXi;+BrFyXlr*S zi-U>)h_!qPg=#zf&2=&nk7>*HYz8sw6nHlLb!o{9AZ}b-TtWgVW4ffjD2d}l^=;Nf zz@Rl8o@<(P(qeZST~akQ9>Mg;$<`DG!^bMr{0UomAXAJjv3GhpOJD!P)Rdt{ng4~r z%!{_6c{@!(#oX73j+))yefxtP5~G1(@?~%fq;TKe0aAw~9s&oH&yEfb&yub5o~B0^rjsPplvVAjj>@AB=z$ zhR;Xy=z+4yZd`kezIo|ybf5QD$4-uqbt?=; z=!TFjs_kn2gC>AxCJtHX>8`i`$kL1p3G?BW^lC0*i}Low5>*85m1r3Yi*jY=0xVW) z0#M~0&~GLH((Wsc0jI%xe1D;**%oXI?6s~?F`$P*X(d80;NagO+Io6=XUw`LbI(0v zD095gr(&y&j6VGJQ8@VmuW9Sr{CutkJ89(3jtlSkW%-@F+)o+Ts>Wb(U%$k%b59w$ z>epzAE7t)8?&-Ea)KGiW7m?b>>0LVZ_V$9mo1kQC;GuIQJym*~V+>~*u|-Cp14=9p zhC5nEDE-OcF&gLe7{@lHV=W8_M~HoVSJ9X3n(Ly3o4iO4HWCUpG3C%Y8jBlByVOL; z7HGKLIRE(_u)ba!f)K6HoUmC4dewsI<3|CS8?`STB#GQB2w$q^pSWoC-x zF>?kN57egwgOWo2byVPUd$thksRz#f`0hM587P1woSGrgy;?=qiRI~u|5E3`M7 z!@9t@(>5t=Gg7Jp6`cFShle1+1HO9{0QH9t9~v1wFDhDs>e=ld zG>U_ZoBaDQ>0e%Z7z|NqvL&{hrle4t^4k6$uoPhRDpFybAVJiUpf*cAJ4Q?7>3!xH@LT}5oe7JTHcs|5r=QhEp|EsXs=v} zScV*1Ds1}NwJ9yrAH_4tXGsPArFNsb009Cu{g1Z*=pLN=ah;av^u9ex5!_rwk#=lDJV#^VOed~(=k_QKNgHNeeo2o!JX z4}ZDJ0E#~lZ|_k80mp{`VbHS*jJ*<6AmZLPNY>}tOxGQ2)_&IF&oEV@IFsvZ&7ro1 z=?9ui1OWCa2H8-BRBx=ot?ef>+(KsbFskCXvog|Gf1 zqPG(^KL*5KJA3>7zCI%(Be|F>CC5K=ZUy`K`Ehb`!hX>9iwq4tM?yhG#g(MYOv+eY zTkBv=lLw2GlAK&6OI!pd0!B|za1V++)GkzA!^7ID&ZVJ*YWLT}{sOX{!%x*^6%{Q1 zMHPzAf|Are^^vV;eMp$(fV8xR=m6ph&|ZNjtEPZl|Amd6{r$?4QW&{qiMoX=e)D5)-e?-=&%c^y(VEYG{~y?N*eXo!wo9 zReK8!y66COy~`^@e8wVYXgcS9ekp2joYK9s-4^a_U}(r|)VSpw_$Q?s+k)_}%5)%ke-N!frtALP1%LO9hUx z#0cP$fhYcOf$G}s`qTl06sRmSG0zPW7&T%IQUMTaYHETC02VxyYAEJ-h4En1_c65n zoFU=iWOram|LMfb5o>h8&7lH?&)Arj)U6=~9J{@XG(Wm#+qqqW>IZK!M0PFgut~riK*|Jj zBSdeCAgKW15Wx;kHh%tJP$r|? ze|$)TNK&TvvO4lBlql=uEJ1O%ahrR2Imn@Ou-ZTj9y!YZW5UJ=8_09>4y>`rWpGNN z7naA@b_`$NKLFG{qsXKOzo(Kja=YOUaspaXVP(pL`nkFR=ob_@)K@zCt$yAnCUm-g z!eSKXhkoe)AO!Z;u>YEfNE#OqNF~LGb5A(Hm5atOFGjYIi(xyzIDqrbAt4y8<`)*` z<>XvY%JNwOsqn_A8)F`j6v?M+3K;GxAkmPgyhhWgfv`Zge*J0=JOfjdf}$cGiE6fF zW<VLS>L?IHlC@cw(^R@g%MF{{e*%yg zul6`SJw@$;9kSkkjx@4>;(!Op*K%Bpu)GP1D{+ZPdLbtRTJlB>H#pCM53Yt3;sbDy zBoL}f^tqk(O!D*c07>Bdtx>*$)p@2d_cUxD3R3V=T)cFna|ZqL2QvVe$|?v!)JbBS zVPL}$P}Y(0IX=8jLU84#M(!J+MN!&!?@lxshlPi4j#s@K8oJKNXaPhPyxzWj3yGam zh=Z3mThwC(mKMAtBNNlRjt&)7)w834Ksc%48V+Ktss9$dDmiyc=?2;26(ZOL1hjLX zLoOB&sI>?LsSrXtIB5a4Te1Khogv}750GXt)lTOX1RmHOK=iu$*Fhx+$|cjaYfiwq z%7$IMf14R;>J(L+F2G4iba}@eB>t4%Tn)v_S~2lN#AAh=1g)U6!1BG<9o_rsQ>*_U zPh~?EX69!Z8D}S(flErFUK<~7ED}i_53ZOYXc|*BH#>yJ8smGNWL>S={~}*iBrZ$% z2$1~sCRB`qpTer2G8AQNXQveC`3<><2|O{rD7rRZgQ@`nhlT}S1dQ7L;t7mnTagHOKGr(gPmnd!HxjxyI ziL*x=tUPtZa4@k?$od3Cb#FFhW(+{#z==vRs_4FyGVLc+r4COrj3SLGk5V(08?P`_9S51+;Fz3Aa(B#{31-hx4Hj%vQF} zBeiJDl_CM;M78=9ka6~?DS~2tiZYN4flU`6@)o$v1CAC`QI^imZ6K;jNwuY{*jQUv zm-v626)Bs!a}Dt~Cn}8_P>hcZjqyy`;onB04cR=M?1b3B+c&Hyj)(?9xv0@6!T}z< z)tN^mvT*#h#pfY{`fnqkd*!@=vIw56l9?Zz9ul+29BQLumX1YD4;R7uYH3sA7tZ}c zBq$!Kz64Z^i9CIB-D zNElyMX_@oUEqGTx?T>0L-ywSC)OZwBaz#($t8oBt1_6 z_rI+*q|;$1dlejd2<K7mkF*rAWuZM-r{H)}vvBUy#8s2%lz9#nfW^jDa$%ICs4^6~O* zxU#*-(2p?G|u-*AC4T$7OjBWQm+?X>-sVfwLAlX-1P3mzc25Nu>UBGB4Gf1r2h z)8}6JGm6il6!q_gjLp8OfK{z@=?BM_*Uul8Hj9?_;Ao@S(W3eH^AU@n&>Egfu24|BDC|abq{uU$O&TAENV9|Z{e-i#qF)`5P7;04 zg&lau&Zkw#RdN_`htI~?A~mF^Cu~xwPNI<{lkxK`%pjzaMyjZ6l^Y`8-lv4b4XDH7 z?ppNd$VjSga9DFC@2m=s-Y?FbYj7#Hg7$C_T}h7%4xf_0m4frLc)j>8YSes%qm zx96!)qTzW_XT*TQyXg-sO zh@ff3gtLb6t=wbqyB#j)2bLaY>0ZGcz^|AMl@ZB*FRQVwUbX&o@D_FZ*(cSouW*#R zhsVdBo;Y43=r|}VvkwReP&Rzs)Knj6EBDD|itk2sWN~lr&o0P+ z1f^W;sb#1C>)HrWZeidPN9T|j*~2q28aXth=;Gy-s#LfW3V^k;vJx~jyZc=84Xs<3 zY-#F~T=VirESc&?hOO%cHTLEPmI@|Tki+a4f>q-5lcMo`xn+%zv$h)OjZ5scguzD> zMu~}ui~YPNZyIff@^W%=fTj~_o`;kIuBZ1cz~rd^!Lf^mZ)_sgQxf_2Vx>Bg`hEA0 zr5S&IS(6tp+gm$&lrd#R(fq2qiJLN3gNETig6GbI^QBlCe?HoIF?^7<-Osj|qRD0Y zsuqmasQeCJ1q9!box60W1v*H;Y=7Z$2*$LorS?IA_s7s~PH6JT@9VRcT)6wx|NGApJE~avFEO-s9f@Z!+>cj~;eeV+ zW9qjITW2`L@x*S8X{1>}@eX~~ z@j)%Bi(#ugCxwI_^UN7BAM&D7l_=2$nilLr4?}k#Nkg>_!iWLyok#5t^jwcm&i*jo z(A>_qKGUjvas1yWo#YB8H6zwCL)x=v<5N@p>gr~|j-gF~#`L!y4c?X~BRM*{W6NG) zs*WNj$vwPdw{o`eYWv`9^TMKC-&{`9V^{Ax@ui#h(Kj`?g8oQB3^#)AdZKFmj&=Z3 zf)BA3h*YC_XyR!F2K66lSE3sSK@7~PvG^QB-~Io=f+Z7zd^1B**CN=jUNr;G35w?X z4h|PDUgYNGovXSRuFnl+cwMV37{J&0qLn>S+W(V$I0z&flEZ_O7>6M3jhDJ1QzO(M z2tubQ_TR4|gh_GCqhNPag)1<_l{*iKx;?vi@|8li?pKJuU~;> zpsMNtFo;GnGBl+CBNDQ*Ho(0MxE0-j2mOq^cg9k5UNfE6k#Dhgoa z3ox&O&kO91d6ktLKqT(ipFiI~uvh;iOXx$UnBI|M^$rwBM#jc)zg(sZKYj?^>%N$Y zZHs`{V`h4~JC^?JaQmLFrkXJ?@_B!=!Y?2pBhun?eP17f{D`R2`=P@yit z^56)R)J63;+P35UgzupJ2rTa!ci6xb4?Oey6a40dK=7LFKOK-Rc<}<osUz&`+{hfAcSH=ysXlGw2ryzng%#RUbZ)~)|IyU%FX<@iglQbyU^*+C5t zf+VNur_q8>Jw0a5q0?bF&*Q*Q9$5@d75ad~+=hmOI0jiDNZi1-T=`gq)pzUbEil=W zH)%l4#T@j43KJS5_4OgBy8LxR>sA^aL80sn3rl|k!^T$bDM++yQw;;NvnHoh?2~b^ zdM4go;By2L1*|a^EU+>mhZ0@9D8SFpB)}xaBQdlKr5p31kigfwo}jQ?AOTbH7?jOp zV-2OHE8zQtd+mY#<~lnHv<(1k-$!h^d}zyx*35fufwjEng@ub9p%=6?HKD$ToeVC7 z-b6OoVIKnoaxV#1gS80^($Kpbm7yzZD;oeNju?8WR_l-HaM;xNINB+Wss@L)Rewqq zzOCE{utU(^TvJ^QMaFCUE|mLX9|a?DDNJu0axo6#U_+^+Kp9q7Q2`F#1~_skJ(|2W zPoa&_%G#O~dVI+Q52D?4^Y+6sVj#u+b+xR*=z|6RALa?&>Igdxs!_2PzY4VYfH(>+byY~&`gmD&DLYPHZPKOXtUbgIC{uBA% zZ0(FP{5!8Kdr2Yo)YjI*gar)cDRAeiAqY_r`ndqNwO*)ruU0h7La!U6#To0l*|4w^ zHtoZ?tcY<5_qA-^W?|icn;{_OfT!|)N$~k|P9S5Uz~r&9y2@qxNiG>0i19uaM?x|o z8yYUrqUak6vx6RXutDe>8UkaGCB`EvDtgVZzC1Vg(Mh*Fe<>*K9?OHdz-;WVj`1=v z`p=fZHr!%!&bj$tUQ;mR$)3r)xW7Rj%`k3&K7y199kQ)PYyrb&U> z;)L1TL6dB?E&q~t6BPUqG@F}!Ay`mILX^}oHl7(Dzv(c}7iyXYunR2}VZn&pGdHh- zj>)5BUFZdfv>nX3gGNzpCe6NUmz^bW<};K2Y6`&zZRc54T|=H|IFB4QrZ^ywki+H64FfL8=R z6#8XpNx0}fvBcyl-Nff*Ni$iY!i9**j8l11wX0Zg`nUS03G+@jp?14k>KK5MO)<_ zd<%ov3@aVDnAGHCXqxX_^nm+hQ+v*S?V3rQT{67MPLZu_9q3V@J|FCjnSvJ)+7!T< z1%VUx*2MewO8}lNAxD>%BDbhxYXS(Jx=zfXu@$dx2k1Kp@6&eQ|LSkcYdAi$;~UBUs7-FQ9#Iyn@5= ziiMxo)*OJyVvr4ES9dcTjCKqc(k=s5nU{BNYDC6I3iEi2wzUmt-pwWY;h z*#mn#tQ*Q`r2RbgN&&z%z8n8 zS*dHR%9;r49Yj%i*2Laa_;(2iKm1U98d%3Ijg8N2#9`V+z)lDD(GIuTcOO1DfKvcy zLW`dd1)eTl02y0w7`wx>xVTp~19U`0M7lO0^8u#SrP+s0dgz~?xONjVQ$9dG zkR%N!Il*B(GU5RXTC*hS3nmP1bf^y*B26blcUkb8O= zJRvV(8$ja^tmtcZt1ST!Lz5FExqG4Xep_#0&;brN5{!7mzcM4y8T`U#4$O8+_7Pp91mw#B8>RG{2=18=pzF(U@OcGJ`1 z&_%$_T@JtqTn(Bf+OKH=i~#njso_7f{;mZHr$emE+H+H~v(Kew=Hwxfk7|Cd; z?XT$PC$8rRK9F*a*5(>`1PcIo02TV?@8Dy(d4MMY#FVukKRhAVrC5CG@4rWTKGd~N z0)P!Gu*YyTJmMhc3c}^8}V2 z9|=64Vhv>HANVbA-awxYJlm)3aN)vHA0)PEB5y{Em$-yiZo(sHD&thAeGa#BJHkK$ zpPxkw?4aEMegpd}&pFRO`;Lq&b8UfKgniS~&m*8tKcqJ^GP*)?@{RuN>kt2D&z?C< zHx3Qgrj=1);SyR3-v%m6sI|+6SDcB1hev1J7_c++p8kS+*+~l`rlUK}N&s)J(d!qn zu|y5wDhddXolx#e z>?)2LSs@ax-+Q1Yg#N&q!f{w=5n$zdviY{hTyGo#)wE1O&n3xC7Oo zW}l}xY++*XY*Gwhy?O8rhokUP_vonZ>$uE5k-v8^2P6z^Cs^sw6$1_fGBC1DO~I0P zfqU}%;f9tC)CCzWMlN5ze5f5pcn%M5SB{w!ji3MAftsYY`odyCSCa;30z4>!>||o@ za}RJ*KnEx=Ulwwp=9>2o;(l6pIlE=6m?9*{#tn?goNM;Bl!kx zi7H(J{%HPRnh2hS6QRJF9B?=vE*_}c5VO;GzS*JBB~m5p=~Dud6F}z<=Y&pMf3I*c z?(FWaL7O31$?0%=I&Al6nHGDJx9`hk#k-(eD|Ih(b8*2VRgU2~K`>F!2L`2io~L1! zer&~r39q{g_;!8FtgU}+Y^VihbuM;wb+v$e34mj$FD5BfR#jO!Un*QZy$cfzbC=iI zrKHZ#LU`v7jV*C$@K<@_#&)<2!9IkNXLuuNz~rGEy5A41%@A+gvL~me2GRxXPmhn* zmLJw)n4+b}<$+EjXv4?bWr8l|hyO-V@!?1D0kZ+!MR>-hRc>>L6z8yXrY+sLA^u1!1&_>eBhB_eVJRzxT- zpF^Gv!oY$;%MB1VtUQR2u0?x7Hv|QZp;Ow$k6>1UzVz{NG^oJE4Nnn5Lvz4% z@Cz1jf(lkb;Rs5bt21!h& z#&(D2NbOofT>&TqJpGU`q`06$?1u+_pnY+U>lWa*0)L@{epbZq{q<(o#XdYT*5X-C zKE5qje4uN5yf(ulh5L4T`V9nmP-da<+Wz%RtGM?U0+y#A^YafNv~RYZ%>sCV$4^1_ zc^Vzf&CdRnzv1zK$>A{Ez|Y~5FCgVYOJ83)*y2>v1%3t-Ne(VVhe6hdMiIix0zFqL z{c&mE+0Eqy!uH|h;i-ek0SauCad2&31$Zly1}9_iGXg6LV5s+(=NP!$pt?nGPWa!L z5?{s5b_m(2DJiG$aHL~;oA>CELi^c((STIVgQ4y*7*-X9}m# zPtPeTECkmHsJO=vxgqpg+poZLysB$0RbanJLcw@MaUSD)BuM&uv{*s=(8FKd01DVtBDO6AW2LF9S~il(&H94y4SFG#lvQ0&Lp^n9+JnPtXXOb`GI) r93aEEW(d;9R{{^vz@9(DxJIUTvv;VC4xW#PDBe~>RmzzCe`^8&<2?5I literal 0 HcmV?d00001 diff --git a/problem_1_data.m b/problem_1_data.m new file mode 100644 index 0000000..562afa4 --- /dev/null +++ b/problem_1_data.m @@ -0,0 +1,23 @@ +clear +clc + +% part a +xa=[1 2 3 4 5]'; +ya=[2.2 2.8 3.6 4.5 5.5]'; + +Za = [ones(size(xa)), xa, xa.^-1]; +[a1,fx1,r2_1] = least_squares(Za,ya); + +% part b +xb=[3 4 5 7 8 9 11 12]'; +yb=[1.6 3.6 4.4 3.4 2.2 2.8 3.8 4.6]'; + +Zb = [ones(size(xb)), xb, xb.^2, xb.^3]; +[a2,fx2,r2_2] = least_squares(Zb,yb); + +% part c +xc=[0.5 1 2 3 4 5 6 7 9]'; +yc=[6 4.4 3.2 2.7 2.2 1.9 1.7 1.4 1.1]'; + +Zc = [exp(-1.5*xc),exp(-0.3*xc), exp(-0.05*xc)]; +[a3,fx3,r2_3] = least_squares(Zc,yc); \ No newline at end of file diff --git a/setdefaults.m b/setdefaults.m new file mode 100644 index 0000000..8c44cd6 --- /dev/null +++ b/setdefaults.m @@ -0,0 +1,3 @@ +set(0, 'defaultAxesFontSize', 16) +set(0,'defaultTextFontSize',14) +set(0,'defaultLineLineWidth',3) \ No newline at end of file