aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorEugeniy Mikhailov <evgmik@gmail.com>2015-08-12 13:36:57 -0400
committerEugeniy Mikhailov <evgmik@gmail.com>2015-08-12 13:36:57 -0400
commit5c95fbddb40ee9b0268014207a0a685a236338e2 (patch)
treeb59750c65a99043dd1a64e9dcd0416c5cfa42f14
parentddab7b9a6c2c54bca2ab76acf4871bb43db62cca (diff)
downloadlinux-matlab-gpib-5c95fbddb40ee9b0268014207a0a685a236338e2.tar.gz
linux-matlab-gpib-5c95fbddb40ee9b0268014207a0a685a236338e2.zip
removed obsolete file
-rw-r--r--HP8596E_take_data.m220
1 files changed, 0 insertions, 220 deletions
diff --git a/HP8596E_take_data.m b/HP8596E_take_data.m
deleted file mode 100644
index 510733b..0000000
--- a/HP8596E_take_data.m
+++ /dev/null
@@ -1,220 +0,0 @@
-% This script reads data from HP8596E spectrum analyzer and saves it to a file
-%
-% Mi Zhang (Revised from the code of Gleb Romanov gromanov@hellok.org)
-% 9/25/2013
-
-clear;
-
-
-
-% Define stuff
-board_index = 0;
-gpib_address = 18;
-data_prefix = 'S';
-data_path = 'Z:\Mi Squeezing\buffered cell\data\';
-run_number_file = 'Z:\Mi Squeezing\85Rb\autofile\runnum.dat';
-
-% Find a GPIB object.
-obj1 = instrfind('Type', 'gpib', 'BoardIndex', board_index, 'PrimaryAddress', gpib_address, 'Tag', '');
-
-% Create the GPIB object if it does not exist
-% otherwise use the object that was found.
-if isempty(obj1)
- %obj1 = gpib('NI', board_index, gpib_address);
- obj1 = lgpib('HP8596E');
-else
- fclose(obj1);
- obj1 = obj1(1);
-end
-
-% Adjust the buffers so the traces fit.
-% Do this before fopen(obj1);
-%obj1.InputBufferSize = 10240;
-%obj1.OutputBufferSize = 10240;
-
-% Connect to instrument object, obj1.
-%fopen(obj1);
-disp('--------------------')
-disp(horzcat('Connected to ',query(obj1, 'ID?')));
-
-% Communicating with instrument object, obj1.
-%
-% You can send commands using:
-% fprintf(obj1, '_command_');
-%
-% Or read stuff using:
-% _data_ = query(obj1, '_command_');
-disp('Reading data...');
-
-% We want our data in dB
-write(obj1, 'O3');
-
-% Read all 3 traces
-trA_string = query(obj1, 'TRA?');
-trB_string = query(obj1, 'TRB?');
-trC_string = query(obj1, 'TRC?');
-
-% Read various stuff
-device_string = query(obj1, 'ID?');
-freq_start_string = query(obj1, 'FA?');
-freq_stop_string = query(obj1, 'FB?');
-freq_center_string = query(obj1, 'CF?');
-freq_span_string = query(obj1, 'SP?');
-amplitude_units_string = query(obj1, 'AUNITS?');
-attenuation_string = query(obj1, 'AT?');
-ref_level_string = query(obj1, 'RL?');
-log_scale_string = query(obj1, 'LG?');
-rbw_string = query(obj1, 'RB?');
-vbw_string = query(obj1, 'VB?');
-sweep_time_string = query(obj1, 'ST?');
-
-% Disconnect from instrument object, obj1.
-clear(obj1);
-disp('Done');
-
-% Convert the traces from strings into vectors
-trA = sscanf(trA_string,'%f,');
-trB = sscanf(trB_string,'%f,');
-trC = sscanf(trC_string,'%f,');
-
-% Transpose the vectors
-trA = trA';
-trB = trB';
-trC = trC';
-
-% Create the frequency trace
-freq_start = sscanf(freq_start_string, '%f');
-freq_stop = sscanf(freq_stop_string, '%f');
-span = freq_stop - freq_start;
-freq = 0:length(trA)-1;
-freq = freq/max(freq);
-freq = freq_start + freq * span;
-
-
-% Plot stuff
-%
-% Open a window
-
-figure1 = figure(1);
-close(figure1);
-figure1 = figure(1);
-%
-% Create axes
-axes1 = axes('Parent',figure1,'YGrid','on','XGrid','on','FontSize',14);
-box(axes1,'on');
-hold(axes1,'all');
-%ylim([-2,12])
-%ylim([-85,-70])
-%
-% Create plot
-plot(freq/1e3,trA,'Color',[1 0 0],'Parent',axes1,'DisplayName','Trace A')
-plot(freq/1e3,trB,'Color',[0 0 0],'Parent',axes1,'DisplayName','Trace B')
-plot(freq/1e3,trC,'Color',[0 1 0],'Parent',axes1,'DisplayName','Trace C')
-%
-% Create xlabel
-xlabel('Detection frequency, kHz','FontSize',14);
-%
-% Create ylabel
-ylabel('Noise power, dB','FontSize',14);
-%
-% Show legend
-legend('show');
-
-
-% Plot with shot noise subtracted
-%
-% Open a window
-figure2 = figure(2);
-close(figure2);
-figure2 = figure(2);
-%
-% Create axes
-axes2 = axes('Parent',figure2,'YGrid','on','XGrid','on','FontSize',14);
-box(axes2,'on');
-hold(axes2,'all');
-%ylim([-2,12])
-%ylim([-85,-70])
-%
-% Create plot
-plot(freq/1e6,trA-trA,'Color',[1 0 0],'Parent',axes2,'DisplayName','Zero')
-plot(freq/1e6,trB-trA,'Color',[0 0 0],'Parent',axes2,'DisplayName','B - A')
-plot(freq/1e6,trC-trA,'Color',[0 1 0],'Parent',axes2,'DisplayName','C - A')
-%
-% Create xlabel
-xlabel('Detection frequency, MHz','FontSize',14);
-%
-% Create ylabel
-ylabel('Noise power, dB','FontSize',14);
-%
-% Show legend
-legend('show');
-%
-drawnow;
-
-
-
-% Open the file containing the run number and read it
-
-run_number_file_handle = fopen(run_number_file,'r');
-run_number = fscanf(run_number_file_handle,'%d');
-fclose(run_number_file_handle);
-% Increment it and write back
-run_number = run_number + 1;
-if run_number > 999
- run_number = 0;
-end;
-run_number_file_handle = fopen(run_number_file,'w');
-fprintf(run_number_file_handle, '%.3d', run_number);
-fclose(run_number_file_handle);
-
-
-
-
-% Get full path of the file to save
-
-save_to_file = horzcat(data_path,'S',num2str(run_number),'D',datestr(date,'yyyymmdd'),'.dat');
-
-% Write the data to a file
-disp(' ');
-disp(horzcat('Saving data to ',save_to_file));
-save_to_file_handle = fopen(save_to_file,'wt');
-fprintf(save_to_file_handle,'%s','# 14 (header lines)');
-fprintf(save_to_file_handle,'\n');
-fprintf(save_to_file_handle,'%s',horzcat('# ', datestr(clock)));
-fprintf(save_to_file_handle,'\n');
-fprintf(save_to_file_handle,'%s',horzcat('# Device:', ' ', device_string));
-fprintf(save_to_file_handle,'%s',horzcat('# Frequency center, Hz', ' ', freq_center_string));
-fprintf(save_to_file_handle,'%s',horzcat('# Frequency span, Hz', ' ', freq_span_string));
-fprintf(save_to_file_handle,'%s',horzcat('# Frequency start, Hz', ' ', freq_start_string));
-fprintf(save_to_file_handle,'%s',horzcat('# Frequency stop, Hz', ' ', freq_stop_string));
-fprintf(save_to_file_handle,'%s',horzcat('# Amplitude units ', amplitude_units_string));
-fprintf(save_to_file_handle,'%s',horzcat('# Attenuation ', attenuation_string));
-fprintf(save_to_file_handle,'%s',horzcat('# Reference level ', ref_level_string));
-fprintf(save_to_file_handle,'%s',horzcat('# Log scale ', log_scale_string));
-fprintf(save_to_file_handle,'%s',horzcat('# Resolution bandwidth, Hz', ' ', rbw_string));
-fprintf(save_to_file_handle,'%s',horzcat('# Video bandwidth, Hz', ' ', vbw_string));
-fprintf(save_to_file_handle,'%s',horzcat('# Sweep time, seconds', ' ', sweep_time_string));
-%
-% An example of how you shouldn't write data:
-%
-% for i = 1:length(freq)
-% fprintf(save_to_file_handle,'%s',sprintf('%f',freq(i)));
-% fprintf(save_to_file_handle,'\t');
-% fprintf(save_to_file_handle,'%s',sprintf('%f',trA(i)));
-% fprintf(save_to_file_handle,'\t');
-% fprintf(save_to_file_handle,'%s',sprintf('%f',trB(i)));
-% fprintf(save_to_file_handle,'\n');
-% end;
-%
-% This is much faster:
-data = [freq; trA; trB; trC;];
-fprintf(save_to_file_handle,'%f\t%f\t%f\t%f\r\n',data);
-% Close the file
-fclose(save_to_file_handle);
-disp('Done');
-
-% Close all opened files
-fclose('all');
-
-% Bring focus back to the command window
-commandwindow;