| 查看: 1631 | 回复: 10 | ||
[求助]
matlab程序问题求解答
|
|
function[c,err,yc]=bisect(f,a,b,delta) % Input - f is the function input as a string 'f' % - a and b are the left and right endpoints % - delta is the tolerance % Output - c is the zero % - yc = f(c) % - err is the error estimate for c ya = feval(f,a); yb = feval(f,b); if ya * yb > 0, end max1 = 1+ round((log(b-a)-log(delta))/log(2)); for k =1:max1 c =(a+b)/2; yc = feval(f,c); if yc==0; a=c; b=c; elseif yb * yc > 0 b=c; yb=yc; else a=c; ya=yc; end if b - a < delta, break, end end c=(a+b)/2; err=abs(b-a); yc=feval(f,c); % 求函数在c处的函数值 for T=200:10:300 K = exp(-9.3932+3204.71/T+0.835931*log(T)... +2.35267*10^-3*T-1.8736*10^-6*T^2+5.1606*10^-10*T^3) f=inline('(4*K-1)*x^2-8*K*x+4*K'); bisect(f,0,1,10^-5) end 主程序调用函数出错 |
回复此楼» 猜你喜欢
朋友圈看到的
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有110人回复
-
明天放榜？
已经有5人回复
-
时间戳变了，能看出什么问题？
已经有11人回复
-
重要消息，中午系统在维护
已经有10人回复
-
两姐妹，一个辛辛苦苦30岁工作当大学教师，另一个23岁就工作当开车教练
已经有7人回复
-
投稿咨询
已经有6人回复
-
时间戳又变了8-15
已经有26人回复
-
今天维护系统维护 祝所有人 高中
已经有9人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教MATLAB中的LMI求解问题，急！！！
已经有7人回复
- 求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复
- matlab出错。。。求大神解答。。。
已经有3人回复
- matlab稀疏矩阵问题求教
已经有3人回复
- 在matlab里运行SVM的问题，期待着您的解答......
已经有5人回复
- 关于matlab线性规划 背包问题的 解法 求救呀~~方程知道 但不知道如何解答····
已经有4人回复
- 求助如何用matlab取出直线上的点
已经有10人回复
- 请帮忙用matlab求矩阵特征值
已经有8人回复
- matlab用多元线性最小二乘法求解多元线性方程的问题
已经有11人回复
- 求助啊！求一程序，用matlab程序做，用蒙特卡洛方法模拟
已经有10人回复
- matlab最小二乘问题求助 怎么定义多变量
已经有3人回复
- matlab出错 求助高手解答
已经有5人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 求助matlab软件压缩包
已经有4人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
- 如何利用MATLAB求最小值
已经有4人回复
- 【求助】求HMM的Matlab程序
已经有4人回复
2楼2013-11-24 17:21:08
3楼2013-11-24 17:46:06
4楼2013-11-24 22:09:26
5楼2013-11-24 22:21:01
6楼2013-11-24 22:22:50
Bruce_lee
铜虫 (著名写手)
- 应助: 5 (幼儿园)
- 金币: 2632.7
- 散金: 424
- 红花: 4
- 帖子: 1277
- 在线: 317小时
- 虫号: 948103
- 注册: 2010-01-25
- 性别: GG
- 专业: 聚合物共混与复合材料
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
yuyi2009: 金币+50, ★★★★★最佳答案, 再次感谢 2013-11-24 23:47:56
感谢参与，应助指数 +1
yuyi2009: 金币+50, ★★★★★最佳答案, 再次感谢 2013-11-24 23:47:56
|
function problem_solve %-------------------------------------------- disp('----------------分割线-----------------'); %-------------------------------------------- function[c,err,yc]=bisect(f,a,b,delta) % Input - f is the function input as a string 'f' % - a and b are the left and right endpoints % - delta is the tolerance % Output - c is the zero % - yc = f(c) % - err is the error estimate for c ya = feval(f,a); yb = feval(f,b); if ya * yb > 0, end max1 = 1+ round((log(b-a)-log(delta))/log(2)); for k =1:max1 c =(a+b)/2; yc = feval(f,c); if yc==0; a=c; b=c; elseif yb * yc > 0 b=c; yb=yc; else a=c; ya=yc; end if b - a < delta, break, end end c=(a+b)/2; err=abs(b-a); yc=feval(f,c); % 求函数在c处的函数值 end %--------------------------------------------- function f=subfun(x,K) f=(4*K-1).*x.^2-8*K.*x+4*K; end %--------------------------------------------- for T=200:10:300 K = exp(-9.3932+3204.71/T+0.835931*log(T)... +2.35267*10^-3*T-1.8736*10^-6*T^2+5.1606*10^-10*T^3) [c,err,yc]=bisect(@(x)subfun(x,K),0,1,10^-5) disp('----------------分割线-----------------'); end end 楼主，希望能帮助到你。 |
7楼2013-11-24 23:37:40
Bruce_lee
铜虫 (著名写手)
- 应助: 5 (幼儿园)
- 金币: 2632.7
- 散金: 424
- 红花: 4
- 帖子: 1277
- 在线: 317小时
- 虫号: 948103
- 注册: 2010-01-25
- 性别: GG
- 专业: 聚合物共混与复合材料
8楼2013-11-24 23:41:39
9楼2013-11-24 23:47:44
Bruce_lee
铜虫 (著名写手)
- 应助: 5 (幼儿园)
- 金币: 2632.7
- 散金: 424
- 红花: 4
- 帖子: 1277
- 在线: 317小时
- 虫号: 948103
- 注册: 2010-01-25
- 性别: GG
- 专业: 聚合物共混与复合材料
10楼2013-11-25 15:34:54
