MATLAB的答案

为什么我的图是空白的?

2个视图(过去30天)
显示旧的评论
艾丽西亚梅森
艾丽西亚梅森 2021年7月14日
  • 0
评论道: 明星黾 2021年7月29日
我试图用这个脚本种子sse绘制一些MEG数据的频谱。m，最初是由奈丁为他2016年研究语言结构的皮质追踪而写的: https://www.nature.com/articles/nn.4186
清晰的
负载(“R1167_test.mat”）
为实验室= 1:尺寸(数据,3)
x (:,:, trl) =重新取样(数据(:,:,实验室),1,1);结束
x =单(x);
保存r1167 - 1000 hz_test x
% %
清晰;
负载(“r1167 - 1000 hz_test.mat”）
fs = 1000;
[b] = fir1(200年,[6]0.5 / (fs / 2));一个= 1;
清晰的xf
为实验室= 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, fromdss,比、压水式反应堆]= dss0 (c0, c1, keep1 keep2);
y =折叠(展开(xf) * todss (: 1:10) * fromdss(1:10,:),大小(xf, 1));
(sa),某人,f) = plot_itc_avg_spectrum (y, fs, [0 5);
次要情节211甘氨胆酸;集(,“xtick”，[1/4 / 1/2 1 2 4]/0.32);网格
持有在;情节(f,意味着(sa), 2) * 3,“k”，“线宽”, 2)
次要情节212甘氨胆酸;集(,“xtick”，[1/4 / 1/2 1 2 4]/0.32);网格
持有在;情节(f,意味着(某人,2)* 3,“k”，“线宽”, 2)
这使用了函数plot_itc_avg_spectrum，如下所示:
函数[itc_sprectrum avg_spectrum f] = plot_itc_avg_spectrum (fs, eeg_data display_frequency_range, f_label)
计算并绘制神经反应谱。
%的输入:
% eeg_data(时间*通道*试用)
% fs(采样率，以Hz为单位)
% display_frequency_range(将要绘制的频率范围)
% f_label(可选，频率被标记)
%输出:
% itc_谱(频率*通道，试验间相位相干谱R^2)，
% avg_spectrum(频率*通道，试验平均响应功率谱)
% f(频率标签)
%的例子:
200年% plot_itc_avg_spectrum (eeg_data [0.5, 4.5], 1:4);
％
丁奈，张文，2016。
% ding_nai@zju.edu.cn
f = 1:尺寸(eeg_data 1); f = f - 1; f = f /大小(eeg_data 1); f = f * fs;
阿兹= 0(大小(eeg_data, 1),大小(eeg_data 3)大小(eeg_data, 2));
为ch = 1:尺寸(eeg_data, 2)
阿兹(:,:,ch) =角(fft(挤压(eeg_data (:, ch,:))));
结束
itc_sprectrum =挤压(pcoh3 (az));
avg_spectrum = abs (fft(平均(eeg_data 3)));
图;
次要情节(211);
情节(f, itc_sprectrum);
xlim (display_frequency_range);
包含(的频率(赫兹)）
ylabel (“inter-trial相位相干性的）
试一试
集(gca),“xtick”, f_label);结束
标题(“itc \ _sprectrum”）;
次要情节(212);
情节(f, avg_spectrum);
xlim (display_frequency_range);
标题(“avg \ _spectrum”）;
包含(的频率(赫兹)）
ylabel (“权力(a.u)。”）
试一试
集(gca),“xtick”, f_label);结束
结束
函数r = pcoh3 (ag)
c = cos (ag);
s =罪(ag);
r = (c, 2)意思。^ 2 +的意思(年代,2)^ 2;
结束
当我单独运行plot_itc_avg_spectrum时，使用一个形式为时间x通道x试验的MatLab数组作为输入，它工作得很好，并为我提供了一个显示数据频谱的图，具有正确标记的轴。
然而，当我运行整个seedssE脚本(带有用于执行去噪源分离(DSS)的各种功能)时，脚本运行时没有任何明显的错误，但产生了一个空白的图。
这两个图都是附加的。
我怎样才能确保运行整个seedssE脚本的结果不是空白的，并且看起来像单独运行plot_itc_avg_spectrum的结果?
5个评论
显示 4年长的评论
明星黾
明星黾 2021年7月29日
我的荣幸!
我还是不知道是什么 “折” 而且 “展开” 在代码中做什么，或者它们是什么。
．

登录发表评论。

登录来回答这个问题。

答案(0)

类别

标签

s manbetx 845


释放

R2019a

社区寻宝

在MATLAB中央找到宝藏，并发现社区可以如何帮助你!

开始狩猎!