引用:
作者: me4047
其实这些命令我都知道！
但是我现在的问题是，我的变量是在for的循环里面，画出来到图总是一个点！
比如说，for h=1:10,
但是画出来 的图就是h=10的时候根本没有别的！
|
你这个小程序写的问题非常多，一个一个来：
function f=ps(M,h,k,sn) % h, k, sn 为函数的形参，但你却在下面再次给它们赋值，那么你把它
% 们写在形参list里又有什么意义呢？
for h = 0:13
for sn = 1:10
y = 0;
for k = 0:M-h
pe = erfc( sqrt(sn) ) / 2;
d = nchoosek(M,k).*( (1-pe).^(M-k) ).*pe.^k;
y = y + d;
end
pm = 1 - y;
z=0;
for k = 0:h
c = nchoosek( M, k );
z = y + c;
end
pf = z.*(1/2).^M;
ps = pm + pf % 你这个函数名是ps，定义了一个返回值f，这里却给函数名赋值？如果这里
% 加上分号，那你就连ps都看不到了。
end
end
=========================================================
下面是我简单修改了的，为了方便你看懂，没有改你的程序结构，只出于解决你的问题略作修改
function [h, f] = ps( M )
h = 0:13; % h 已知，这里直接定义为数组，下面采用h(i)的方式顺序访问h中的元素。
j = 1; % 返回值f是个数组，长度我不知道，定义个 j 进行loop自增。
for i = 1:length(h)
for sn = 1:10
y = 0;
for k = 0:M-h(i)
pe = erfc( sqrt(sn) ) / 2;
d = nchoosek(M,k).*( (1-pe).^(M-k) ).*pe.^k;
y = y + d;
end
pm = 1 - y;
z=0;
for k = 0:h(i)
c = nchoosek( M, k );
z = y + c;
end
pf = z.*(1/2).^M;
f(j) = pm + pf; % 给返回值 f 赋值，f 必然是个数组
j = j + 1;
end
end
%% 完毕，在命令行里输入[ h, f ] = ps( M的值 ) 即可运行，返回 h 和 f。但是，这里依然是没法绘图的，因为由程序可以看出 1个 h 会有 10 个f 输出，我不了解你的程序意图，不敢说我这里逻辑是对的，但程序本身已经可以运行并返回你想要的两组值，你仔细想想修改一下，plot出来就可以。