全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 MATLAB等数学软件专版
2304 14
2016-05-14
悬赏 1800 个论坛币 未解决
下面是我附的代码,给定t值,求w值,一个t对应两个w解,然后将w结果带入方程检验,但是发现有误差。wj2(2315,:)有一个值不为零,想问一下,有谁会用二分法求该方程的解,或者提高该方程的运算精度!

%给定m,n,rou,alpha,opw为p(rou,w)取得最小值,w值,求w两解
clc;clear
m=500;n=5;rou=1;alpha=0.0023;
s=rou*m*(n-1);
ucl=chi2inv(1-alpha/2,n-1);
lcl=chi2inv(alpha/2,n-1);
opw=rou*m*(n-1)^2*((log(ucl/lcl))/(ucl-lcl));w=zeros(2315,2);
y=1/(1-chi2cdf(opw*ucl/s,n-1)+chi2cdf(opw*lcl/s,n-1))
for t=1.1:541.1
w(t,:)=fsolve(@(w) 1-chi2cdf(w*ucl/s,n-1)+chi2cdf(w*lcl/s,n-1)-1/t,[1 opw],optimset('Display','off'));
%检验结果w对与否
wj1=zeros(2315,1);
wj2=zeros(2315,1);
wj1(2315,:)=1-chi2cdf(w(t,1)*ucl/s,n-1)+chi2cdf(w(t,1)*lcl/s,n-1)-1/t
wj2(2315,:)=1-chi2cdf(w(t,2)*ucl/s,n-1)+chi2cdf(w(t,2)*lcl/s,n-1)-1/t
end


6W8OV[P)A~6HXJUIKO(2JML.png

原图尺寸 7.51 KB

6W8OV[P)A~6HXJUIKO(2JML.png

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2016-5-14 22:08:45
有偿,会的可以谈!
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2016-5-15 03:09:18
我写了一个函数,你把 w(t,:)=... 的那一行换成 w(t,:)=f1(t,m,n,rho,alpha,prec); 就可以了,而prec就是你需要的精确度,比如1e-6或者1e-12

附件列表

f1.zip

大小:376 Bytes

只需: 1 个论坛币  马上下载

本附件包括:

  • f1.m

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2016-5-15 03:11:30
请看我的答复:

https://bbs.pinggu.org/forum.php?mod=redirect&goto=findpost&ptid=4611737&pid=36345645&from^^uid=5975757
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2016-5-15 11:13:03
貌似你要求解的方程只有一个根, 而不是像你说的一个t对应两个w, 下面是t=1.1时, 对应函数的图像, 可以发现其是单调递减的, 只有一个零点. untitled.jpg
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2016-5-15 15:53:20
enxizheng 发表于 2016-5-15 11:13
貌似你要求解的方程只有一个根, 而不是像你说的一个t对应两个w, 下面是t=1.1时, 对应函数的图像, 可以发现其 ...
函数我证明过了,两个解类似于抛物线形式,还有个解在t比较小时10000以上
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群