| 查看: 1447 | 回复: 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 主程序调用函数出错 |
回复此楼» 猜你喜欢
【考研调剂】化学专业 281分，一志愿四川大学，诚心求调剂
已经有11人回复
-
一志愿深大，0703化学，总分302，求调剂
已经有4人回复
-
材料学硕333求调剂
已经有3人回复
-
材料求调剂
已经有3人回复
-
333求调剂
已经有5人回复
-
求助
已经有6人回复
-
278求调剂
已经有9人回复
-
297求调剂
已经有12人回复
-
0703化学297求调剂
已经有3人回复
-
材料学硕301分求调剂
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教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人回复
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
2楼2013-11-24 17:21:08
3楼2013-11-24 17:46:06
4楼2013-11-24 22:09:26
