clear all;[wptr,wrect]=Screen('OpenWindow',1,0,[50,50,1500,720]);%创建一个窗口DrawFormattedText(wptr,'Dong 2020nniHarbour.','center','center',255);%文字内容和格式WaitSecs(2);%显示2sScreen('Flip',wptr);%显示到屏幕上 Screen('CloseAll');%关闭窗口

clc;InitializePsychSound(1);[waveform,Fs] = audioread('di.wav');numChannels = size(waveform, 2);pahandle = PsychPortAudio('Open', 2, [], 1, Fs, numChannels);PsychPortAudio('FillBuffer', pahandle, waveform' );startTime = PsychPortAudio('Start',pahandle,1);WaitSecs(1);%做一个延时，否则下一句代码执行，会立即停止声音的播放PsychPortAudio('Close');%可有可无，最好有

%% 画矩形PsychDefaultSetup(2);[window, windowRect] = PsychImaging('OpenWindow', 0,0,[0,0,1920,1200]);%矩形参数设置[xCenter, yCenter] = RectCenter(windowRect);baseRect = [0 0 200 200];centeredRect = CenterRectOnPointd(baseRect, xCenter, yCenter);rectColor = [1 0 0];Screen('FillRect', window, rectColor, centeredRect);%在屏幕上显示以上内容Screen('Flip', window);%按键终止程序KbStrokeWait;sca;

%% 根据按下的上下箭头反馈对应的文字%% keyboard arrow%初始化ptbPsychDefaultSetup(2);[wptr, windowRect] = PsychImaging('OpenWindow', 0,0,[0,0,1920,1200]);%设置按键upKey = KbName('UpArrow');downKey = KbName('DownArrow');%操作提示feedbacktext='up:↑             down:↓nnwhen screen is blacknselect in 2 seconds';justSaySomething(wptr, feedbacktext);WaitSecs(2);Screen('Flip',wptr);%两秒钟时间做出按键选择t=GetSecs;while GetSecs

% MATLAB API V1.16.00%说明：API使用手册提供的demo中是使用while函数，当采集到指定长度的脑电信号之后停止,因为matlab本身是一行一行代码完成的，那么在实际和范式结合的时候，不能满足同时呈现范式以及采集当前刺激下脑电的记录，因此使用timer函数进行信号的采集clc;%初始化PTB设置，因为懒得写延时的代码，用的PTB的waitsecs函数PsychDefaultSetup(2);%因为用到timer，初始化全局变量global gds_interfaceglobal eegdatagds_interface = gtecDeviceInterface();eegdata = [];%初始化空矩阵用于存放eegdata%端口设置gds_interface.IPAddressHost = '127.0.0.1';gds_interface.IPAddressLocal = '127.0.0.1';gds_interface.LocalPort = 50224;gds_interface.HostPort = 50223;%连接设备connected_devices = gds_interface.GetConnectedDevices();%确认设备为USBampgusbamp_config = gUSBampDeviceConfiguration();gusbamp_config.Name = connected_devices(1,1).Name;gds_interface.DeviceConfigurations = gusbamp_config;%获取可用通道数目available_channels = gds_interface.GetAvailableChannels();%设置采样率和扫描数gusbamp_config.SamplingRate = 256;gusbamp_config.NumberOfScans = 8;%通道设置，滤波等for i=1:size(gusbamp_config.Channels,2)    if (available_channels(1,i))    	gusbamp_config.Channels(1,i).Available = true;        gusbamp_config.Channels(1,i).Acquire = true;        % do not use filters        gusbamp_config.Channels(1,i).BandpassFilterIndex = 49;        gusbamp_config.Channels(1,i).NotchFilterIndex =2;        % do not use a bipolar channe        gusbamp_config.Channels(1,i).BipolarChannel = 0;    endend%应用以上设置gds_interface.DeviceConfigurations = gusbamp_config;gds_interface.SetConfiguration();%设置定时器，用于连续读取缓冲区的脑电信号t=timer('Name','geteegdata','TimerFcn',{@geteegdata},'Period',0.3,'ExecutionMode','fixedSpacing');%开始采集（缓冲区开始读取）gds_interface.StartDataAcquisition();%计时器开始工作，记录脑电start(t);WaitSecs(4);%4s%停止计时器stop(t);%缓冲区停止获取数据gds_interface.StopDataAcquisition();%删除计时器delete(t);%解除设备占用delete(gds_interface);clear gds_interface;clear gusbamp_config;%save eegdata; %保存脑电数据以离线处理%画图看一下plot(eegdata(:,1));%通道1的xlabel('channel 1');ylabel('a');%计时器回调函数function geteegdata(~,~)         global gds_interface         global eegdata         global data         [~ , data] = gds_interface.GetData();         eegdata=[eegdata;data];end