{[ promptMessage ]}

Bookmark it

{[ promptMessage ]}

LAB REPORT

LAB REPORT - 2.1 function[data count =...

This preview shows pages 1–3. Sign up to view the full content.

2.1 function [data, count] = readbuoydata(datafile) fid = fopen(datafile, 'r' ); %Open datafile tline = fgetl(fid); %Read next line of file tline = fgetl(fid); %Read next line of file [A,count] = fscanf(fid, '%d %d %d %d %d %f %f %d %f %f' ,[10 inf]); %Read file and convert into array A. Count is te number of successful reads data.date = datenum([A(1:5,:); zeros(1,size(A,2))]')'; %Convert date/time in specified vectors data.Hs = A(6,:); % Assign 6 th column of array as Hs data.Tp = A(7,:); % Assign 7 th column of array as Hs data.Dp = A(8,:); % Assign 8 th column of array as Hs data.Ta = A(9,:); % Assign 9 th column of array as Hs data.SST = A(10,:); % Assign 10 th column of array as Hs fclose(fid); 2.2 subplot(2,1,1); plot(data.date,data.Hs); datetick( 'x' ,6) xlabel( 'Time' ) ylabel( 'Period' ) subplot(2,1,2); plot(data.date,data.Tp); xlabel( 'Time' ) ylabel( 'Period' ) datetick( 'x' ,6)

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
2.3 function y=mafilt(N,x) for n= N:length(X), y(n)=sum(X(n-2:n))/N; end end 2.4 avgdata=mafilt(5,data.Tp); subplot(3,1,1) plot(data.Tp, 'color' ,[0.8 0.8 0.8], 'LineWidth' ,2); hold on ; plot(avgdata, 'k' ); hold off ; avgdata=mafilt(21,data.Tp); subplot(3,1,2) plot(data.Tp, 'color' ,[0.8 0.8 0.8], 'LineWidth' ,2); hold on ;
This is the end of the preview. Sign up to access the rest of the document.

{[ snackBarMessage ]}