建模作业MAABWord文件下载.docx
- 文档编号:16041609
- 上传时间:2022-11-17
- 格式:DOCX
- 页数:6
- 大小:33.95KB
建模作业MAABWord文件下载.docx
《建模作业MAABWord文件下载.docx》由会员分享,可在线阅读,更多相关《建模作业MAABWord文件下载.docx(6页珍藏版)》请在冰豆网上搜索。
f=sin(x)+cos(x)
2.如何创建符号方程和符号微分方程?
、
创建符号方程
equation=‘a*x^2+b*x+c=0’
创建符号微分方程
diffeq=‘Dy-y=x’
3.符号和数值之间转换都有那几种方式?
符号运算得到的是精确的数值解,有时候要对得到的解析解进行数值转换。
主要有三个函数:
digits,vpa,subs(变量替换函数)。
Digits(D):
函数设置有效数字个数为D的近似解精度.
Vpa(s):
符号表达式s在digits函数设置下的精度的近似解.
Vpa(s,D):
符号表达式s在digits(D)精度下的近似解.
Subs(s,old,new):
将符号表达式中old的换成new变量.
4.如何进行极限、微分、差分计算?
符号极限
limit(f,x,a)%f在x?
a时的极限
limit(f)%f在x?
0时的极限
limit(f,x,a,right)%右极限
limit(f,x,a,left)%左极限
微分函数diff
symsxy;
diff(sin(x^2))%求函数的一阶导数
diff(sin(x*y),x)
diff(sin(x*y),x,2)%关于x求二阶导数
5.如何进行积分运算?
(1)不定积分
f1=int(cos(x))%关于x积分
f2=int(cos(x*t),t)%关于t积分
(2)定积分
x=int(x^2,0,1)
6.如何求解非线性方程组?
(1).非线性方程求解:
solve(‘fun’),求得解析解
x=solve(‘x^2+3*x+4=0’)%得解析解
x1=vpa(x,6)%化为数值解
(2).非线性方程组求解fsolve(‘fun’,x0)
fun由M文件给出函数,x0为初值,是一种迭代解法.
7.如何进行微分方程求解?
求解常微分方程的解析解,用函数
dsolve(‘equation’)
8.编写如下问题的M文件
7.3.1MATLAB中实现符号计算功能有哪几种途径?
(1)通过调用Matlab自己开发的各种功能函数进行常用的符号运算。
如符号矩阵的运算、符号微分方程求解等等。
(2)Matlab语言中的符号计算功能已经很强大了,但为了给一些特殊专业的人员提供方便,Matlab还保留着Maple.m和mpa.m两个函数与Maple的接口,以便实现更多的功能。
3)对那些用惯了计算器的用户来说,Matlab同样是最佳的选择,因为Matlab还提供了符号函数计算器(FunctionCalculator)。
计算器上不超过两个符号函数的基本运算和微积分运算的功能,而且还有函数可视化图形窗。
(1)
,
clear;
clc;
symsx;
f=(x^2-6*x+8)/(x^2-5*x+4);
a=limit(f,x,4);
disp(a);
2/3
(2)
f=abs(x)/x;
a=limit(f,x,0);
NaN极限不存在
(3)
f=(sqrt(1+x^2)-1)/x;
(4)
.
f=x+1/x;
a=limit(f,x,inf);
Inf无穷大
symsxa;
f=1/(x+a);
int(f,x)
ans=log(a+x)
(2)
symsx;
f=(1-3*x)^(1/3);
ans=-(1-3*x)^(4/3)/4
f=1/(sin(x)*cos(x))^2;
ans=-2*cot(2*x)(4)
f=x^2/sqrt(a^2+x^2);
ans=(x*(a^2+x^2)^(1/2))/2-(a^2*log(x+(a^2+x^2)^(1/2)))/2
symsabcx;
x=solve('
a*x^2+b*x+c=0'
'
x'
);
disp(x);
-(b+(b^2-4*a*c)^(1/2))/(2*a)
-(b-(b^2-4*a*c)^(1/2))/(2*a)
2*sin(3*x-pi/4)=1'
(5*pi)/36
(13*pi)/36
sin(x)-sqrt(3)*cos(x)=sqrt
(2)'
-log((2^(1/2)*(-1-1i))/(3^(1/2)+1i))*1i
-log((2^(1/2)*(-1+1i))/(3^(1/2)+1i))*1i
x^2+10*(x-1)*sqrt(x)+14*x+1=0'
(2^(1/2)-1)^2
(17^(1/2)-4)^2
[x,y]=dsolve('
Dx=3*x+4*y'
Dy=5*x-7*y'
x(0)=0,y(0)=1'
t'
)
x=-exp(t*(3*5^(1/2)-2))*((3*5^(1/2))/5+1)*(5^(1/2)/6-1/2)-(5^(1/2)*exp(-t*(3*5^(1/2)+2))*(3*5^(1/2)+5)*((3*5^(1/2))/5-1))/30
y=(5^(1/2)*exp(-t*(3*5^(1/2)+2))*(3*5^(1/2)+5))/30-exp(t*(3*5^(1/2)-2))*(5^(1/2)/6-1/2)
x3.y3+3x3的图形.
ezplot('
5*x^3*y^3+3*x^3'
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 建模 作业 MAAB