% 产生新的样本点 new_idx = zeros((1-p0)*N, 1); %-------------- for i = 1:M for j = 1:(1-p0)*N idx = chains(i, j); new_idx(j) = [new_idx(j); idx]; end new_idx = cell( (1-p0)*N,1); end无法执行赋值,因为左侧和右侧的元素数目不同。
时间: 2023-12-17 17:05:23 浏览: 30
这段代码中存在两个问题:
1. 在执行 `new_idx(j) = [new_idx(j); idx];` 时,`new_idx(j)` 是一个标量,而右侧 `[new_idx(j); idx]` 是一个向量。这样无法将一个向量赋值给一个标量,需要将 `new_idx(j)` 改为一个向量,比如使用 `new_idx{j} = [new_idx{j}; idx];`。
2. 在循环中,每次都将 `new_idx` 初始化为一个空的 `cell` 数组,这样会导致前面循环中已经添加的元素丢失。应该在循环外部先初始化 `new_idx`,然后在循环内部只进行向 `cell` 数组的添加。
修改后的代码如下:
```
new_idx = cell((1-p0)*N, 1);
for j = 1:(1-p0)*N
for i = 1:M
idx = chains(i, j);
new_idx{j} = [new_idx{j}; idx];
end
end
```
相关问题
new_idx = zeros((1-p0)*N, 1); %----------------以下有问题--------------- for i = 1:M for j = 1:(1-p0)*N idx = chains(i, j); new_idx(j) = [new_idx(j); idx]; end new_idx = cell( (1-p0)*N,1); end 无法执行赋值,因为左侧和右侧的元素数目不同。
这段代码有一些错误,可能是由于变量类型不匹配导致的。你可以尝试以下修改:
```
new_idx = zeros((1-p0)*N, 1);
for j = 1:(1-p0)*N
idx = chains(1, j);
new_idx(j) = idx;
end
new_idx = num2cell(new_idx);
```
这里我假设 `chains` 是一个矩阵,第一维度有 `M` 个元素,第二维度有 `(1-p0)*N` 个元素。我把 `for` 循环的内部改成了一个循环,因为 `chains` 矩阵的第一维度似乎没有用到。我还把 `new_idx` 初始化为一个数值数组,并且在循环中直接把 `idx` 赋值给 `new_idx(j)`。最后,我把 `new_idx` 转换成了一个 cell 数组,因为你原来的代码中用到了 cell 数组的语法。
请注意,这个修改可能还需要根据你的实际情况进行调整。如果你需要更详细的帮助,请提供更多的代码和上下文信息,这样我才能更好地理解你的问题。
帮我修改一下这个matlab代码的错误% 读取语音信号 [x, fs] = audioread("C:\Users\ASUS\Desktop\20230607_210020.wav"); % 设置参数 win_len = 320; % 窗长 overlap = 0.5; % 帧重叠比例 win_type = 'hamming'; % 窗类型,可以选择'rectangle'或'hamming' % 分帧加窗win = window(win_type, win_len); n_overlap = win_len * overlap; n_frame = floor((length(x) - win_len) / n_overlap) + 1; frames = zeros(win_len, n_frame); for i = 1:n_frame idx_start = (i-1) * n_overlap + 1; idx_end = idx_start + win_len - 1; frames(:,i) = x(idx_start:idx_end) .* win; end % 计算短时自相关函数 time_delay = -win_len+1:win_len-1; R = zeros(win_len*2-1, n_frame); for i = 1:n_frame R(:,i) = xcorr(frames(:,i), 'biased'); end % 归一化处理 for i = 1:n_frame R(:,i) = R(:,i) / R(win_len,i); end % 绘制时域波形 t = linspace(0, length(x)/fs, length(x)); figure;plot(t, x); title('时域波形'); xlabel('时间(秒)'); ylabel('幅值'); % 绘制短时自相关函数 figure; imagesc(1:n_frame, time_delay, R); axis('xy') title('短时自相关函数'); xlabel('帧序号'); ylabel('时间延迟');
这段 MATLAB 代码存在以下问题:
1. 反斜杠字符 '\' 在 MATLAB 中是转义字符,需要用两个反斜杠 '\\' 来表示一个反斜杠,否则会导致路径错误。
2. 函数 window() 和 xcorr() 的第一个参数应该是窗长度 win_len,而不是数据 frames(:,i)。
3. imagesc() 函数缺少逗号分隔符。
下面是修改后的代码:
```
% 读取语音信号
[x, fs] = audioread("C:\\Users\\ASUS\\Desktop\\20230607_210020.wav");
% 设置参数
win_len = 320; % 窗长
overlap = 0.5; % 帧重叠比例
win_type = 'hamming'; % 窗类型,可以选择'rectangle'或'hamming'
% 分帧加窗
win = window(win_type, win_len);
n_overlap = win_len * overlap;
n_frame = floor((length(x) - win_len) / n_overlap) + 1;
frames = zeros(win_len, n_frame);
for i = 1:n_frame
idx_start = (i-1) * n_overlap + 1;
idx_end = idx_start + win_len - 1;
frames(:,i) = x(idx_start:idx_end) .* win;
end
% 计算短时自相关函数
time_delay = -win_len+1:win_len-1;
R = zeros(win_len*2-1, n_frame);
for i = 1:n_frame
R(:,i) = xcorr(frames(:,i), win_len-1, 'biased');
end
% 归一化处理
for i = 1:n_frame
R(:,i) = R(:,i) / R(win_len,i);
end
% 绘制时域波形
t = linspace(0, length(x)/fs, length(x));
figure; plot(t, x);
title('时域波形');
xlabel('时间(秒)');
ylabel('幅值');
% 绘制短时自相关函数
figure; imagesc(1:n_frame, time_delay, R);
axis('xy');
title('短时自相关函数');
xlabel('帧序号');
ylabel('时间延迟');
```
希望能对您有所帮助!