%-- 4/14/07 11:20 AM --%
cp1 = csvread(‘JJRAL_10_3104to31204.csv’, 1, 4, [1 4 10 4]);
cp1 = csvread(JJRAL_10_3104to31204.csv, 1, 4, [1 4 10 4]);
cp1 = csvread('JJRAL_10_3104to31204.csv', 1, 4, [1 4 10 4]);
cp1 = csvread('C:\Documents and Settings\asesha01\Desktop\JJRAL_10_3104to31204.csv', 1, 4, [1 4 10 4]);
plot cp1
plot cp1;
plot;
plot
x clear
clear
clear;
cp2 = csvread('C:\Documents and Settings\asesha01\Desktop\JJRAL_10_3104to31204.csv', 1, 4, [1 4 10 4]);
size(cp2)
cp2
cp2 = csvread('C:\Documents and Settings\asesha01\Desktop\JJRAL_50_3104to31204.csv', 1, 4, [1 4 50 4]);
cp2 = csvread('C:\Documents and Settings\asesha01\Desktop\JJRAL_50_3104to51004.csv', 1, 4, [1 4 50 4]);
cp2
for i=11:50, y2[i - 10] = 1/11*(cp2(i) + cp2(i - 1) + cp2(i-2)+cp2(i-3)+cp2(i-4)+cp2(i-5)+cp2(i-6)+cp2(i-7)+cp2(i-8)+cp2(i-9)+cp2(i-10)); end
for i=11:50, y2(i - 10) = 1/11*(cp2(i) + cp2(i - 1) + cp2(i-2)+cp2(i-3)+cp2(i-4)+cp2(i-5)+cp2(i-6)+cp2(i-7)+cp2(i-8)+cp2(i-9)+cp2(i-10)); end
y2
y2 = [zeros(1, 10) y2]
n = [1:50]
plot(n, cp2, 'o'); plot(n, y2, '*');
help plot
plot(n, cp2, 'o', n, y2, '*');
for i=11:50, y2(i - 10) = 1/11*(cp2(i) + cp2(i - 1) + cp2(i-2)+cp2(i-3)+cp2(i-4)+cp2(i-5)+cp2(i-6)+cp2(i-7)+cp2(i-8)+cp2(i-9)+cp2(i-10)); end
y2s = y2
plot(n, cp2, 'o', n, y2s, '*');
y2s
y2s = 0
for i=11:50, y2s(i - 10) = 1/11*(cp2(i) + cp2(i - 1) + cp2(i-2)+cp2(i-3)+cp2(i-4)+cp2(i-5)+cp2(i-6)+cp2(i-7)+cp2(i-8)+cp2(i-9)+cp2(i-10)); end
plot(n, cp2, 'o', n, y2s, '*');
y2 = 0
for i=11:50, y2(i - 10) = 1/11*(cp2(i) + cp2(i - 1) + cp2(i-2)+cp2(i-3)+cp2(i-4)+cp2(i-5)+cp2(i-6)+cp2(i-7)+cp2(i-8)+cp2(i-9)+cp2(i-10)); end
y2 = [zeros(1, 10) y2]
plot(n, cp2, 'o', n, y2, '*');
y2s = [y2s zeros(1, 10)]
plot(n, cp2, 'o', n, y2, '*');
plot(n, cp2, 'o', n, y2s, '*');
y2s = y2s(1:40)
y2s = [zeros(1, 5) y2s zeros(1, 5)]
plot(n, cp2, 'o', n, y2s, '*');
y2s = y2s(2, 50);
y2s = y2s(2:50);
y2s
y2s = [y2s zeros(1,1)]
plot(n, cp2, 'o', n, y2s, '*');
y2s = y2s(1:49);
y2s = [zeros(1,1) y2s]
plot(n, cp2, 'o', n, y2s, '*');
plot(n, cp2, 'o', n, y2, '*');
cp_ds = csvread('C:\Documents and Settings\asesha01\Desktop\JJRAL_50_3104to31204.csv', 1, 4, [1 4 50 4]);
cp_ds = csvread('C:\Documents and Settings\asesha01\Desktop\JJRAL_50_3104to51004.csv', 1, 4, [1 4 50 4]);
for i=1:50, xr(i) = cp_ds(i) - cp_ds(1) + ((cp_ds(1) - cp_ds(50))/49)*(i - 1); end;
xr
plot(n, xr)
function Xk = dft(x)
[N,M] = size(x);
if M ~= 1 % makes sure that x is a column vector
x = x';
N = M;
end
Xk = zeros(N, 1);
n = 0:N - 1;
for k = 0:N-1
Xk(k + 1) = exp(-j*2*pi*k*n/N)*x;
end
dft
Xr = dft(xr);
plot(n, Xr)
plot(n, abs(Xr))
plot(n, Xr)
plot(n, abs(Xr))
xr_hat(n) = (2/N)*((real(Xr(2))*cos((2*pi/N)*n)) - (imag(Xr(2))*sin((2*pi/N)*n)));
N = 50
xr_hat(n) = (2/N)*((real(Xr(2))*cos((2*pi/N)*n)) - (imag(Xr(2))*sin((2*pi/N)*n)));
plot(n,xr,'-',n,xr_hat,'--')
for i=1:50, cp_ds_hat(i) = xr_hat(i) + cp_ds(1) - ((cp_ds(1) - cp_ds(50))/49)*(i - 1); end;
plot(n,cp_ds,'-',n,cp_ds_hat,'--')
cp_ds_hat
34.2978 - 34.5227
help slope