查看: 277|回復: 0|關注: 0
打印 上一主題 下一主題

[未答復] 構造一個函數文件,是一個似然函數,但里面的積分過程不會處理。

[復制鏈接]

新手

5 麥片

財富積分


050


2

主題

4

帖子

0

最佳答案
跳轉到指定樓層
1#
發表于 2020-1-29 15:02:31 | 只看該作者 回帖獎勵 |倒序瀏覽 |閱讀模式
%編寫的似然函數文件
function p=p(x)
k=[3,3,3,3,3,3,3,3,3,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,5,5,5,5,5,5,5,6,6,6,6,7,7,8];%累乘的元素
b=log((x+1)./x)./x;%默認量
pc=1;
syms w
for i=1:length(k)    %關鍵項
f=w.^k(i).*exp(-(x+1).*w);   %關鍵項
pc=pc.*((x+1).^(-k(i))+x.*exp(x.*b)./fact(k(i)).*int(f,w,0,b));    %關鍵項
end

s=min(k);
pn=0;
for j=1:s-1   %關鍵項
     g=w.^j.*exp(-(x+1).*w);%關鍵項
     pn=pn+((x+1).^(-j)+x.*exp(x.*b)./fact(j).*int(g,w,0,b));%關鍵項
end

p=pc.*exp(-x)./(1-pn).^length(k);
end


構造這個函數文件后,想測試一下,于是"x=linspace(0,5,100);plot(x,p(x))"想畫個圖都無法實現,就連基本的p(0.7)vsp(1.7)比大小也無法返回邏輯值。其中上面的w是構造函數的過渡項,后面積分就沒了,應該是只剩下x的函數。我覺得的關鍵應該是我插入的兩個積分,我們來是想用[email protected](w)...構造一個函數,但是x是變量,積分里有變量也無法實現。

回復主題 已獲打賞: 0 積分

舉報

您需要登錄后才可以回帖 登錄 | 注冊

本版積分規則

關閉

站長推薦上一條 /4 下一條

快速回復 返回頂部 返回列表
哪一款德州扑克还能玩 选4开奖结果上海彩控 辉煌app7606 车联网功能 网络上怎么赚钱 经典单机四人麻将 在家创业小本办厂 nba集锦 宝博棋牌游戏下载地址 20选8快乐十分开奖结果山西 炒股平台下载 遇乐棋牌大厅安装 六合历史开奖记录 麻将来了血流换四张玩法 北京赛车pk10方案 意甲最新积分榜 能猫四川麻将血战到