初学,MATLAB方程画图的问题,请多指教
来源:百度知道 编辑:UC知道 时间:2024/06/07 01:17:26
通过数学推倒得出一个公式,想用MATLAB画出公式对应的曲线
我的程序
syms n pd k N ii M Nv r pv j
M=20
N=16
k=2
Nv=14
p0=symsum((pd/k)^n/sym('n!'),n,0,((N-ii)/k)-1)+pd^((N-ii)/k)/(sym('((N-ii)/k)!')*k^((N-ii)/k))*symsum(pd/(N-ii)^(n-((N-ii)/k)),n,((N-ii)/k),M)
pbb=(pv^ii/sym('ii!')/symsum(pv^r/sym('r!'),r,0,Nv))*(pd^M/((N-ii)^(M-((N-ii)/k))*sym('((N-ii)/k)!')*k^((N-ii)/k)))*p0^(-1)
pb=symsum(pbb,ii,0,Nv)
ls=symsum((pd/k)^j/sum('(j-1)!'),j,1,((N-ii)/k)-1)+symsum(j*(pd^j)*(p0^(-1))/((N-ii)^(j-((N-ii)/k))*sum('((N-ii)/k)!')*k^((N-ii)/k)),j,((N-ii)/k),M)
td=ls/((pd/0.04)*(1-pd))-1/0.04
ezplot(td)
运行出现错误,哪里的问题呢,请高手指教,谢谢
式中的pv,pd,ii为变量,需要画出其中一个作为变量,另两个为固定值时的图形
例如pv=20,pd=10,ii=0:14时的图形
我的程序
syms n pd k N ii M Nv r pv j
M=20
N=16
k=2
Nv=14
p0=symsum((pd/k)^n/sym('n!'),n,0,((N-ii)/k)-1)+pd^((N-ii)/k)/(sym('((N-ii)/k)!')*k^((N-ii)/k))*symsum(pd/(N-ii)^(n-((N-ii)/k)),n,((N-ii)/k),M)
pbb=(pv^ii/sym('ii!')/symsum(pv^r/sym('r!'),r,0,Nv))*(pd^M/((N-ii)^(M-((N-ii)/k))*sym('((N-ii)/k)!')*k^((N-ii)/k)))*p0^(-1)
pb=symsum(pbb,ii,0,Nv)
ls=symsum((pd/k)^j/sum('(j-1)!'),j,1,((N-ii)/k)-1)+symsum(j*(pd^j)*(p0^(-1))/((N-ii)^(j-((N-ii)/k))*sum('((N-ii)/k)!')*k^((N-ii)/k)),j,((N-ii)/k),M)
td=ls/((pd/0.04)*(1-pd))-1/0.04
ezplot(td)
运行出现错误,哪里的问题呢,请高手指教,谢谢
式中的pv,pd,ii为变量,需要画出其中一个作为变量,另两个为固定值时的图形
例如pv=20,pd=10,ii=0:14时的图形
ls=symsum((pd/k)^j/sum('(j-1)!'),j,1,((N-ii)/k)-1)+symsum(j*(pd^j)*(p0^(-1))/((N-ii)^(j-((N-ii)/k))*sum('((N-ii)/k)!')*k^((N-ii)/k)),j,((N-ii)/k),M)
这句里面的pd是什么啊,好像没给出其表达式。
============================
补充:
你这里的解析表达式,这个问题太复杂了。
你最好把错误也补上