MATLAB答案

为什么我的情节出现空白？

1视图（过去30天）
显示较旧的评论
艾丽西娅·梅森
艾丽西娅·梅森 2021年7月14日
  • 0.
评论： 星际漫游者 2021年7月29日
我试图使用脚本seedssE.m为一些MEG数据绘制频谱，该脚本最初由Nai Ding为2016年的语言结构皮层跟踪研究撰写： https://www.nature.com/articles/nn.4186
清楚的
加载（'r1167_test.mat'的）
为了trl=1：大小（数据，3）
x（：，：，trl）=重组（数据（：，：，trl），1,1）;结尾
x=单个（x）；
拯救R1167-1000Hz_test X.
%%
清除;
加载（'r1167-1000hz_test.mat'的）
FS = 1000;
[b] =fir1（200[0.56]/（fs/2））；a=1；
清楚的XF.
为了trl = 1：尺寸（x，3）
xf（：，：，trl）=过滤器（b，a，x（：，：，trl））；
结尾
xf=xf（fs+[1:0.32*fs*4*10]，：，：）；
x0=展开（xf）；
x1=平均值（xf，3）；
c0 = x0'* x0;
c1=x1'*x1；
KEEP2 = 10. ^ - 12; KEEP1 = [];
[TODSS，FROMDS，比率，PWR] = DSS0（C0，C1，KEEP1，KEEP2）;
y=折叠（展开（xf）*托兹（：，1:10）*从DSS（1:10，：），大小（xf，1））；
[sa，sb，f]=plot_-itc_-avg_谱（y，fs[05]）；
子地块211;set（gca，“克斯蒂克”，[1/4 1/2 1 2 4] / 0.32）;网格
抓住在…上;图（F，平均值（SA，2）* 3，“k”那“线宽”,2)
子地块212;set（gca，“克斯蒂克”，[1/4 1/2 1 2 4] / 0.32）;网格
抓住在…上;图（f，平均值（sb，2）*3，“k”那“线宽”,2)
这使得使用函数plot_itc_avg_spectum，完整如下所示：
作用[ITC_SPRECTRUM，AVG_SPECTRUM，F] = PLOT_ITC_AVG_SPECTRUM（EEG_DATA，FS，DISPLAY_FREQUENCY_RANGE，F_LABEL）
％计算和绘制神经响应谱
%输入：
%eeg_数据（时间*通道*试验）
％FS（Hz中的抽样率）
％display_frequency_range（将绘制的频率范围）
%f_标签（可选，正在标记的频率）
％ 输出：
%itc_谱（频率*通道，试验间相位相干谱R^2），
％avg_spectrum（频率*通道，在试验中平均响应的功率谱）
%f（频率标签）
%例如：
％plot_itc_avg_spectrum（eeg_data，200，[0.5,4.5]，1：4）;
％
％Nai Ding＆Wen Zhang，2016
％ding_nai@zju.edu.cn.
f=1：大小（脑电图数据，1）；f=f-1；f=f/大小（脑电图数据，1）；f=f*fs；
AZ = zeros（大小（eeg_data，1），size（eeg_data，3），size（eeg_data，2））;
为了ch=1：大小（eeg_数据，2）
az（：，：，ch）=角度（fft（挤压（eeg_数据（：，ch，：））））；
结尾
itc_sprectrum=挤压（pcoh3（az））；
avg_谱=abs（fft（平均值（eeg_数据，3））；
图形
子图（211）;
地块（f、itc_光谱）；
xlim（display_frequency_range）;
xlabel('频率（Hz）'的）
伊拉贝尔('试验相干阶段连贯'的）
尝试
set（gca，“克斯蒂克”，f_标签）；结尾
头衔('ITC \ _sprectrum');
子图（212）;
绘图（F，AVG_Spectrum）;
xlim（display_frequency_range）;
头衔(“平均频谱”);
xlabel('频率（Hz）'的）
伊拉贝尔(“电源（a.u.）”的）
尝试
set（gca，“克斯蒂克”，f_标签）；结尾
结尾
作用r = pcoh3（ag）
c = cos（ag）;
s=sin（ag）；
r=平均值（c，2）。^2+平均值（s，2）。^2；
结尾
当我自己运行plot_itc_avg_频谱时，使用MatLab阵列的form time x channel x trials作为输入，它工作正常，并给出一个显示数据频谱的曲线图，该数据具有正确标记的轴。
但是，当我运行整个seedssE脚本时，使用设计用于执行去噪源分离（DSS）的各种函数，脚本运行时没有任何明显错误，但会生成一个空白图。
这两幅图均附呈。
如何确保从运行整个种子脚本的绘图不空白，看起来像运行Plot_itc_avg_spectrum自己的曲线？
5条点评
显示＆nbsp4旧的评论
星际漫游者
星际漫游者 2021年7月29日
我的荣幸！
我还是不知道是什么 “折叠” 和 “展开” 在代码中执行，甚至是它们是什么。
。

登录评论。

登录来回答这个问题。

答案（0）

类别

标签

s manbetx 845


释放

R2019a

社区寻宝

找到Matlab Central中的宝藏，并发现社区如何帮助您！

开始打猎吧！