24小时热门版块排行榜    

查看: 2587  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

化名为哀

铜虫 (初入文坛)

[求助] 怎么将此for循环改为矩阵形式?

M=2;       % M进制调制

BW=25e+3;
delta_f=1e+3;
BW1=(BW-2*delta_f)/2
Td=1/BW1   % 码元周期
fd=1/Td
fs=15.5e6;   % 采样频率
f1=1e+6;
f2=f1+delta_f;
BW=2*BW1+2*delta_f
fc=[f1,f2];    % 载频
counter=Td*fs;   
T=20e-3;    % 模拟信号长度
N=T*fs     % 总的样本点数
Ns=floor(N/counter)     % 码元数
code0=rand(1,Ns);
code=fix(code0*M); % 产生随机码元
step = 1/fs;
t= 0: step: Td*Ns-step;
length(t)
for m=1:length(t)
   a=floor(t(m)/(Td));
   y(m)=sin(2*pi*fc(code(a+1)+1)*t(m));
   phase_shift=exp(-i*2*pi*fc(code(a+1)+1));
end;

% m=[1:length(t)];
% y(m)=sin(2*pi*fc(code(floor(t(m)/(Td))+1)+1)*t(m));
% phase_shift=exp(-i*2*pi*fc(code(floor(t(m)/(Td))+1)+1));

程序如上,可正常运行,但考虑到for循环耗时太久,所以讲for循环改为下面注释的部分,改之后提示inner矩阵维度不匹配,求高人的指点如何修改?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mei3014

金虫 (正式写手)

木有头屑

【答案】应助回帖


xiegangmai: 金币+1, 辛苦了! 2012-11-27 22:51:07
% for m=1:length(t)
%    a=floor(t(m)/(Td));
%    y(m)=sin(2*pi*fc(code(a+1)+1)*t(m));
%    phase_shift=exp(-i*2*pi*fc(code(a+1)+1));
% end;

m=[1:length(t)];
y(m)=sin(2*pi*fc(code(floor(t(m)/(Td))+1)+1).*t(m));
phase_shift=exp(-i*2*pi*fc(code(floor(t(m)/(Td))+1)+1));

» 本帖已获得的红花(最新10朵)

一勤天下无难事,处世无难忍最高。
5楼2012-11-27 12:02:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

zkw666666

木虫 (正式写手)

益虫。。。。

★ ★
xiegangmai: 金币+2, 谢谢参与 2012-11-26 09:29:25
a=floor(t/(Td));
    y=sin(2*pi*fc(code(a+1)+1).*t);
    phase_shift=exp(-i*2*pi*fc(code(a+1)+1));
2楼2012-11-25 22:07:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zkw666666

木虫 (正式写手)

益虫。。。。

【答案】应助回帖

感谢参与,应助指数 +1
Elapsed time is 0.308414 seconds.
3楼2012-11-25 22:09:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mei3014

金虫 (正式写手)

木有头屑

【答案】应助回帖


感谢参与,应助指数 +1
xiegangmai: 金币+1, 谢谢参与 2012-11-27 22:50:58
% y(m)=sin(2*pi*fc(code(floor(t(m)/(Td))+1)+1).*t(m));
就是矩阵之间运算的时候缺个点。我已经运行过了,将*改为.*就ok了。
一勤天下无难事,处世无难忍最高。
4楼2012-11-27 12:01:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[基金申请] 被我言中:新模板不强调格式了,假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 能源材料化学课题组招收硕士研究生8-10名 +3 脱颖而出 2026-03-16 6/300 2026-03-17 21:19 by z1z2z3879
[考研] 302求调剂 +9 负心者当诛 2026-03-11 9/450 2026-03-17 17:13 by ruiyingmiao
[考研] 26考研求调剂 +6 丶宏Sir 2026-03-13 6/300 2026-03-17 16:13 by 醉在风里
[考研] 材料与化工304求B区调剂 +7 邱gl 2026-03-11 8/400 2026-03-17 09:36 by 努力学习赚彩礼
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 材料与化工 323 英一+数二+物化,一志愿:哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 308求调剂 +5 是Lupa啊 2026-03-11 5/250 2026-03-13 22:13 by JourneyLucky
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[考研] 一志愿西南交大,材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
信息提示
请填处理意见