热门
最新
红包
立Flag
投票
同城
我的
发布
yunweigo
3 年前
truezwolfer
保持热爱共赴山海
下一条：
有朋友解释一下这几行MATLAB代码不？
syms x Q Qm jieguo1 jieguo2;f=1/(4.74*sqrt(2*pi)).*exp(-1/2*((x-11.73)/4.74).^2);eq=int(f,x,-inf,Q)-0.77;Q=solve(eq,Q);VQ=vpa(Q);Q=fix(VQ);disp(Q);eq1=int(f,x,-inf,Qm)-0.99;Qm=solve(eq1,Qm);VQm=vpa(Qm);Qm=fix(VQm);disp(Qm);m=fix(Q);n=fix(Qm);h=1:1:n+1;f1=1:1:n+1;h(1)=int(f,-inf,0);f1(1)=h(1);for i=2:n+1 h(i)=int(f,-inf,i); f1(i)=h(i)-h(i-1);endjieguo1=0;jieguo2=0;for i=1:n+1 jieguo1=jieguo1+15*(m-i+1)*f1(i);endfor i=1:m+2 jieguo2=jieguo2+15*(m+1-i+1)*f1(i);endfor i=m+1:n+1 jieguo1=jieguo1+50*(i-1-m)*f1(i);endfor i=m+2:n+1 jieguo2=jieguo2+50*(i-1-m-1)*f1(i);endjieguo1=vpa(jieguo1);disp(fix(jieguo1))jieguo2=vpa(jieguo2);disp(fix(jieguo2))
立即登录
syms x Q Qm jieguo1 jieguo2;f=1/(4.74*sqrt(2*pi)).*exp(-1/2*((x-11.73)/4.74).^2);eq=int(f,x,-inf,Q)-0.77;Q=solve(eq,Q);VQ=vpa(Q);Q=fix(VQ);disp(Q);eq1=int(f,x,-inf,Qm)-0.99;Qm=solve(eq1,Qm);VQm=vpa(Qm);Qm=fix(VQm);disp(Qm);m=fix(Q);n=fix(Qm);h=1:1:n+1;f1=1:1:n+1;h(1)=int(f,-inf,0);f1(1)=h(1);for i=2:n+1 h(i)=int(f,-inf,i); f1(i)=h(i)-h(i-1);endjieguo1=0;jieguo2=0;for i=1:n+1 jieguo1=jieguo1+15*(m-i+1)*f1(i);endfor i=1:m+2 jieguo2=jieguo2+15*(m+1-i+1)*f1(i);endfor i=m+1:n+1 jieguo1=jieguo1+50*(i-1-m)*f1(i);endfor i=m+2:n+1 jieguo2=jieguo2+50*(i-1-m-1)*f1(i);endjieguo1=vpa(jieguo1);disp(fix(jieguo1))jieguo2=vpa(jieguo2);disp(fix(jieguo2))