Question

求matlab高手解答

定义函数
function
err_rate=PrimePe(snr_in_dB)

E=1;

snr=10^(snr_in_dB/10);

sgma=sqrt((E/snr)/2);

%桢个数

N=10;

%信息序列

frame_length=260;

signal=randint(1,N*frame_length);

%不经过卷积编码

%经信道加高斯燥声

u=2*signal-1;

u=sign(u+sgma*randn(size(u))+eps);

u=(u+1)/2;

%误码率

err_rate=(length(find(signal~=u)))/length(signal);

主函数
snr_in_db=5;
err_rate=PrimePe(snr_in_dB);
plot(snr_in_db,err_rate);

运行后出现

??? Error using ==> mpower
Inputs must be a scalar and a square matrix.

Error in ==> PrimePe at 3
snr=10^(snr_in_dB/10);

求高手解答
定义函数都改为点乘和点除了
已改为u=sign(u+sgma.*randn(size(u))+eps);

主函数改为
snr_in_db=0：1：10;
err_rate=PrimePe(snr_in_dB);
plot(snr_in_db,err_rate);

运算结果
>> snr_in_dB=0:1:10;
err_rate=PrimePe(snr_in_dB);
plot(snr_in_db,err_rate);
??? Error using ==> times
Matrix dimensions must agree.
Error in ==> PrimePe at 13
u=sign(u+sgma.*randn(size(u))+eps);

求解释

Answer 1

snr=10^(snr_in_dB/10);

这一句有问题，如果snr_in_dB/10＝A ＝ (Aij)是一个矩阵，我猜你想做下面两件事之一（感觉你是想做第一件）

1）返回一个矩阵B＝(Bij)，其中Bij = 10^Aij，这时你必须用运算符 .^
其中.表示按元素运算，即写成snr = 10 .^ (snr_in_dB/10)
2）如果A是一个方阵，那么10^A是有定义的，在数学上，
10^A = e^{ln10 * A}，而e^{B} 定义为e^B = 1+B+B^2/2!+...+B^n/n!+...
所以如果你没用.^这个运算，matlab就理解成为上面这种情况，但是要求A是一个方阵，因而报错。

追问

那么下面的相关的乘除都要改成点乘和点除？ 但是u=sign(u+sgma*randn(size(u))+eps);始终报错，求解

追答

如果是按元素运算，是的。

从你的代码看，sgma也是一个矩阵，所以也要用点乘。
不过你最好把错误信息贴出来

追问

以补充问题了，求指点

追答

错误信息中说得很清楚了，
??? Error using ==> times
Matrix dimensions must agree.
Error in ==> PrimePe at 13
u=sign(u+sgma.*randn(size(u))+eps);

就是说你的乘法（times）中，两个矩阵的维数不一样。
在按元素进行乘法时，sgma与u的维数必须一样。而你的u的维数是(1,N*frame_length)，sgma的维数与snr_in_dB相同，所以维数不一致。

追问

那么该怎么改

追答

我之前没理解你想干什么，现在知道了。如果你是想输入一列snr_in_dB，对每个snr_in_dB输出一个err_rate，一个比较简单的方法是写一个循环：
u = u*signal - 1;%% 从这里开始加入一个循环

err_rate = zeros(size(snr_in_dB));
for i = 1:length(snr_in_dB)
u1=sign(u+sgma(i)*randn(size(u))+eps); %% u1是一个临时变量，这时只用到sgma的第i分量
u1=(u1+1)/2;
% sgma(i)对应的误码率
err_rate(i)=(length(find(signal~=u1)))/length(signal);
end