大家好! 小弟我是matlab的新手, 最近在撰写程式过程中遇到一点问题, 曾试着先爬文或从相关教学书籍甚至google网站寻找解答, 但一直无法解决, 因此来PO文,希望有好心大大可以帮忙解惑! 其实我主要要求解的程式如下列网址 http://imgur.com/a/mval1 而我的程式码(M档案)如下: %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function [result,minimalvalue,Tvalue] = project20160901RFP(q) format long g T = 1:1:10 ; Cp = 100 ; Cf = 500 ; result = zeros(1,length(T)) ; for i = 1:length(T) syms y x result(1,i) = int(( ( Cp + int(fun3(x,q)*q*fun2(x),x,0,y) ) / y ) * fun1(y,q) * fun5(y),y,0,T(i)) + ( ( Cp + int(fun3(x,q)*q*fun2(x),x,0,T(i)) ) / T(i) ) * fun1(T(i),q) * fun6(T(i)) + int(( ( Cf + int(fun3(x,q)*q*fun2(x),x,0,y) ) / y ) * fun6(y) * fun1(y,q) * (1-q) * fun2(y),y,0,T(i)) ; end [minimalvalue, I] = min(result(:)) ; [order] = ind2sub(size(result), I) ; Tvalue = T(order) ; function Fpbary = fun1(y,q) syms t Fpbary = exp( (-1) * int((1-q)*fun2(t),t,0,y) ) ; function rofy = fun2(y) beta = 2 ; theta = 5 ; rofy = (beta/theta)*((y/theta)^(beta-1)) ; function hofy = fun3(y,q) syms t mean = 100 ; sigma = 20 ; cinf = 1000 ; constant = 0.2 ; delta = ( ( icdf('norm',(cdf('norm',0,mean,sigma) + q),mean,sigma) ) / cinf ) ; hofy = ( ((1/q)*int(t*fun4(t),t,0,delta*cinf)) + constant*y ) ; function aoft = fun4(t) mean = 100 ; sigma = 20 ; aoft = (1/(sigma*((2*pi)^(1/2))))*exp((-1)*((t-mean)^2)/(2*(sigma^2))) ; function gofy = fun5(y) k = 2 ; lamda = 5 ; gofy = ( k / (lamda^k) ) * ( y^(k-1) ) * exp(-1 * ( (y/lamda)^k ) ) ; function Gbarofy = fun6(y) k = 2 ; lamda = 5 ; Gbarofy = exp(-1 * ( (y/lamda)^k ) ) ; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 跑完程式後,会有计算出来的相关数值,也会有下列的讯息: Warning: Explicit integral could not be found. 我知道这讯息的意思是积分无 close-form solution,但我的问题如下: 1. 这样的话 matlab所算出来的值是否正确? (我知道可以徒手再去计算检验,但计算式真的不容易徒手去算) 2. 如果所算出来的值正确,那怎会出现上述Waring的讯息? 3. 承2.,如果算出来的值不正确,那我该如何修改此程式码? 我的程式码有点冗长,或许不易阅读,真的很不好意思! 希望有好心的大大愿意帮忙解惑,由衷感谢! --

※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 123.110.47.54 ※ 文章网址: https://webptt.com/cn.aspx?n=bbs/MATLAB/M.1472828593.A.CE0.html ※ 编辑: nago (123.110.47.54), 09/02/2016 23:04:55