24小时热门版块排行榜

返回列表

waitingywy

新虫 (小有名气)

应助: 1 (幼儿园)
金币: 1223
散金: 168
红花: 3
帖子: 247
在线: 64.7小时
虫号: 2220691
注册: 2013-01-04
性别: GG
专业: 零件成形制造

[求助] MATLAB出图小程序找错误，求帮忙已有1人参与

这段代码无法出图，不是很懂，求大神指教

bats_codes_out=cell(5,M);
sucess_rates=[];
m=10;
sucess_rate_sum=0;
fail_rate_sum=0;
for i=1:1:m
n=10;                                        % round((1+epsilong)*K/M);
mu_n=randsample([1:length(mu)],n,true,mu);
D=max(mu_n);
d_min=ceil(D/K);
bats_codes=[];
for j=1:n
      bats_codes_out=bats_encode(d_min,j,q,codes,mu_n,K,g,L,M);
      bats_codes_bec=becpacket(M,e,bats_codes_out,g);
      bats_codes_in=bats_codes_inner(q,M,bats_codes_bec,g);
      bats_codes_bec_bec=becpacket(M,e,bats_codes_in,g);
      bats_codes=[bats_codes,bats_codes_bec_bec];
end
[fail_rate,de_codes,x]=bats_decodes(bats_codes,M,n,K,L,g);
sucess_rate(i)=1-fail_rate
X(i)=x;
sucess_rate_sum=sucess_rate_sum+sucess_rate(i);
fail_rate_sum=fail_rate_sum+fail_rate;
end
sucess_aver=sucess_rate_sum/m;
fail_rate_aver=fail_rate_sum/m;

suc_one=0;
for i=1:m
if sucess_rate(i)==1
      suc_one=suc_one+1;
end
end

% for i=1:K
% if ~(isempty(de_codes{2,i}))
% j=j+1;
% end
% end

% s1=1:length(mu);
% subplot(121);
% plot(s1,mu);
% s2=1:length(mu1);
% subplot(122);
% plot(s2,mu1);

回复此楼

» 猜你喜欢

反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有125人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复

十年

1楼 2020-05-09 18:04:17

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

waitingywy

新虫 (小有名气)

应助: 1 (幼儿园)
金币: 1223
散金: 168
红花: 3
帖子: 247
在线: 64.7小时
虫号: 2220691
注册: 2013-01-04
性别: GG
专业: 零件成形制造

clear all;
clc;
g=8;
q=2^g;
K=100;
L=2;
e=0.1;
M=16;
codes=randsrc(L,K,[0:q-1]);
codes=gf(codes,g);
mu=importdata('simDegreeM16K100p0.1.mat');
for i=1:size(mu,1)
if i>K
      mu(K)=mu(K)+mu(i);
      mu(i)=0;
end
end
sum(mu)
bats_codes_out=cell(5,M);
sucess_rates=[];
m=10;
sucess_rate_sum=0;
fail_rate_sum=0;
for i=1:1:m
n=10;                                        % round((1+epsilong)*K/M);
mu_n=randsample([1:length(mu)],n,true,mu);
D=max(mu_n);
d_min=ceil(D/K);
bats_codes=[];
for j=1:n
      bats_codes_out=bats_encode(d_min,j,q,codes,mu_n,K,g,L,M);
      bats_codes_bec=becpacket(M,e,bats_codes_out,g);
      bats_codes_in=bats_codes_inner(q,M,bats_codes_bec,g);
      bats_codes_bec_bec=becpacket(M,e,bats_codes_in,g);
      bats_codes=[bats_codes,bats_codes_bec_bec];
end
[fail_rate,de_codes,x]=bats_decodes(bats_codes,M,n,K,L,g);
sucess_rate(i)=1-fail_rate
X(i)=x;
sucess_rate_sum=sucess_rate_sum+sucess_rate(i);
fail_rate_sum=fail_rate_sum+fail_rate;
end
sucess_aver=sucess_rate_sum/m;
fail_rate_aver=fail_rate_sum/m;

suc_one=0;
for i=1:m
if sucess_rate(i)==1
      suc_one=suc_one+1;
end
end

% for i=1:K
% if ~(isempty(de_codes{2,i}))
% j=j+1;
% end
% end

% s1=1:length(mu);
% subplot(121);
% plot(s1,mu);
% s2=1:length(mu1);
% subplot(122);
% plot(s2,mu1);