![]() |
[求助]画不出图来!
请现看看我写的程序
function f=ps(M,h,k,sn) 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 end end 我费劲心血写出来的,不要笑话我啊:sweat: 我现在虽然能写出来ps的答案,但是画不出图啊! 请问怎么样才可以画出一个ps和h的图像呢 ? 感谢先 |
楼主,你那是作图程序吗?
|
是个公式
不过也要画这个公式的图象啊! |
2 个附件
... ...
看来你对maltab绘图还没什么认识,给你个 基础入门的绘图教程吧。中文的。 压缩包解压后里面有个index.html,打开这个就是目录,点击目录可以浏览。 |
[QUOTE=remnant;13550]... ...
看来你对maltab绘图还没什么认识,给你个 基础入门的绘图教程吧。中文的。 压缩包解压后里面有个index.html,打开这个就是目录,点击目录可以浏览。[/QUOTE] 还没看,先谢谢! |
其实这些命令我都知道!
但是我现在的问题是,我的变量是在for的循环里面,画出来到图总是一个点! 比如说,for h=1:10, 但是画出来 的图就是h=10的时候根本没有别的! |
[QUOTE=me4047;13580]其实这些命令我都知道!
但是我现在的问题是,我的变量是在for的循环里面,画出来到图总是一个点! 比如说,for h=1:10, 但是画出来 的图就是h=10的时候根本没有别的![/QUOTE] 你这个小程序写的问题非常多,一个一个来: 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出来就可以。 |
[QUOTE=remnant;13582]你这个小程序写的问题非常多,一个一个来:
function f=ps(M,h,k,sn) % h, k, sn 为函数的形参,但你却在下面再次给它们赋值,那么你把它 % 们写在形参list里又有什么意义呢? for ...[/QUOTE] :biggrin: 先感谢, 再看帖! :biggrin: |
老大
看了半天我也没看明白哪里出问题了,算来算去,不管M是多少,答案都是一样的啊! 冥思苦想中................ |
这个要看你的题目是什么了,也就是你的程序要完成什么function.
从你的程序里看到,每次取一个h,然后sn 从1循环到10。每个sn计算出1个f。 |
所有时间均为北京时间。现在的时间是 06:49。 |
Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.