matlab中的符号运算与微积分_图文_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
评价文档：
matlab中的符号运算与微积分
上传于||文档简介
&&m​a​t​l​a​b​中​的​符​号​运​算​与​微​积​分
大小：442.00KB
登录百度文库，专享文档复制特权，财富值每天免费拿！
你可能喜欢第3章 MATLAB的符号运算_微分方程求解_符号代数方程_图文_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
第3章 MATLAB的符号运算_微分方程求解_符号代数方程
上传于||文档简介
&&较​详​细​介​绍​了​m​a​t​l​a​b​的​符​号​运​算​,​包​括​【​符​号​运​算​入​门​】​,​【​符​号​算​术​运​算​、​符​号​矩​阵​的​加​、​减​运​算​】​,​【​符​号​微​积​分​运​算​】​【​符​号​微​分​方​程​,​代​数​方​程​求​解​】​【​e​z​p​l​o​t​画​图​】​【​Z​变​换​,​傅​里​叶​变​换​,​拉​式​变​换​,​及​其​各​自​的​反​变​换​】​【​m​a​p​l​e​介​绍​】​【​M​a​t​l​a​b​符​号​函​数​图​形​计​算​器​F​U​N​T​O​O​L​】
阅读已结束，如果下载本文需要使用1下载券
想免费下载本文？
定制HR最喜欢的简历
下载文档到电脑，查找使用更方便
还剩20页未读，继续阅读
定制HR最喜欢的简历
你可能喜欢MATLAB切比雪夫逼近代码问题_百度知道
MATLAB切比雪夫逼近代码问题
function f=Chebyshev(y,k,x0)T(1:k+1) =T(1) = 1;T(2) =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)/c(2)=2*int(subs(y,findsym(sym(y)),sym('t'))*T(2)/sqrt(1-t^2),t,-1,1)/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)/
f=f+c(i)*T(i);
f=vpa(f,6);
if(i==k+1)
if(nargin==3)
f=subs(f,'t',x0);
f=vpa(f,6);
endend这是我的功能函数代码，我在MATLAB里面运行了以下程序： f=Chebyshev('1/(2-x)',6)结果就报错Warning: Explicit integral could not be found.
??? The following error occurred converting from sym to double:Error using ==& mupadmexError in MuPAD command: DOUBLE cannot convert the input expression into adouble array.If the input expression contains a symbolic variable, use the VPA functioninstead.Error in ==& Chebyshev at 8c(1)=int(subs(y,findsym(sym(y)),sym('t'))*T(1)/sqrt(1-t^2),t,-1,1)/求大神指教，这是书上的源码反复核对了几遍都没有找到和书上有什么不同，但是就是有问题。
提问者采纳
程序没有问题，可以运行，只是编写得有些繁琐，且结果并未化简，将m文件修改如下： function f=Chebyshev(y,k,x0)T(1)=sym(1);T(2)=t;c(1)=int(subs(y,findsym(sym(y)),t)*T(1)/sqrt(1-t^2),t,-1,1)/c(2)=2*int(subs(y,findsym(sym(y)),t)*T(2)/sqrt(1-t^2),t,-1,1)/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)),t)*T(i)/sqrt(1-t^2),t,-1,1)/
f=f+c(i)*T(i);
f=simple(vpa(f,6));
if(i==k+1)
if(nargin==3)
f=simple(vpa(subs(f,'t',x0),6));
endend输入：clearf=Chebyshev('1/(2-x)',6)输出：f = .8941*t+.117457*t^2+.*t^3+.*t^4+.*t^5输入：clearf=Chebyshev('1/(2-x)',6,'z')输出：f = .500477珐姬粹肯诔厩达询惮墨+.258941*z+.117457*z^2+.*z^3+.*z^4+.*z^5
提问者评价
谢谢啦，我已经知道问题了，我的版本用的是2012b的，其中c的定义应该改成=0*他，而不是0.0。
其他类似问题
为您推荐：
切比雪夫的相关知识
等待您来回答
下载知道APP
随时随地咨询
出门在外也不愁怎么用matlab编辑一个微分方程组,比如：dx=y-0.5x dy=z-x dz=-x+sgn(x);其中sgn(x)是符号函数用ode45编辑的程序是怎么样的？谢谢哪个大侠给我说一下
function DX=YD(t,X)x=X(1);y=X(2);z=X(3);dx=y-0.5*x;dy=z-x;dz=-x+sgn(x);DX=[dx,dy,dz];end
嗯，这个我之前已经弄出来了，我现在又遇到一个问题了，比如编写了个M文件：
function dxdt=fun(t,x)
f(x1)=tanh(20*x1+40)+tanh(20*x1)+tanh(20*x1-40);
10*x2-10*f(x1);
然后再命令窗口中输入了：
>> t_end=1000;
>> x0=[0;1;0];
还需要对f(x1)这个函数赋初值，但是我弄的初值都是错误的，你能给我说怎么算这个的初值吗？
初值的含义是微分方程组的初值，只要初值确定了，解曲线就清楚了。f(x) 是定义的一个函数还是中间变量？如果是中间变量则不需要对f(x1)赋值。有点乱，我不知你在算什么？是rkf积分器吗？
为您推荐：
其他类似问题
扫描下载二维码用MATLAB符号工具箱编程求常微分方程的通解_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
用MATLAB符号工具箱编程求常微分方程的通解
上传于||文档简介
&&不​错
阅读已结束，如果下载本文需要使用1下载券
想免费下载本文？
定制HR最喜欢的简历
下载文档到电脑，查找使用更方便
还剩1页未读，继续阅读
定制HR最喜欢的简历
你可能喜欢}