北大心理學系.Psychtoolbox1

1、PsychtooboxLihan Chen2015.4.28MATLAB 圖像刺激 圖像處理 圖像基礎 圖像的變換與操縱 計算生成圖像圖像基礎圖像基礎像素、顏色、索引圖像圖像基礎像素 數(shù)字圖像 離散的坐標：像素 離散的強度：0255 (黑白)圖像基礎顏色 RGB顏色系統(tǒng) 每個像素的紅、綠、藍均由uint8表示 (0255) alpha通道：透明度 24位/32位真彩色imtool( photo.jpg);用MATLAB存取并顯示圖片讀取并顯示圖片 im = imread(tmp.png); size(im) imshow(im)寫入圖片 imwrite(im, tmp.jpg);a=ones(

2、400,40),zeros(400,40); %先生成一根白條和一根黑條a1=repmat(a,1,5); % 將這個組合復制5遍a2=xor(a1,a1);%將這個黑白條紋與自己的轉置做異或邏輯運算imshow(a2) %棋盤格就生成了實例：紅藍3D與雙眼競爭 im1 = imread(Koala.jpg); im2 = imread(Penguins.jpg);% 用第二幅圖替代第一幅圖的紅色通道 im1(:,:,1) = im2(:,:,1); % 把綠色通道清零 im1(:,:,2) = 0; imshow(im1)image與colormap image函數(shù)是MATLAB最基本的圖像

3、顯示函數(shù)，可以繪制索引圖像，即每個像素的值對應顏色查找表中的索引 colormap：定義圖像顯示用的顏色查找表 imagesc將數(shù)據(jù)scale后繪制成圖（例如繪制相關矩陣）image與imshow Imshow僅用于顯示由RGB或灰度值定義的圖像（image也可以） 無論是哪個函數(shù)，若圖像是以uint8表示的，則取值范圍為0255，若以double表示，則取值范圍是01image 繪制數(shù)據(jù) 246810123456789100.50.550.60.650.70.750.80.850.90.951imagesc & colormap 繪制相關矩陣a = rand(10);for i =

