| 查看: 1914 | 回复: 0 | |||
[求助]
切比雪夫函数逼近的问题
|
|
本人是初学者，一本教材上有这个函数的代码，但是在matlaB中运行会出问题，请高手指点一下 代码是 function f = Chebyshev(y,k,x0) %用切比雪夫多项式逼近已知函数 %已知函数：y %逼近已知函数所需项数：k %逼近点的x坐标：x0 %求得的切比雪夫逼近多项式或在x0处的逼近值：f syms t k; T(1:k+1) = t; T(1) = 1; T(2) = t; c(1:k+1) = 0.0; c(1)=int(subs(y,findsym(sym(y)),sym('t'))*T(1)/sqrt(1-t^2),t,-1,1)/pi; c(2)=2*int(subs(y,findsym(sym(y)),sym('t'))*T(2)/sqrt(1-t^2),t,-1,1)/pi; f = c(1)+c(2)*t; for i=3:k+1 T(i) = 2*t*T(i-1)-T(i-2); c(i) = 2*int(subs(y,findsym(sym(y)),sym('t'))*T(i)/sqrt(1-t^2),t,-1,1)/2; f = f + c(i)*T(i); f = vpa(f,6); if(i==k+1) if(nargin == 3) f = subs(f,'t',x0); else f = vpa(f,6); end end end 输入后显示的错误是： ??? Error using ==> sym.colon at 27 Cannot compute the number of steps from 1 to k + 1 by 1. Error in ==> Chebyshev at 9 T(1:1:k+1) = t; 求高手指点 |
回复此楼» 猜你喜欢
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
物理学I论文润色/翻译怎么收费?
已经有288人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
找到一些相关的精华帖子，希望有用哦~
【分享】Excel 2003 函数应用大全—教你搞定函数问题【已搜索无重复】
已经有205人回复
- 【求助】用graythresh函数遇到的问题
已经有8人回复
- 【求助】怎么使用切比雪夫逼近函数
已经有4人回复
- 【求助】Mittag-leffler函数问题
已经有3人回复
- 【求助】请教一个隐函数的求导问题
已经有10人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
