24小时热门版块排行榜    

查看: 2586  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 281求调剂(0805) +3 烟汐忆海 2026-03-16 8/400 2026-03-18 00:47 by 烟汐忆海
[考研] 302求调剂 +6 呼呼呼。。。。 2026-03-17 6/300 2026-03-17 23:56 by 星空星月
[考研] 344求调剂 +4 knight344 2026-03-16 4/200 2026-03-17 17:27 by ruiyingmiao
[考研] 考研化学学硕调剂,一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 工科材料085601 279求调剂 +3 困于星晨 2026-03-17 3/150 2026-03-17 14:08 by ms629
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 0703化学调剂,求各位老师收留 +8 秋有木北 2026-03-14 8/400 2026-03-16 15:21 by 哦哦123
[考研] 085600材料与化工 求调剂 +13 enenenhui 2026-03-13 14/700 2026-03-16 15:19 by 了了了了。。
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 265求调剂 +4 威化饼07 2026-03-12 4/200 2026-03-14 17:23 by userper
[考研] 材料与化工(0856)304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 求材料调剂 +5 隔壁陈先生 2026-03-12 5/250 2026-03-13 22:03 by 星空星月
[考研] 材料与化工085600调剂求老师收留 +9 jiaanl 2026-03-11 9/450 2026-03-13 20:22 by JourneyLucky
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 工科调剂 +4 Jiang191123! 2026-03-11 4/200 2026-03-13 15:15 by Miko19
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
信息提示
请填处理意见