Введение

В этом примере у вас есть схема R-C, состоящая из 1 К$\Omega$ резистор последовательно с 0.1 $\mu$F конденсатор. Схема R-C присоединена к устройству ADALM1000 с Каналом устройства, обеспечивающего стимул напряжения, состоящий из сигнала щебета. Выход Канала A соединяется с резистором, и земля соединяется с конденсатором. Образуйте канал B устройства используется, чтобы измерить напряжение через конденсатор. Следующая принципиальная схема показывает настройку измерения.

Можно использовать стимул и измеренные формы волны, чтобы оценить передаточную функцию схемы R-C и сравнить измеренный отклик с теоретическим ответом схемы R-C.

Узнайте устройства, подключенные с вашей системой

daqlist("adi")

ans=1×4 table
    DeviceID              Description                 Model              DeviceInfo       
    ________    _______________________________    ___________    ________________________

     "SMU1"     "Analog Devices Inc. ADALM1000"    "ADALM1000"    [1×1 daq.adi.DeviceInfo]

Создайте DataAcquisition для устройства ADALM1000

ADIDaq = daq("adi")

ADIDaq = 
DataAcquisition using Analog Devices Inc. hardware:

                     Running: 0
                        Rate: 100000
           NumScansAvailable: 0
            NumScansAcquired: 0
              NumScansQueued: 0
    NumScansOutputByHardware: 0
                   RateLimit: []

Show channels
Show properties and methods

Добавьте каналы источника и измерения напряжения

Устройство ADALM1000 способно к определению источника и измерению напряжения одновременно на отдельных каналах. Настройте устройство в этом режиме.

К исходному напряжению добавьте канал аналогового выхода с устройством ID SMU1 и ID канала A, и установите его тип измерения на Напряжение.

addoutput(ADIDaq,'smu1','a','Voltage');

Чтобы измерить напряжение, добавьте канал аналогового входа с устройством ID SMU1 и ID канала B, и установите его тип измерения на Напряжение.

addinput(ADIDaq,'smu1','b','Voltage');

Подтвердите настройку каналов.

ADIDaq.Channels

ans=1×2 object
    Index    Type    Device    Channel      Measurement Type             Range            Name  
    _____    ____    ______    _______    _____________________    _________________    ________

      1      "ao"    "SMU1"      "A"      "Voltage (SingleEnd)"    "0 to +5.0 Volts"    "SMU1_A"
      2      "ai"    "SMU1"      "B"      "Voltage (SingleEnd)"    "0 to +5.0 Volts"    "SMU1_B"

Задайте и визуализируйте форму волны щебета для стимулирования схемы

Используйте форму волны щебета амплитуды на 1 вольт, располагающейся в частоте от 20 Гц до 20 кГц для стимулирования схемы. Щебет происходит в период 1 секунды.

Fs = ADIDaq.Rate;
T = 0:1/Fs:1;
ExcitationSignal = chirp(T,20,1,20e3,'linear');

Добавьте смещение DC 2 В, чтобы гарантировать, что выходное напряжение устройства всегда выше 0 В.

Offset = 2;
ExcitationSignal = ExcitationSignal + Offset;

Визуализируйте сигнал стимула во временном интервале.

figure(1)
plot(T, ExcitationSignal,'Blue')
xlim([0 0.15])
xlabel('Time (s)')
ylabel('Magnitude (V)')
title('Time domain plot of stimulus signal')

Визуализируйте сигнал стимула в частотном диапазоне.

figure(2)
spectrogram(ExcitationSignal,1024,1000,1024,Fs,'yaxis')
title('Frequency domain view of stimulus signal')

Стимулируйте схему и одновременно измерьте частотную характеристику

Сгенерируйте устройство, выход и мера измеряют напряжение через конденсатор одновременно на другом канале.

MeasuredTable = readwrite(ADIDaq,ExcitationSignal');
MeasuredSignal = MeasuredTable{:,1};

Постройте стимул и измеренные сигналы

figure(1)
hold on;
plot(T,MeasuredSignal,'Red');
xlim([0 0.15])
ylim([1 3])
title('Time domain plot of stimulus and measured signal')
legend('Excitation Signal','Measured Signal')

figure(3)
spectrogram(MeasuredSignal,1024,1000,1024,Fs,'yaxis')
title('Frequency domain view of the the measured signal')

Вычислите передаточную функцию схемы

Сравните измеренный сигнал и сигнал стимула вычислить передаточную функцию схемы R-C и построить ответ величины.

Удалите смещение DC перед обработкой.

MeasuredSignal = MeasuredSignal - mean(MeasuredSignal);
ExcitationSignal = ExcitationSignal - Offset;
[TFxy,Freq] = tfestimate(ExcitationSignal,MeasuredSignal,[],[],[],Fs);
Mag = abs(TFxy);

Сравните предполагаемую передаточную функцию с теоретическим ответом величины.

R = 1000;    % Resistance (ohms)
C = 0.1e-6;  % Capacitance (farads)
TFMagTheory = abs(1./(1 + (1i*2*pi*Freq*C*R)));

figure(4);
semilogy(Freq,TFMagTheory,Freq,Mag);
xlim([0 20e3])
xlabel('Frequency (Hz)')
ylim([0.05 1.1])
ylabel('Magnitude')
grid on
legend('Theoretical frequency response','Measured frequency response')
title({'Magnitude response of the theoretical'; 'and estimated transfer functions'});

Очистите DataAcquisition и фигуры

clear ADIDaq
close all