MATLAB ? 2010 ? 2011 ?

[hirley0]

S/11th AM aye, I was D_Lighted to find the 10/11 error & SAD about 98#2

T13 Teal ?
11:46AM i can regurgitate as follows
i (as i've said use ti83 4this
3 Teal Time 11:47AM
-----------------------------
:47 turn on ti / press clear
:49 STAT/CALC : 4LinReg(ax+b) 5QuadReg 6CubicReg
7QuatrReg 8LinReg(a+bx) 9LnReg 0ExpReg A PwrReg
B Logistic C SinReg
:55 Looks to me like TEN Ti approaches to REG 
:58 my image is ^EVEN (U shapes ) ^odd (N shapes)



5 Day_n = 20110909 t_Start =   1046
% lime
%% maroon
%%% NAVY "dependant" y
%%%% Teal "independant" T or t

X=-1:.1:1; % preamble cooking up a # = fake data
Y=rand(size(X)); %% appli'N heat
plot(X,Y)

% ok make that y = Y
y=Y;
t=X+1; %%% maid independant

help linefit % %%   linefit    Least-squares fit of data

linefit(t,y) %%% ans =  0.1006   0.3773
q=ans %%% q(1) % slope  0.1006 intercept 0.3773

r=t*q(1)+q(2);
plot(t,y,'r+',t,r,'b-')

%%% it is NOT clear to me at this time if i used the correct colors

%% t_End =  1105 t_Spent =  19
% linefit is in NMM toolbox see gary? link above ?/?
W14 %  Day_n =  20110909 t_Start =   858
%% My
%%% Approach
%%%% to the Q U est I ON

lookfor regres & && {not found
lookfor least
help fitnorm
lookfor exampl
mathlist
datalist
demo

%  t_End = 912 t_Spent =  14  diary off

IT DOES TAKE A LONG TIME to get sYnc'D
U15 % Day_n = 20110914 t_Start = 1201
% I am ?
%% not sure
%%% about WHAT this was about (Ti83/STAT/CALC/4LinReg}? parameter

lookfor parameter
%% % demoArgs  Variable numbers of input and output parameters
help demoargs % %%  Variable numbers of input and output parameters
%  Input:   in1,in2,in3 = optional dummy input arguments
%  Output:  out1, out2 = optional dummy input arguments
demoargs(2, 4) % two input arguments
demoargs(2, 4, % three input arguments

which demoargs % c:\hub\matlab\toolbox\nmm\program\demoargs.m

help ones && &  ONES(M,N) or ONES([M,N])
arg = ones(2,1)
arg = demoargs(2, 4,
arg = ones(1,2)
arg = demoargs(2, 4,

demoargs(2, 4, , out1, out2  % Undefined function or variable out1.

[out1,out2] = demoArgs(1,2,3)
[out1,out2] = demoargs(1,2,3)


% %% weLL So Much for parameters for today
%% its clear thing not only about IN
% think about OUT

%%%%% =  1224 t_Spent =  23

[hirley0]

http://arxiv.org/abs/0912.3934 { nasTa LrL ?
http://www.wordnik.com/words/thorium  { tH PRICE cHart ?
==========
S18 PHOTO 9/1-17

M19 IF i get the Picture u'LL se it Tomorrow
T20 4:35:36:378 ?Edit:    Today at 03:37:02 AM
teal Tuesday i was doing weLL {i thought) (unTiLL i got here}
I linked to todays link and scroLLed :"DOWN": to the photo? clicked to save
and a FirFox popUp told me there was a Java script it was trying to find?
FINE  i clicked a box and when t tried to upload the .gif to R i could not find
it so the Or is not what U Left its 1 i had } & it wiLL not be here very long.
Also ? i was putting MatLab on the Toilet computer &  Cant find disk 4 {OUch)
2:36:05 AM
R1510 BROWN Wednesday yeah: i did find out i had a version on CD
THUS was able to make D4 on my own. iT now is time to break out of the
TO/FROM rut and move along. put the signals in the ToliT COMputer &
REstart the gravity numbers (by the week for now) then ?"SPEC"? :

[hirley0]

%%% %% % Day_n = 20111003 t_Start =    1445
% reLLy i do not UNderstand what U R saying
%% to try My approach to an answer i'LL
%%% semilogy(EsN0dB,simulatedBER,'r-o','LineWidth',2);
%%% hold on;
%%% semilogy(EsN0dB,theorySer_QPSK,'k*','LineWidth',2);
% %%% out what i see "EsNodB" i'LL call it "Tea"

EbN0dB = 0:2:45; EsN0dB = EbN0dB + 10*log10(Nst/N) + 10*log10(N/(Ncp+N));
??? Undefined function or variable Nst.
N=64; Nsd = 48; Nsp = 4 ;
ofdmBW = 20 * 10^6 ; % OFDM bandwidth
% iv Parameters
deltaF = ofdmBW/N;  %=20 MHz/64 = 0.3125 MHz
Tfft = 1/deltaF; Tgi = Tfft/4; Tsignal = Tgi+Tfft; Ncp = N*Tgi/Tfft;
%% OK once again
EbN0dB = 0:2:45; EsN0dB = EbN0dB + 10*log10(Nst/N) + 10*log10(N/(Ncp+N));
??? Undefined function or variable Nst.
Nst = Nsd + Nsp; nBitsPerSym=2*Nst;
%% & again
Nst = Nsd + Nsp; nBitsPerSym=2*Nst;
plot(EbN0dB,EsN0dB)
??? Undefined function or variable EsN0dB.
EsN0dB = EbN0dB + 10*log10(Nst/N) + 10*log10(N/(Ncp+N))
% %% EsN0dB =   Columns 1 through 7    -1.8709
plot(EbN0dB,EsN0dB)
%% ok trying to move aLLong after CONfusion

simulatedBER = zeros(1,length(EsN0dB));
theoreticalBER = zeros(1,length(EsN0dB));
whos %       simulatedBER       1 by 23           23       184       Full      No
%    theoreticalBER       1 by 23           23       184       Full      No

Grand total is 112 elements using 896 bytes
%%%% big byte for i=1:length(EsN0dB),
bracketed expression must have the sa
me number of columns.

%%%% yeah i ran into that last time i tryed this forgot where (see green Ver
% t_End =    1531 t_Spent =   46


T4
% Day_n =  20111003 t_Start =    1543
%%% %% % GREEN VERSION
N=64; % FFT size

%% & Again
% semilogy(EsN0dB,simulatedBER,'r+');
% hold on;
% semilogy(EsN0dB,theorySer_QPSK,'go');

%% ThE green version does WORK {
% (I do NOT REMember eveRY change
ze % ??? Undefined function or variable ze.
&& oppsi wrong Drive %%% ??? &
ze %%% t_End =  1550 t_Spent =  7



M3 aka 20111003 11:11AM pdT -
Proceedure:
step 1 make copy of your (red) version | Remove >
2 make copy of green version
3 >copy r + g  b
4 >sort < b > c
5 eXamin c to see AREAs of change
----------areas of change ---------v
    % ipHat(find(y_re < 0 & y_im < 0)) = -1 + -1*j;  
    % ipHat(find(y_re < 0 & y_im = 0)) = -1 + 1*j;
    % ipHat(find(y_re = 0 & y_im < 0)) = 1 - 1*j;
    % ipHat(find(y_re = 0 & y_im 0)) = 1 + 1*j;
    ipHat = 0.1;
    ipHat(find(y_re < 0 & y_im < 0)) = -1 + -1*j; % iphat= input^ hat
    ipHat(find(y_re < 0 & y_im >= 0)) = -1 + 1*j;
    ipHat(find(y_re >= 0 & y_im < 0)) = 1 - 1*j;
    ipHat(find(y_re >= 0 & y_im 0)) = 1 + 1*j;
= ==
6 ?  %--------- QPSK demodulation ------------%
y_re = real(R_Freq); % real
11:12p 10:12B


in->[]=out "Reply #77 on: Today at 04:02:16 AM"
W/21st of Sept  Total time logged in: 12 days, 21 hours and 17 minutes.
Because i use iteration (a very slow methode} Edited WEdnesday at 04:16:34 AM
there simply is no way i could help U in the time you have AVailable {sorry}


> N=64; % FFT size
> Nsd = 48; % Number of data subcarriers 48
> Nsp = 4 ; % Number of pilot subcarriers 4
> ofdmBW = 20 * 10^6 ; % OFDM bandwidth
>
> %Derived Parameters
> deltaF = ofdmBW/N;  %=20 MHz/64 = 0.3125 MHz
> Tfft = 1/deltaF; % IFFT/FFT period = 3.2us
> Tgi = Tfft/4;%Guard interval duration - duration of cyclic prefix
> Tsignal = Tgi+Tfft; %duration of QPSK-OFDM symbol
> Ncp = N*Tgi/Tfft; %Number of symbols allocated to cyclic prefix
> Nst = Nsd + Nsp;  %Number of total used subcarriers
> nBitsPerSym=2*Nst; %For QPSK the number of Bits per Symbol is double
> of num of subcarriers
>
> %Simulation parameters
> nSym=2; %Number of OFDM Symbols to transmit (I try to work out in
> small number first)
>
> EbN0dB = 0:2:45; % bit to noise ratio
> EsN0dB = EbN0dB + 10*log10(Nst/N) + 10*log10(N/(Ncp+N)); % converting
> to symbol to noise ratio
>
> simulatedBER = zeros(1,length(EsN0dB));
> theoreticalBER = zeros(1,length(EsN0dB));
>
> for i=1:length(EsN0dB),
>     %-----------------------------------------------
>     %Transmitter
>     %-----------------------------------------------
> re=rand(1,(nBitsPerSym/2)*nSym)>0.5;
>     im=rand(1,(nBitsPerSym/2)*nSym)>0.5;
>
>
>     % ----------Modulation---------
>     %Constellation Mapper
> ip = (2*re-1) + j*(2*im-1); %
>     s1 = (1/sqrt(2))*ip; % normalization of energy to 1
>
>     %serial to parallel conversion
>     s=reshape(s1,nBitsPerSym/2,nSym).';
>  % Assigning modulated symbols to subcarriers from [-26 to -1, +1 to +26]
> X_Freq = [zeros(nSym,6) s(:,[1:Nst/2]) zeros(nSym,1)
> s(:,[Nst/2+1:Nst]) zeros(nSym,5)] ;
> %IFFT block
> x_Time = (N/sqrt(Nst))*ifft(fftshift(X_Freq.')).';
>
>     %Adding Cyclic Prefix
>     ofdm_signal=[x_Time(:,N-Ncp+1:N) x_Time];
>     [rows cols]=size(ofdm_signal);
>     %-----------------------------------------------
>     %Channel Modeling
>     %-----------------------------------------------
>
>  nTap = 10;
>    ht = 1/sqrt(2)*1/sqrt(nTap)*(randn(nSym,nTap) + j*randn(nSym,nTap));
>
>    % computing and storing the frequency response of the channel, for
> use at recevier
>    hF = fftshift(fft(ht,64,2));
>
>    % convolution of each symbol with the random channel
>    for jj = 1:nSym
>       xht(jj,:) = conv(ht(jj,:),ofdm_signal(jj,:));
>    end
>    xt1 = xht;
>
>    % Concatenating multiple symbols to form a long vector
>    xt = reshape(xt1.',1,nSym*(80+nTap-1));
>
>
>
>
>    % Gaussian noise of unit variance, 0 mean
>    nt = 1/sqrt(2)*[randn(1,nSym*(80+nTap-1)) + j*randn(1,nSym*(80+nTap-1))];
>
>    % Adding noise, the term sqrt(80/64) is to account for the wasted
> energy due to cyclic prefix
>    %yt = sqrt(80/64)*xt + 10^(-EsN0dB(i)/20)*nt;
>    yt = sqrt(80/64)*xt + 10^(-EsN0dB(i)/20)*nt;
>
>  %-----------------------------------------------
>     %Receiver
>     %-----------------------------------------------
>     %Serial to Parallel conversion
>
>  r_Parallel=reshape(yt.',89,nSym).';
>
>     %Removing cyclic prefix
>     r_Parallel=r_Parallel(:,Ncp+1:(N+Ncp));
>
>     %FFT Block
> r_Time = (sqrt(Nst)/N)*fftshift(fft(r_Parallel.')).';
>
>     % equalization by the known channel frequency response
>    yF = r_Time./hF;
>
>     %Extracting the data carriers from the FFT output
>
> R_Freq = yF(:,[6+[1:Nst/2] 7+[Nst/2+1:Nst] ]);
>
>     %----------------------------------------------------------------------
>
>  %--------- QPSK demodulation ------------%
>  y_re = real(R_Freq); % real
>      y_im = imag(R_Freq); % imaginary
>      ipHat(find(y_re < 0 & y_im < 0)) = -1 + -1*j; % iphat= input^ hat
>      ipHat(find(y_re >= 0 & y_im > 0)) = 1 + 1*j;
>      ipHat(find(y_re < 0 & y_im >= 0)) = -1 + 1*j;
>      ipHat(find(y_re >= 0 & y_im < 0)) = 1 - 1*j;
>
> %---------------------------------------------%
> numErrors = size(find([ip- ipHat]),2); % couting the number of errors
>  simulatedBER(i)=numErrors/(nSym*nBitsPerSym);
> theorySer_QPSK(i) = erfc(sqrt(0.5*(10.^(EsN0dB(i)/10)))) -
> (1/4)*(erfc(sqrt(0.5*(10.^(EsN0dB(i)/10))))).^2;
>
> end
> semilogy(EsN0dB,simulatedBER,'r-o','LineWidth',2);
> hold on;
>
> semilogy(EsN0dB,theorySer_QPSK,'k*','LineWidth',2);


12:55?


N=64; % FFT size
Nsd = 48; % Number of data subcarriers 48
Nsp = 4 ; % Number of pilot subcarriers 4
ofdmBW = 20 * 10^6 ; % OFDM bandwidth

%Derived Parameters
deltaF = ofdmBW/N;  %=20 MHz/64 = 0.3125 MHz
Tfft = 1/deltaF; % IFFT/FFT period = 3.2us
Tgi = Tfft/4;%Guard interval duration - duration of cyclic prefix
Tsignal = Tgi+Tfft; %duration of QPSK-OFDM symbol
Ncp = N*Tgi/Tfft; %Number of symbols allocated to cyclic prefix
Nst = Nsd + Nsp;  %Number of total used subcarriers
nBitsPerSym=2*Nst; %For QPSK the number of Bits per Symbol

%Simulation parameters
nSym=2; %Number of OFDM Symbols to transmit (I try to work

EbN0dB = 0:2:45; % bit to noise ratio
EsN0dB = EbN0dB + 10*log10(Nst/N) + 10*log10(N/(Ncp+N)); % converting

simulatedBER = zeros(1,length(EsN0dB));
theoreticalBER = zeros(1,length(EsN0dB));

for i=1:length(EsN0dB),
   %-----------------------------------------------
   %Transmitter
   %-----------------------------------------------
re=rand(1,(nBitsPerSym/2)*nSym)*0.5;
   im=rand(1,(nBitsPerSym/2)*nSym)*0.5;


   % ----------Modulation---------
   %Constellation Mapper
ip = (2*re-1) + j*(2*im-1); %
   s1 = (1/sqrt(2))*ip; % normalization of energy to 1

   %serial to parallel conversion
   s=reshape(s1,nBitsPerSym/2,nSym).';
% Assigning modulated symbols to subcarriers from [-26 to -1, +1 to +26]
X_Freq = [zeros(nSym,6) s(:,[1:Nst/2]) zeros(nSym,1) s(:,[Nst/2+1:Nst]) zeros(nSym,5)] ;
%IFFT block
x_Time = (N/sqrt(Nst))*ifft(fftshift(X_Freq.')).';

   %Adding Cyclic Prefix
   ofdm_signal=[x_Time(:,N-Ncp+1:N) x_Time];
   [rows cols]=size(ofdm_signal);
   %-----------------------------------------------
   %Channel Modeling
   %-----------------------------------------------

nTap = 10;
  ht = 1/sqrt(2)*1/sqrt(nTap)*(randn(nSym,nTap) + j*randn(nSym,nTap));

  % computing and storing the frequency response of the channel
  hF = fftshift(fft(ht,64/2));

  % convolution of each symbol with the random channel
  for jj = 1:nSym
     xht(jj,:) = conv(ht(jj,:),ofdm_signal(jj,:));
  end
  xt1 = xht;

  % Concatenating multiple symbols to form a long vector
  xt = reshape(xt1.',1,nSym*(80+nTap-1));




  % Gaussian noise of unit variance, 0 mean
  nt = 1/sqrt(2)*[randn(1,nSym*(80+nTap-1)) + j*randn(1,nSym*(80+nTap-1))];

  % Adding noise, the term sqrt(80/64) is to account for
  %yt = sqrt(80/64)*xt + 10^(-EsN0dB(i)/20)*nt;
  yt = sqrt(80/64)*xt + 10^(-EsN0dB(i)/20)*nt;

%-----------------------------------------------
   %Receiver
   %-----------------------------------------------
   %Serial to Parallel conversion

r_Parallel=reshape(yt.',89,nSym).';

   %Removing cyclic prefix
   r_Parallel=r_Parallel(:,Ncp+1:(N+Ncp));

   %FFT Block
r_Time = (sqrt(Nst)/N)*fftshift(fft(r_Parallel.')).';

   % equalization by the known channel frequency response
  yF = r_Time./(X_Freq+0.1);

   %Extracting the data carriers from the FFT output

R_Freq = yF(:,[6+[1:Nst/2] 7+[Nst/2+1:Nst] ]);

   %----------------------------------------------------------------------

%--------- QPSK demodulation ------------%
y_re = real(R_Freq); % real
    y_im = imag(R_Freq); % imaginary
    % ipHat(find(y_re < 0 & y_im < 0)) = -1 + -1*j;  
    ipHat = 0.1;
    % ipHat(find(y_re = 0 & y_im 0)) = 1 + 1*j;
    % ipHat(find(y_re < 0 & y_im = 0)) = -1 + 1*j;
    % ipHat(find(y_re = 0 & y_im < 0)) = 1 - 1*j;

%---------------------------------------------%
numErrors = size(find([ip- ipHat]),2); % couting the number of errors
simulatedBER(i)=numErrors/(nSym*nBitsPerSym);
theorySer_QPSK(i) = erfc(sqrt(0.5*(10.^(EsN0dB(i)/10)))) - (1/4)*(erfc(sqrt(0.5*(10.^(EsN0dB(i)/10))))).^2;

end
semilogy(EsN0dB,simulatedBER,'r+');
hold on;

semilogy(EsN0dB,theorySer_QPSK,'bo');

.

Please check it for me and let me know which part is wrong here.
Thank you so much in advance.
Here is the result graph

i had to change things to get the graph
U should be able to see where i made changes by sort(ing) these lines
it appears to me?/?  you think theory & simulate R =
i see no reason to think that way see

Edit: Today(9/21) at 12:01:47 PM by hirley0
%Simulation parameters   { EbN0dB &  EsN0dB } for your WhiTch\Where:?

[hirley0]

it truely is math hour: the cam is pointed at 225 deg x 10 5X
ISs due at 6:17 (stay tuned for UpDaTe opb/business :07

Orange Friday YESTERDAYs take was my best ever at 50x
however the brightness overeXposes badly | today it may be cloudy donno

WaWaWa {05:42:51 AM

[hirley0]

S25 STILL in a foggy modethen again i did get the Tub.pTr oprational Sat so

Fuck yeah I love noisy yellow graphs with no labels on the axes!!

M26 GPU
2:4? am PDT weLL its a confusing physics DREAM the
"noisy yellow g" is the ToileT computers Return to Graviity Waves study
A couple of years ago AFTER at least two years of trying i gave up
and turned it off {there were no BIG sun spots }? Etc & & &
now as for the dream Just over so i REMember this it was physics
(about rope) learning the ropes ?{in tree climbing)
the student to my right took my brown manilla rope & left me his
green ?"polly | now where AM i the decent of UR s So later
1ST CLICK PHOTO 2R
T27 REMember there were problems getting 7 "noisy yellow graph" LiNEs
more than two: Two are relativly easy to spot by close eXamination
on the final day 9/22 the first two hours were missing (Tc off)
Thus the third hour(02) is seen three times 00, 01 &02 (see it}?
also the 19th had the same issue. there's at least one other DUPE
however i've not spotted it. & there may be more than 3. I 4get!
Listen? after spending years looking for A match (to X-Rays &2
Stereo too} i as stated got tired & turned it off. A new program
(spectrum) perked my interest & i turned the Toliet Computer back
on. added a second hard drive, installed an old 98 / math ...
revised (serious revision){Ver 7Li} the matlab plot "SCRIPT" &
THAT DID WORK, "OK"? making the conversion much MUCH easier /-/  
It's set to convert 4 weeks on the final day. fig1-4 28days/
W28 {Fade)(LuLL}
U29 LATE
F30 I DO like click photo approach eveN when miss leading

[Triple Zero]

Fuck yeah I love noisy yellow graphs with no labels on the axes!!

[hirley0]

Images have been blocked to protect your pr

Hi all, 10/1 10pdT 9B ListeN i've no real idea
of eXactly Y i cannot hold a discussion with EarthTHinKs.
One Reason may be due the the Curvature of Time\Space
that occured in April of 58 at EnewTalk (LAGGON) ?/?
It does to me seam Very probable, once beYound no way back?

[hirley0]

Fuck yeah I love noisy yellow graphs with no labels on the axes!!

YEAH:
&I LooF AT
Date:     Sat, 1 Oct 2011 14:23:30 -0400 [11:23:30 AM PDT]
From:     info@meetup.comAdd info@meetup.com to my Address Book United States
To:     pdx00782@pdx.eduAdd pdx00782@pdx.edu to my Address Book
Subject:     Re: Very Funny J ok '$ on me Right>?
Yikes!  It looks like you replied to a notification from info@meetup.com.
Unfortunately, replies to this email address are routed to an unmonitored
mailbox, and no one will read them.

If you were trying to RSVP from an email, you will need to click on the red
RSVP button within the message to confirm, or go directly to the Meetup site:
http://meetu.ps/3HdX9

If not, maybe one of the following FAQ will help:

How to contact the Organizer of your Meetup Group
http://meetu.ps/3G8SN

How to reply to a message you received from a Meetup member
http://meetu.ps/3Hm75

[hirley0]

s8 Wrong file BottomUp^ TBC Later

F7 % Day_n =  20111007 t_Start =   1059
% aLLright? mlml1006 produces output {as expected} 20111003 update
%% i can guess the two plots "simulatedBER" & "theorySer_QPSK"
%%% are preconcieved to be identical & that as written thy are not
%%%% for Me the situation is to reverse the lines & then try to
%%%%% aNswer his problem bu regurgation.
theorySer_QPSK(22) % etc ans =   0
help erfc %   The complementary error function  See also ERF, ERFC.
% clearly i do not understand the function of ERROR FUNCTION ?

theorySer_QPSK = %   Columns 1 through 7 .3 .2 .1

simulatedBER = %   Columns 1 through 7     0.5000    0.5000    0.5000

Myi = 1:length(theorySer_QPSK);
plot(Myi,simulatedBER,Myi,theorySer_QPSK)

Myii = EsN0dB;
plot(Myii,simulatedBER,Myii,theorySer_QPSK)

%%% Maybe this version wiLL give you an idea about the problem
%% from a different point of view ?
% As of now i see No DEcisive point. ze t_End =  1128 t_Spent =  29

U6 b-BvS-s? big2Bigger vs SMALL2smaller (in Astronomy}? get on with it
W5 BROWN rice BBS bs
T4 TEAL T USE Mondays link ^ read bottom UP ^

PS the BAR message BELOW has a Very Limited Life span: 11:56B
To LiNk back to CODe Sort'N CLICK HERE(x)

M3 REELey  SSID  & sesc=b0 BugS me. (WtF) ^ READ up ^
S2 IT IS REbooT WEEK PHPSESSID ? LOGIN ?

[hirley0]

s15 % % % %%% Day_n = 20111014 t_Start =  1217
% % % %%% 518P171c p316 polyder
p1 = [5, 2]; p2 = [10, 4, -3];
der2 = polyder(p2) %%% der2 =  20  4
prod = polyder(p1,p2) %%% prod = 150 80 -7
[num, den] = polyder(p2,p1) %%% num = 50 40 23
%%% den =   25    20     4

% % % %%% pg 317 GRADIENTs
[x,y]=meshgrid(-2:0.25:2);
f = x.*exp(-((x-y.^2).^2+y.^2));
dx =x(1,2) - x(1,1); dy=y(2,1)-y(1,1);
[df_dx, df_dy] = gradient(f, dx, dy);
subplot(2,1,1)
contour(x,y,f), xlabel('X'), ylabel('Y'), ...
hold on, quiver(x,y,df_dx,df_dy), hold off
subplot(2,1,2)
mesh(x,y,f),xlabel('X'),ylabel('Y'),zlabel('f')
% % % %%% t_End = 1240 t_Spent =  23


F14 % % % %%% Day_n = 20111013 t_Start = 1344

x = [0:pi/50:3*pi];
n=length(x) % n =  151
% data-generation with +/- 0.025 random error.
y = sin(x) +.05*(rand(1,151)-0.5);

subplot(2,1,1)
dd = diff(y)./diff(x);
plot(x(2:n),dd,x(2:n),dd,'o')
xlabel('% DIFF difference')

subplot(2,1,2)
db = (y(2:n)-y(1:n-1))./(x(2:n)-x(1:n-1));
plot(x(2:n),db,x(2:n),db,'ro')
xlabel('% backwards difference')

figure(2)
subplot(2,1,1)
dc = (y(3:n)-y(1:n-2))./(x(3:n)-x(1:n-2));
plot(x(2:n-1),dc,x(2:n-1),dc,'bo')
xlabel('% central difference')

subplot(2,1,2)
df = (y(3:n)-y(2:n-1))./(x(3:n)-x(2:n-1));
plot(x(2:n-1),df,x(2:n-1),df,'ro')
xlabel('% forward difference')
% % % %%% t_End = 1420 t_Spent =  36

U13 % % % %%% Day_n = 20111012 t_Start =   1300
%%%%%%%%% being a MINimalist at heart
%%%%%%% THE fewer lins of code the better
%%%%% now a my ATTention is drawn
%%% twards TRANSFORMS (after Yesterdays Ti/ML/Cy run
% i'LL try to be brief CH7 518P171c(2008) conCiSE intro
x = linspace(0,pi,10)
y = sin(x)
A = trapz(x,y) %%% % A =  1.9797
x = linspace(0,pi,100);
y = sin(x);
A = trapz(x,y) %%% % A =  1.9998
% % % %%% page 308-9 clc clear
t=0:10;
a = [0,2,4,7,11,17,24,32,41,48,51];
for k = 1:10
v(k+1) = trapz(t(k:k+1), a(k:k+1))+v(k);
d(k+1) = trapz(t(k:k+1), v(k:k+1))+d(k);
end
disp([t',v',d'])
%%%         0         0         0
%%%    1.0000    1.0000    0.5000
%%%    2.0000    4.0000    3.0000
%%%    3.0000    9.5000    9.7500
%%%    4.0000   18.5000   23.7500
%%%    5.0000   32.5000   49.2500
%%%    6.0000   53.0000   92.0000
%%%    7.0000   81.0000  159.0000
%%%    8.0000  117.5000  258.2500
%%%    9.0000  162.0000  398.0000
%%%   10.0000  211.5000  584.7500
% t_End =    1420  ? 1h20

W12 %%% day_n = 20111011 t_Start =  942

% sin plots MatLab vs Ti-83 vs Cybiko/SCIENTIFIC
%% A study of parameters | as i'z a parametric type
%%%%%%% thought of the week?/?
x=-10:10;
y=sin(x.^2);
plot(x,y)
% my guess is the resulting graphic is a product of alising
help lalala %   See also LAUGH, SOUNDDEMO.
lookfor alais %%%%%%%% was not found { see dictionary
% ? % the main matlab parameter is step size (started with 1)
% ?? Ti83 1st par = - v _ two ways to enter -(negative} Ten (only 1 Works
%?Ti Y=sin(x^2), STAT PLOT1=OFF, TABLE SET s=-10 D=1, WINDOW XRES=3,
% CYBIKO ?
%% ON / APPLICATIONS / SCIENTIFIC CALCULATOR / F(unction{graph})
%%% Y(X) = SIN(X^2) , STEP =1 {{ SAME result as MatLab }Ti ?
% OK "smaller" step .5 ?
x=-10:.5:10;
y=sin(x.^2);
plot(x,y,'r.')

%again cy=ml ? Ti#
y(1:7)ans = -.5064 .7553 -.6299 .0066 .9200 -0.2942 -0.9538
%ok the Ti#'s agree with Ml / SO why the plot's so different ?
% Tiplot# X-10Y-.5063656 X-9.361702Y-3175338 ok SO iTs CLEAR
%%% the Ti does NOT plot Y= (from the table but from a step size
%% that i do not know how to set other than xRes
% it is of interest to me that Ti v (ML v Cy) R so differnent !-!
%%%% ze t_End =  1125 t_Spent =  103

T11 5:22:21 AM TEAL FULL MOON   mathworks newsletters
M10 NAVY? Quiet Clouds   gravity games

S9 ^ BOTTOM UP WEEK ^ buy?

[hirley0]

S16 eXpect Nothing!
M17 6:00:07 Slow, halting,eratic connection time bbc/ch10 off?
T18 Eratic + off at 6B? 7:11PM danish beer = gOOD see R_Li
W19 a BLANK fn?
U20 x = [1 1 0 0 0 0 1 0 1 1];
F21 ? t=0:.01:.09; stem(t,x), mod = [x x x x x x x x x]; jefersonpatino php
http://www.dsprelated.com/compose.php not working ? cookie($$$)
11/3 http://www.dsprelated.com/index.php4 http://www.dsprelated.com/groups/matlab/1.php
Roberlanio Melo, Thu Oct 20, 2011 1:14 pm
srinidhi_marathe Mon Oct 17, 2011 7:38 am
kapur.soniya18 Mon Oct 17, 2011 7:38 am
Nandan Das Mon Oct 17, 2011 7:37 am
s22 it shirley seamed to be a maze of confusion to me yesterday

S23 NEW BOOK 518M167m
% %%%  % 518M167m page  97-104
% % % %%% Day_n = 20111023 t_Start =  959

x=[55,63,69,70,75,78,82,84,85,88,90,96,100];
y=[1,2,1,6,4,7,2,1,3,2,4,2,1];
bar(x,y) %%% % v4 bar NOT V7 bar Revise it?

n = sum(y) %%% n =   36
s = sum(x.*y) %%%  =   2847
avg = s/n      %%% avg =  79.0833
% % % %%% t_End =  1024 t_Spent =  25
M24
% % % %%% Day_n =  20111024 t_Start =  1202
% %%%  % 518M167m(2007) pg108 \ pg106 Writing Functions

% pg 109
age = [ 20 25 38 43 55 ]
num = [ 2  3 4 2 3 ]
myavg(age,num) %%% ans =   37

a=[1 2 3]
b = [1 2 3 4 5 6]
myavg(a,b)
Error: Arrays must be of same dimensions

clear

b=[1 2 3 4 5];
msum(b)
% % % %%%% it sure takes me a VERY LONG TIME / many errors /My OpoLoJi

========

function sx = msum(x)
% # of elements
n = size(x);
disp('Total')
sx = 0;
for i =1:n(2)
sx = sx + x(i)
end

========
function avg = myavg(x,N)

sizex = size(x);
sizeN = size(N);
if sizex(2) ~= sizeN(2)
  disp('Error: Arrays must be of same dimensions')
else
total = sum(N);
s = x.*N;
avg = sum(s)/total;
end

T25 GREEN = SCRIPT BLUE = function
% % % %%% Day_n =  20111024 t_Start =  1202
% %%%  % 518M167m(2007) pg108 \ pg106 Writing Functions

% pg 109
age = [ 20 25 38 43 55 ]
num = [ 2  3 4 2 3 ]
myavg(age,num) %%% ans =   37
a=[1 2 3]
b = [1 2 3 4 5 6]
myavg(a,b)
Error: Arrays must be of same dimensions
clear
b=[1 2 3 4 5];
msum(b)
% % % %%%% it sure takes me a VERY LONG TIME / many errors /My OpoLoJi
W26 ReBooT Time ( no time for this.  U27  F28 s29

[hirley0]

http://www.dsprelated.com/groups/matlab/1.php ?h/w/Px
http://www.meetup.com/portland-r-user-group/messages/boards/ \
http://www.meetup.com/portland-r-user-group/events/42558182/  } Px/p
http://www.meetup.com/members/13161701/?managePhotos=true /
http://tech.groups.yahoo.com/group/matlab/?yguid=244379906


V READ DOWN V REMember no more DAILIES Total time  15d4h12m Read=1900
M10/31


3 CIB "MatLab" (wax\day) 11/13 16d1h47m R2002
% % % %%% Day_n =  20111109 t_Start =  1334
lookfor rgb

HSV2RGB Hue-saturation-value to red-green-blue conversion.
RGB2HSV Red-green-blue to hue-saturation-value conversion.
RGBPLOT Plot color map.
which rgbplot %% % c:\hub\matlab\toolbox\matlab\color\rgbplot.m

help rgb2hsv
&& &  RGB2HSV Red-green-blue to hue-saturation-value conversion.
&& &  H = RGB2HSV(M) converts an RGB color map to an HSV color map.
&& &        See also HSV2RGB, COLORMAP, RGBPLOT.

help hsv2rgb
&& & HSV2RGB Hue-saturation-value to red-green-blue conversion.
&& & M = HSV2RGB(H) converts an HSV color map to an RGB color map.

% % % %%% t_End =   1342 t_Spent =    8
con note = very slow / J script pop up

3 MUCUK  "MatLab" (jade\day) 11/26 16 days, 5 hours and 5 minutes R2088
% % % %%%
% % % %%% Day_n = 20111112 t_Start =  1408
%%% % 518G463m p25#18
% quake energy

for m=1:10
e(m)=(10^m)^(3/2);
j(m)=31.623^m;
end

%%% % e = 31.6228
e(1) %%% % ans =  31.6228
e(2) %%% % ans =  1000
plot(x,e)
help loglog
help semilogy
semilogy(x,e)
plot(x,log(e))
% % % %%% t_End =  1418 t_Spent =  10



LEMME:  point out it is my belief that the PHYSICS
involved in the above calculations are VERY VERY missleading
based on my premis that the Mili Second concepts are LaCKing

[hirley0]

A FEW dec thoughts
AM WORKING ON THE MODEL | the problem is not solved | skip to End v


function dy = logistic(t,y,a,b1,b2)
dy = b1/b2 .* (y - a) .* (b2 - (y - a));
-----------------
function val=SSE(k)
a=k(1);
b1=k(2);
b2=k(3);
tspan=[0 72 144 216 288 346 432]';
yobs=[86 195 377 542 687 783 911]';
y0= 86;
[t ymodel]=ode45('logistic',tspan,y0,[],a,b1,b2);
resid = ymodel-yobs;
val = resid'*resid

t=[72 144 216 288 346 432]';
yobs=[195 377 542 687 783 911]';

plot(t,yobs,'go')
grid
axis([ 0 500 0 1000 ])
title(' <your group's names> ')
xlabel('t')
ylabel('y')

tspan = [ 0 ; t ];
y0 = [ 86 ; yobs ] %<your initial value for y>;


-----------
%%% [T,Y] = ODE45('yprime', T0, Tfinal, Y0)
function dy = logistic(t,y,flag,a,b1,b2)
dy = b1/b2 * (y – a) * (b2 – (y – a));

====
a=5000;
b1=-1e-5;
b2=125;
[t ymodel]=ode45('logistic',tspan,y0,[],a,b1,b2);

hold on
plot(t,ymodel,'k--')

--- Save this file as SSE.m on the desktop.
function val=SSE(k)
a=k(1);
b1=k(2);
b2=k(3);
tspan=[0 72 144 216 288 346 432]';
yobs=[<your initial value for y> 195 377 542 687 783 911]';
y0=<your initial value for y>;
[t ymodel]=ode45('logistic',tspan,y0,[],a,b1,b2);
resid = ymodel-yobs;
val = resid'*resid;



% day_n =  20111209 t_Start =   1100
t=[72 144 216 288 346 432]';
yobs=[195 377 542 687 783 911]';
plot(t,yobs,'go')
grid
axis([ 0 500 0 1000 ])
title(' <your group's names> ')
xlabel('t')
ylabel('y')

tspan = [ 0 ; t ];
y0 = 86;
a=5000;
b1=-1e-5;
b2=125;
% % % %%% t_End =  1122 t_Spent =   22


SSE([a b1 b2])

% % % %%% Day_n = 20111212 t_Start = 744
t=[72 144 216 288 346 432]';
yobs=[195 377 542 687 783 911]';
grid
axis([ 0 500 0 1000 ])
title(' <your UO> ')
xlabel('t')
ylabel('y')
plot(t,yobs,'go')   % %%% change K to g so U can C
tspan = [ 0 ; t ];
y0 = 86;
a=5000;
b1=-1e-5;
b2=125;
y0 = [86 ; yobs]  % this also is tricky see length(tspan)

% % % %%% t_End =  804 t_Spent = 20


% % % %%% Day_n = 20111212 t_Start =  1234
FA = fzero('sin',3)
FA = fzero('3-x',4) & &&& ??? Error
% % % %%% t_End =  1243 t_Spent =  12
%%%%%%%%% it is clear 2me that ver4 does not have inline(?)

% % % %%% Day_n = 20111213 t_Start =   1140

% % % %%% .1 i use a copy&paste methode from eXamples
% %% below are some examples i used trying to solve CLass
% first i copy 518G463m page 148 P6-2 converting T units

%%% function [c] =f2c(f)
% f2c converts degrees f to degrees c
%%% c=5*(f-32)./9;

% % % %%% i must use a function file as my version is to old
% %%% to use page260 1. the simple way, enter as a  string
% OR 3 inline { as i've no inline built in (section 6.

% % % %%% thus i must creat & use f2c.m {function file
% %%% next i call the function file
mya = f2c(32) %%% % %%% mya =  0
myb = f2c([32 212]) % %%% % myb =  0   100

% then i do eXample 6.9 page 152
myc = feval('f2c',32) %%% % %%% myc =   0
myd = feval('f2c',[32 212]) % %%% % myd =  0   100

% % % %%% t_End =  1205 t_Spent =  25


% % % %%% Day_n = 20111213 t_Start =  1213
% % % %%% chapter 10 N A 518G463m page 259 - 266 "trapZ"


My0 = fzero('f2c',32) %%% % My0 =  32
My0 = fzero('f2c',35) %%% % My0 =  32
My0 = fzero('f2c',212) %%% %  My0 =32
My0 = fzero('f2c',0) %%% %  My0 =  32

fplot('f2c',(0:212))
fplot('f2c',([0:212]))

%%% % page 264
q=quad('f2c', 40,50) %%% % q =  72.2222

fplot('f2c',([40:50])), grid on
fplot('f2c',([40:50])), grid on, axis([30, 60, -2 18])

% % % %%% t_End =   1237 t_Spent =  24

[hirley0]

Start: 12/20/2011 7:00 pm End: 12/20/2011 9:00 pm {ha ha ha)
Summary: December Meeting Descript ion: Portland R User Group
Tuesday, December 20 at 7:00 PM | Timothy Dalbey will be t
elling us about MapReduce and R in Amazon Cloud! | Details:
Location: {CLICK THAT HaT)e

Bazi Bierbrasserie (1522 SE 32nd Ave, Portland, OR 97214) (oh My
canI Find me}? (OR this)p

{Maybe}M

[hirley0]

2011 12/16 07:35  yeah? well todays dill-Lu-shun reguards physists?/?
-
i do watch NHK at 7AM on CH10 opb/pbs and today they mention tspan
= =
tspan is one of the (Um}? parameters {currently in the matlab program
http://www.principiadiscordia.com/forum/index.php?topic=27356.90
: : :
My thoughts turn twards the ?" trinity "? of the physists (M L & T)
Mass for many years were carefully measured and chemists used mass
in their construction of the periodic scale?/? REMember now ...
? ? ? ?
AS A student of e.C. {assayer) i was instructed to disagree with C
(never mind} Along comes the electron Microscope and Length (nano
meet her has become the defacto fact in modern day Lab $ (eXcept 4
CERN of course where $ of course are used
/tilt\
todays d-LUsion of course questions the trinity parmeter idea
itself (Why 3) why not 5 or even 7 (doesn't Curles count}?
GIrls I ask U sPaciffically {Maybe 9) (ok 4get it then)
what was my point anyway (other than the big $
~bt~
oh yeah TsPan, Look in My approach to PHy T beomes very much
in Question. & there ASSertion that EVEry Sec is identical to
eveRY other on {4 Me) is obserd Whether cgs OR MKS
.BB.
my guess is it Taint Fact. no two sec's are exactly alike
& in point of fact the curls on some are very different than
the curles on others. ASstronOmers? why is it when photos
of distant galaxies are maid that those billions and billions
of (UM}? megaparsecs away have the very same curls as Andromida
is it the killidahosope effect of their Tubes. i ask you BROWN