应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 怎么将此for循环改为矩阵形式？
61/1返回列表
查看: 2581  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

化名为哀

铜虫 (初入文坛)

[求助] 怎么将此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矩阵维度不匹配，求高人的指点如何修改？
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-11-25 16:56:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖

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));
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

化名为哀
一勤天下无难事，处世无难忍最高。
5楼2012-11-27 12:02:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

化名为哀

铜虫 (初入文坛)
送鲜花一朵
引用回帖:
5楼: Originally posted by mei3014 at 2012-11-27 12:02:04
% 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=;
y(m)=sin(2*pi*fc(code(floor(t(m)/(Td))+1)+1).* ...

已解决 谢谢
回复此楼
6楼2012-11-27 14:10:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 化名为哀 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 机械专硕325，寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 一志愿华中师范071000，325求调剂 +6 RuitingC 2026-03-12 6/300 2026-03-16 14:50 by 可淡不可忘
[考研] 309求调剂 +5 花与叶@ 2026-03-10 5/250 2026-03-16 14:13 by 哦哦123
[考研] 304求调剂 +6 小熊joy 2026-03-14 6/300 2026-03-16 12:59 by Iveryant
[考研] 0703 物理化学调剂 +3 我可以上岸的对� 2026-03-13 5/250 2026-03-16 10:50 by 我可以上岸的对�
[考研] 290求调剂 +5 孔志浩 2026-03-12 10/500 2026-03-16 09:01 by 余晖&
[考研] 311求调剂 +6 冬十三 2026-03-15 6/300 2026-03-16 08:00 by wang_dand
[考研] 255求调剂 +3 李嘉慧， 2026-03-12 4/200 2026-03-14 16:58 by 有只狸奴
[考研] 267一志愿南京工业大学0817化工求调剂 +5 SUICHILD 2026-03-12 5/250 2026-03-14 14:53 by jean5056
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 0856材料与化工309分求调剂 +6 ZyZy…… 2026-03-10 6/300 2026-03-14 00:38 by JourneyLucky
[考研] 求调剂，一志愿江南大学环境工程085701 +3 Djdjj12 2026-03-10 4/200 2026-03-14 00:31 by JourneyLucky
[考研] 279求调剂 +3 Dizzy123@ 2026-03-10 3/150 2026-03-13 23:02 by JourneyLucky
[考研] 285化工学硕求调剂（081700） +6 柴郡猫_ 2026-03-12 6/300 2026-03-13 20:46 by hmn_wj
[考研] 085600材料与化工 309分请求调剂 +7 dtdxzxx 2026-03-12 8/400 2026-03-13 14:43 by jxchenghu
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生，2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
[考研] 化工0817调剂 +8 灿若星晨 2026-03-10 8/400 2026-03-10 22:44 by 星空星月
信息提示
关闭
请填处理意见
关闭