摘要：
1.在得到头动评估的数据基础上，可以汇出头动变化的折线图。运行本程序需要下载相应的实例数据headmove_data.txt头动是fMRI的一个重要问题，该程序画出了某个被试进行扫描时头动的具体变化情况%%%% 导入数据rp_ascan001=load('headmove_data.txt');figuresubplot(2,1,1)plot(rp_ascan001(:,1:3))xlabel('image');legend({'x translation','y translation','z translation&# 阅读全文
摘要：
1.多行注释单行注释是加%%{若干语句%}2.快捷键多行注释: 选中要注释的若干语句, 编辑器菜单Text-<Comment, 或者快捷键Ctrl+R取消注释: 选中要取消注释的语句, 编辑器菜单Text-<Uncomment, 或者快捷键Ctrl+T 3.利用逻辑判断，相当于“解释运行时”注释掉... 阅读全文
摘要：
本程序意在解释这样几个问题：完整版代码在本文的最后。1.实验的设计如何转换成设计矩阵？2.设计矩阵的每列表示一个刺激条件，如何确定它们？3.如何根据设计矩阵和每个体素的信号求得该体素对刺激的敏感性？程序详解：1.构造hrf hrf_small = [ 0 4 2 -1 0 ];figure(1);clf; plot(0:4,hrf_small,'o-'); grid on; xlabel('Time (in units of TRs, 4s long... 阅读全文
摘要：
这里包含了这样一个数据集：slice_data.mat。这个数据集中包含的mri数据是：64*64*25。共有25个slice。每个slice的分辨率是64*64。程序非常简短：load slice_data.mat %转载数据for i=1:25 %总共有25个切片 subplot(5,5,i) slice=scan1(:,:,i); imagesc(slice'); %这里对图像矩阵过了转置。不然头是歪的。end本程序讲解了：fMRI每1秒钟对大脑进行25层全脑扫描究竟发生了什么？得到的数据究竟什么样？如过要重复的做一件事：显示单层图像，用循环语句如何实... 阅读全文
摘要：
MATLAB支持工作区的保存。用户可以将工作区或工作区中的变量以文件的形式保存，以备在需要时再次导入。保存工作区可以通过菜单进行，也可以通过命令窗口进行。数据导出1. 保存整个工作区 选择File菜单中的Save Workspace As…命令，或者单击工作区浏览器工具栏中的Save，可以将工作区... 阅读全文
摘要：
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%clear all;clc;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% haemodynamic response function: BOLD signal has stereotyped shape every% time a stimulus hits it. We usually use the spm_hrf.m in SPM toolbox to% produce HRF. Here, we use the 29 vector 阅读全文
摘要：
使用过新浪博客的人都知道，新浪博客的首页有访问量统计功能，迁移到博客园之后发现博客园却没有这项功能，所幸博客园在后台管理的设置选项中有一个公告栏和设置页首页脚代码功能，使用起来非常灵活和方便。借此我们可以通过网络上提过的其他强大的插件来为获取更加详细的访问统计信息（博客园这难道是借鉴了传统unix中... 阅读全文