4、2:10a(:,i) = a(:,i-1)+rand(10,1)*0.5;endimagesc(corr(a)colorbarcolormap autumn圖像的變換與處理圖像的變換與處理調整曲線、窗口化像素亮度的變換操作 (PS中的“曲線”)實例：去除掃描時背面透過來的影 im = imread(book-000123.png); g = rgb2gray(im); d = double(g); y = uint8(d*1.7-85); imshow(y); imwrite(y, out.png);圖像的“窗口”操縱每個像素的亮度：乘以一個系數(shù)生成高斯窗口窗口化圖像 imagesc(GaWi

5、ndow) tmp = double(im).*GaWindow; Windowed = uint8(tmp); imshow(Windowed)鄰居處理鄰居處理高斯模糊、邊緣檢測空間濾波 定義一個矩陣(稱為濾波器filter) 將該矩陣覆蓋在每個像素點上 對所有被覆蓋的點 將原圖像矩陣和該矩陣的對應點相乘 求和變成新圖像的點 由鄰居決定新的像素點的強度高斯模糊 改變窗口大小，令x = y = 10, stdev = 5，生成一個“高斯窗口” 讓窗口矩陣的和等于1GaFilter = GaWindow/sum(sum(GaWindow); 用imfilter或filter2或conv2對原圖進

6、行濾波imf = imfilter(im, GaFilter);imshow(imf)高斯模糊的圖片邊緣檢測原始圖像原始圖像 grayedge(gray)imfilter(gray, fspecial(sobel)更多的濾波器 用fspecial生成二維的特殊濾波器 average 簡單平均 gaussian 高斯濾波 motion 模擬攝像機抖動 sobel 檢測邊緣算子 doc fspecial計算生成圖像計算生成圖像正弦光柵, Gabor Patch, 生成正弦光柵 X, Y=meshgrid(1:50,1:100); imshow(sin(X) imshow(sin(X+Y) imsh

7、ow(sin(X+Y)/2) imshow(sin(X+2*Y)/2) imshow(0.5+0.5*sin(3*X+4*Y)/6)完整代碼x,y=meshgrid(1:Boundx, 1:Boundy); % define the costant for the cosine wo=SAngle/180*pi; % conversion from deg to radPhRad=PhGrad/180*pi;% conversion from deg to radf=Fcxi/xysize(1); % frequency conversion ax = f*cos(wo); % convert

8、 the spatial frequency along yby = f*sin(wo); % convert the spatial frequency along x % Calculate the Sinusoid %(+1 is to obtatin all values being more than 1)ImOut=sin(2*pi*(ax*x+by*y)+PhRad)+1;生成Gabor% 先生成正弦光柵 X,Y=meshgrid(-50:50,-50:50); im = sin(X+2*Y)/3); imshow(im)% 生成窗口“蓋住”光柵 Ga = exp(-(X/20)

9、.2) - (Y/20).2); imshow(Ga.*im) Psychtoolbox函數(shù)介紹（圖像刺激的實現(xiàn)） SCREEN函數(shù)與視覺刺激函數(shù)與視覺刺激產(chǎn)生視覺刺激的三種方式 導入圖片(jpg,. tif,) 在PTB里畫圖 文本圖片（text figures)Screen函數(shù) Screen是PTB與視頻顯示有關的函數(shù)的集合 都有哪些子函數(shù)？ Command window輸入Screen Screen(函數(shù)字符串，參數(shù)列表） /Screen windowPtr, rect=Screen(Openwindow,windowPtrOrScr

10、eenNumber, , color, rect, pixelSize,numberOfBuffers ,stereomode , multisample , imagingmode , specialFlags , , clientRect);使用Screen函數(shù) 常用的子函數(shù)、使用順序及功能 查詢有多少屏幕 Screens 開窗口（創(chuàng)建主頁面） OpenWindow 繪制視覺刺激 畫線 DrawLine, 畫點 DrawDots, 貼圖 DrawTexture 頁面切換 Flip 關閉 CloseAll 附注 ShowCursor; HideCursor;Screen函數(shù)：Screens子

11、函數(shù) Screens子函數(shù)：查詢有多少屏幕 FAQ: 有多少屏幕？ 若有多個顯示器，一般主屏幕用來控制，編號最大的屏幕用來呈現(xiàn)刺激 此步可省略，可以直接開0號窗口 screens=Screen(Screens); screenNumber=max(screens);Screen函數(shù)： OpenWindow子函數(shù) OpenWindow用來打開一個窗口(“畫布”) w, rect = Screen(OpenWindow, screenNumber); 返回值 w: 該窗口的“句柄值”，類似于“指針” 以后對這個窗口進行任何操作都要用這個指針 rect: 該窗口的“矩形”，由左上、右下兩個點的座標定

12、義，比如0 0 1024 768 更多參數(shù) 背景顏色、窗口大??？ Screen OpenWindow?繪制視覺刺激 畫什么？ 線條？DrawLine 點？DrawDots 多邊形？FillPoly / FramePoly 橢圓？FillOval / FrameOval 文字？DrawText 照片？MakeTexture + DrawTexture 效果？ 位置？ 透明？Screen函數(shù)： DrawDots子函數(shù) 畫點：DrawDots Screen(DrawDots, w, xy, size, color, center,type); 用來畫一個或一堆點，它的參數(shù)們： w: 窗口的“指針”

13、xy: 兩行的行向量，每個點的x、y座標各占一行 size: 點的寬度，以像素計 color: (每個)點的顏色( R G B 或 R G B A ) center: 點座標的相對中心 type: 方點，圓點，還是抗鋸齒的圓點Screen函數(shù)： DrawLine子函數(shù) 畫線：DrawLine Screen(DrawLine, w, color, fromH, fromV, toH, toV, penWidth); 用來畫一條線，它的參數(shù)們： w: 窗口的“指針” color: 線的顏色 from/to H/V: 起點/終點的x/y座標 penWidth: 線條的粗細Screen函數(shù)： Make

14、Texture子函數(shù) Texture -由圖片轉換而來的“紋理”，便于快速貼圖 怎么把圖片轉成紋理？ image = imread(photo.jpg); textureIndex=Screen(MakeTexture, w, image); 返回值textureIndex就是這個紋理的“指針” DrawTexture子函數(shù) 把紋理（即照片）畫到屏幕上 Screen(DrawTexture, w, textureIndex); 把照片畫到哪里？ 源矩陣 src_rect=Screen (Rect,tex); 畫紋理的哪部分（默認全畫）？ 擺在中心 ctr_rect=CenterRect(src

15、_rect,wRect); 偏移 dst_rect=OffsetRect(dst_rect,-150,150); 畫到偏移完的目標上 Screen(DrawTexture, w, tex, src_rect, dst_rect);Screen函數(shù)： Flip子函數(shù) Psychtoolbox有兩個緩沖區(qū)(“畫板”) 屏幕上顯示的只是其中一面 顯示的同時可以在下面的畫板中繼續(xù)畫 需要更新的時候，直接切換，速度快 把畫板切換上去 馬上就翻：Screen(Flip,w); T時刻(或之后)再翻：vbl = Screen(Flip,w, T); 返回值vbl是真正刷屏時刻的時間戳(GetSecs)Scr

16、een函數(shù)： CloseAll子函數(shù) 為了避免程序運行時出錯，而無法執(zhí)行到CloseAll就卡在中間關不了全屏窗口 try screens=Screen(Screens); screenNumber=max(screens); w, rect = Screen(OpenWindow, screenNumber); Screen (CloseAll); catch Screen (CloseAll); rethrow(lasterror); end呈現(xiàn)透明刺激 Alpha通道 R G B A中第四個通道 0為完全不透明，255為完全透明 設置Blend Function Screen(Blend

17、Function, w, GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA); “去掉圖片中不要的部分”？Psychtoolbox-聲音刺激與反應錄入 聲音基礎與產(chǎn)生聲音聲音基礎與產(chǎn)生聲音采樣率、位數(shù)聲音基礎 音頻信號的表示 1或2行(列)采樣得到的行(列)向量 時間上的離散：采樣率 存儲上的離散：位數(shù) (在MATLAB中是double) 聲道：單聲道和立體聲 采樣率：每秒離散信號的采樣個數(shù) 44100Hz 位數(shù)：用多少位的數(shù)字表示信號對正弦波的4位采樣用MATLAB存取并播放聲音讀取并播放聲音 wav, sf, nbits = wavread(sound.wav);

18、size(wav) t = size(wav,1)/sf % 時長 = 采樣點個數(shù)/采樣率 wavplay(wav, sf)寫入聲音 wavwrite(wav, sf, nbits, sound.wav);播放聲音：wavplay 改變采樣率 改變播放時長，并變調 wavplay(wav,8000) wavplay(wav,80000) “異步”播放 指定async參數(shù)可以使聲音在“后臺”播放 否則wavplay會等待播完，然后才執(zhí)行下一句 利用async可以“疊加”兩個聲音 wavplay(wav, 44100, async)wavplay 與 sound 類似imshow和image，MA

19、TLAB中這兩個函數(shù)都可以播放一段聲音，用法類似 只有wavplay可以“異步”播放 soundsc將數(shù)據(jù)scale到-1 +1區(qū)間后播放（如果有意更改聲音的音量，不要用這個函數(shù)）生成噪音 生成一段單聲道噪音就是生成一列隨機數(shù) sf = 44100; % 采樣率 t = 1.5; % 持續(xù)時長 noise = rand(sf*t, 1); wavplay(noise, sf)plot(noise(1:200)02040608010012014016018020001生成純音 生成一段單聲道純音就是生成一列正弦波 sf = 44100; %

20、采樣率 t = 1.5; % 持續(xù)時長 f = 261.6; % 純音頻率(C4 = 261.6Hz) tmp = linspace(0, t, sf*t); % 計算正弦用的采樣矩陣 tone = sin(2 * pi * f * tmp); wavplay(tone, sf)plot(tone(1:200)020406080100120140160180200-1-0.8-0.6-0.4-0.60.81生成復合音疊加若干純音sf = 44100; % 采樣率t = 1.5; % 持續(xù)時長f = 261.6; % 基準頻率(C4 = 261.6Hz)tmp = linsp

21、ace(0, t, sf*t);tone1 = sin(2 * pi * f * tmp);tone2 = sin(2 * pi * 2*f * tmp);tone3 = sin(2 * pi * 4*f * tmp);harmonic = tone1 + tone2 + tone3; % 疊加并歸一化harmonic = harmonic / max(abs(harmonic);wavplay(harmonic, sf)plot(harmonic(1:200)020406080100120140160180200-1-0.8-0.6-0.4-0.60.81生成一段聲音拼接

22、若干純音 sf = 44100; % 采樣率 t = 0.5; % 持續(xù)時長 f1 = 261.6; f2 = 293.6; f3 = 329.6; % (C4, D4, E4) tmp = linspace(0, t, sf*t); do = sin(2 * pi * f1 * tmp); re = sin(2 * pi * f2 * tmp); mi = sin(2 * pi * f3 * tmp); silence = zeros(1, sf*t); % 靜音 wavplay(do re mi silence do re mi re do, sf)聲音處理的幾個一般函數(shù) 純音: Make

23、Beepbeep, samplingRate1= MakeBeep(freq,duration, samplingRate) freq-音調 duration-音頻長度（秒） samplingRate- 采樣頻率Beep-音頻數(shù)據(jù)返回值 samplingRate1- 返回采樣頻率 Beeper 格式 Beeper (frequency, loudness, duration) loudness “01” Snd(command, signal,rate,samplesize) command集成- Open, Play, Close, Wait, Quiet, DefaultRate, IsP

24、laying signal-單、雙通道 rate-播放頻率 samplesize-比特率(8或16） Snd(Play, sin(0:5000) Snd(Close) freq = 44100;tone(1,:) = 0.9 * MakeBeep(1000, 0.04, 44100); tone(2,:) = tone(1,:);toneb(1,:) = 0.9 * MakeBeep(500, 0.04, 44100); toneb(2,:) = toneb(1,:); Snd(Play, tone, freq, 16); Snd(Quiet)Snd(Close) 音頻數(shù)據(jù)播放: sound

25、Sound(y, Fs) Sound(y,Fs,bits) Fs-采樣頻率 bits- 比特數(shù) 音頻播放器：audioplayer player= audioplayer(Y,Fs,nBits) player= audioplayer(recorder, ID) 創(chuàng)建錄音機對象：audiorecorder recorder=audiorecorder(Fs, nBits, nChannels, ID)錄制一段聲音錄制一段聲音recorder=audiorecorder(44100, 8, 1, 0);disp(Start Recording);Recordblocking(recorder,

26、6); % record 6 secsdisp(End of Recording);play(recorder);audiodata=getaudiodata(recorder);plot(audiodata);產(chǎn)生其他波形 chirp: 產(chǎn)生信號的頻率隨時間的增長而變化 t=0:0.005:2; % 2 secs 5000 Hz SF y=chirp(t,220,2,440); wavplay(y,5000) sawtooth: 鋸齒波 fliplr, flipud: 把矩陣左右或上下顛倒 repmat: 重復矩陣頻譜圖“通過按鍵音還原手機號碼” 頻譜圖 t=0:0.0001:2; % 2

27、secs 10000Hz y=chirp(t,220,2,4400); wavplay(y,10000) spectrogram(y,128,120,128,1e5)立體聲、ITD和ILD 立體聲由兩列-1/+1之間的數(shù)表示 ITD(兩耳時間差)和ILD(兩耳強度差)是人進行聲源定位的兩個線索 HRTF(頭相關傳輸函數(shù))刻畫人耳接收空間中某一點傳來的聲音的特征 /projects/HRIRrepository聲音的操縱聲音的操縱包絡與濾波音量控制 數(shù)字音頻中，0dB是最大的音量，在此基礎上可以減小音量，如-10dB wav * 10(db/20

28、) 可將wav音量降低為db分貝，db0聲音的淡入和淡出 在聲音前后一小段加上“門控” gatedur = .01; gate = cos(linspace(pi, 2*pi, fs*gatedur); gate = (gate + 1) / 2; offsetgate = fliplr(gate); sustain = ones(1, (length(tone)-2*length(gate); envelope = gate, sustain, offsetgate; smoothed_tone = envelope .* tone;聲音的包絡(envelope) 獲取音頻設備信息: aud

29、iodevinfo devinfo= audiodevinfo Name= audiodevinfo(IO, ID); ID= audiodevinfo(IO,rate, bits, chans)Low latency, precisely timed sound with PsychPortAudio Based on PortAudio, a free, open-source, cross-platform audio library for realtime audio: http:/ Features: Multi-channel playback, recording and fu

30、ll-duplex feedback operation Low-latency sound output on MacOS/X, and MS-Windows (ASIO- required) Precisely timed sound onset on MacOS/X, and MS-Windows (ASIO requried)Low latency sound on Windows needs The standard sound subsystem of Microsoft Windows is not capable of reliably timed sound playback

31、 and capture with low latency If you want low latency sound with high timing precision on Microsoft Windows, youll need a sound card with native ASIO driver support and a special plugin for Psychtoolbox 采用ASIO(AudioStream Input Output)技術，可以減少系統(tǒng)對音頻流信號的延遲，增強聲卡硬件的處理能力。Delta 1010 聲卡聲卡Low latency,timed s

32、ound with PsychPortAudioSimple interface, allows for asynchronous sound output at ascheduled system time, e.g., time tvisualonset:1. Padevice = PsychPortAudio(Open, deviceid, mode, 2, 96000);2. Mysound = 0.9 * MakeBeep(1000, 0.1, 96000);3. PsychPortAudio(FillBuffer, Padevice, Mysound);4. PsychPortAu

33、dio(Start, Padevice, 5, tvisualonset);5. Visonset = Screen(Flip, window, tvisualonset 0.004);6. .whatever.7. Audioonset = PsychPortAudio(Stop, Padevice);8. PsychPortAudio(Close ,Padevice); Auto-selects settings for low latency, but overrides possible. High precision, low latency with standard hardwa

34、re on OS/X.PsychPortAudio-ASIO driverInitializePsychSound(1);freq = 96000;ISOI = 0.12; % 120 msStimdur=0.03; % auditory durationlatbias = (64 / freq); %hardware delay pahandle = PsychPortAudio(Open, ,2,freq);% Tell driver about hardwares inherent latencyprelat = PsychPortAudio(LatencyBias, pahandle,

35、 latbias);postlat = PsychPortAudio(LatencyBias, pahandle);%1: single sound; 2: two sounds; 3: multiple soundstone20 = 0.9 * MakeBeep(500,0.02,freq); %tone with 20 mstone(1,:) = 0.9 * MakeBeep(500, ISOI+Stimdur , freq); tone(2,1:length(tone20) = tone20; tone(2,end-length(tone20)+1:end) = tone20;for i

36、=0:2:6 range = 1+i*length(tone20): length(tone20)*(i+1); tone(3, range) = tone20;end PsychPortAudio(FillBuffer, pahandle, repmat(tone(1,:),2,1); PsychPortAudio(Start, pahandle, 1, 0); % PsychPortAudio(Stop, pahandle);使用前裝載聲使用前裝載聲卡驅動程序卡驅動程序List10.4Timingtry w=Screen(OpenWindow,0); WaitSecs(10); % 等待時

37、間，以秒為單位等待時間，以秒為單位 Screen(CloseAll);catch Screen(Close, w) rethrow(lasterror)endtry Screen(OpenWindow,0); t0= GetSecs; WaitSecs(5); t1= GetSecs; Screen(CloseAll ); t_elapsed = t1-t0; （反應時）（反應時）catch Screen(Close, w); rethrow(lasterror)endtry screenNum=0;res=1280 1024;clrdepth=32;wPtr,rect=Screen(Open

38、Window,screenNum,0,0 0 res(1) res(2), clrdepth);black=BlackIndex(wPtr);white=WhiteIndex(wPtr);Screen(FillRect,wPtr,black);priorityLevel = MaxPriority(wPtr);Priority(priorityLevel);Screen(wPtr, Flip);Priority(0); Screen(CloseAll);catch Screen(Close, wPtr) rethrow(lasterror)EndPriority() - Switch Matl

39、ab process to realtime-scheduling mode.Priority() - Switch Matlab process to realtime-scheduling mode. 優(yōu)先級優(yōu)先級Basic commands for system control and timingT = GetSecs （獲取時間值）（獲取時間值） Query time with microsecond resolution. Uses highest resolution system clock for measurement of time.WaitSecs(duration)

40、Wait for a specified amount of time duration. On MacOS/X: Accurate to 0.3 milliseconds on average. On MS-Windows: Accurate to KbDemoKbWaitdisp(Testing KbWait: hit any key. Just once.);startSecs = GetSecs;timeSecs = KbWait;keyIsDown, t, keyCode = KbCheck; str=KbName(keyCode), typed at time, num2str(t

41、imeSecs - startSecs), secondsdisp(str);Key Presses Collecting keypresses using pause and input Inputdlg promptParameters = Subject Name, Age, Gender (F or M?),Handedness (L or R);defaultParameters = , , R;Subinfo = inputdlg(promptParameters, Subject Info , 1, defaultParameters); Input resp_num= inpu

42、t (press a number key); resp_char = input (press a number key);disp(Using input command)resp=x;while resp=a & resp=bresp=input(press a or b . , s);endrespBear in mind that the command window needs to be in front for the key-press to be available to Matlab, so input doesnt work well if you are us

43、ing Psychtoolbox. In that case you will need to use GetChar, CharAvail, or KbCheck instead.% CharAvailWaitSecs(1)disp(Using CharAvail command)disp (Wait 2 sec to see if you pressed a key during that time);FlushEventsticwhile toc0 break; end end answer = find(buttons,1);MouseHideCursor;ShowCursor;Set

44、Mouse;GetClicks List10.14 List10.15List10.16dioIn=digitalio(parallel,LPT1);dioOut=digitalio(parallel,LPT1);addline(dioOut,0:7,out);addline(dioIn,11:12,in); getvalue(dioIn); if sum(getvalue(dioIn)=1 response=0; elseif sum(getvalue(dioIn)=0 response=1; elseif sum(getvalue(dioIn)=2 response=2; end Keyb

45、oard/Mouse Manipulations Adjustment with keyboard- Kingdom (1997, Perception)-see KeyboardMani.m Mller-Lyer illusion (Mller-Lyer, 1889) see mousemanipulation.mList10.17 Psychtoolbox實驗設計 內(nèi)容 實驗設計步驟和編程概要 刺激的隨機化 文件操作、結果保存 實驗設計實例 GUI圖形界面實驗步驟 實驗前關閉屏幕保護程序。 環(huán)境設置（圖形界面） 指導語的呈現(xiàn) 輸入被試個人資料（圖形界面） 進行實驗 記錄結果和數(shù)據(jù)初步整理刺

46、激的隨機化 reshape fullfact rand/seed randperm gentrials mod 被試內(nèi)、被試間和重復次數(shù)應用舉例 視覺搜索(Visual Search) 適應后效(Adaptation aftereffect) 協(xié)同運動(Coherent motion) 變化盲(Change Blindness) 快速序列視覺呈現(xiàn)(RSVP) 視聽整合：audiovisual Ternus, what you see is what you hear, bouncing illusion1、視覺搜索Set size: 4 levels(4,8,12,16)Conjunction

47、 vs. feature search;Target present vs. absent參考visualSearch, visualsearch2011文件夾，以及Ma et al. (2011)-Nature Neurosci程序框架和函數(shù)runVisualSearch.msettingsVisualSearch.mmakeSearchDisplay.m2、 適應后效-Adaptation to natural facial categories (Nature 2004：428, 558-561)程序框架-主函數(shù)和調用函數(shù)faceAdaptationMain.mgetParams.mdo

48、FaceAdapt.mdoFaceAdaptTrial.mopenFaceAdpt.mshowImage.mgetResponse.mSmallBlank.mfunction showImage(imName, wPtr, dur)img=imread(imName, JPG);textureIndex=Screen(MakeTexture, wPtr, double(img);Screen(DrawTexture, wPtr, textureIndex);Screen(wPtr, Flip);WaitSecs(dur)function response=getResponse(params)

49、FlushEvents;sound(params.beep);response=GetChar;function doFaceAdapt(params, wPtr, dur, time)lastImage=0; for repeat=1:(time/dur)whichImage=ceil(rand*params.numAdaptImages);while whichImage=lastImage whichImage=ceil(rand*params.numAdaptImages); endimName=params.adaptstr,num2str(whichImage);showImage

50、(imName, wPtr, dur);lastImage=whichImage;end完整的程序clear all; close all;cd(/Users/ariel/PTBTutorial/FaceExperiment/)addpath(pwd);addpath(morph);params = getParams;response = ;Screen(Preference, VisualDebuglevel, 3)wPtr rect = openFaceAdapt;doFaceAdapt (params, wPtr, params.preAdaptImageDuration, param

51、s.preAdaptTime);for trialCounter=1:params.nTrials responsetrialCounter = doFaceAdaptTrial (wPtr, params);if trialCounterparams.nTrials doFaceAdapt(params, wPtr, params.reAdaptImageDuration, params.reAdaptTime);endend Screen(CloseAll);ShowCursorfunction params=getParamsparams.preAdaptTime=10; params.preAdaptImageDuration=2; params.numAdaptImages=10; params.numTestImages=49; %There are 49 morph imagesparams.testImageDuration=0.5; %Presentation of the test image (in secs)params.reAdaptTime=6; %Top up adaptationpar