用matlab程序题：求解下列方程数值解并画出曲线
结果：
代码：
clear allclc
f=@(x,y)([y(2); 0.357*y(1)-0.1905*y(1)*y(2)]);[x,Y]=ode45(f,[0 100],[3 0]);plot(x,Y(:,1),x,Y(:,2)),legend('y','dy/dt')xlabel('x')
如何用matlab求解微分方程并画图，可以先用dsolve（）或ode（）求出其微分方程（组）的解析解或数值解，然后用plot（）绘制其图形。
例如：解微分方程 y'=y-2t/y，y(0)=1，0<t<4
1、用dsolve（）求解，代码及结果如下
2、用ode45（）求解，代码及结果如下
3、当然喽，使用dsolve（）或ode（）求解要根据题意去分析，来决定用那个函数。一般来说，用ode45求解微分方程（组）的数值解用点比较多。
syms y(x)
Dy = diff(y); D2y = diff(y,2);
y = dsolve((1+x^2)*D2y == 2*x*Dy,y(0)==1,Dy(0)==3);
x0 = -1:0.01:1;
y0 = eval(subs(y,x,x0));
plot(x0,y0)
y-x图像
MATLAB程序设计题两道,要求写出计算公式,程序和程序运行结果
clearclc 第一题a=1;%请自定义acmd=['x^2+y^2+',num2str(a),'*x-',num2str(a),'*sqrt(x^2+y^2)'];%将含参数a的方程转化为字符串形式，否则ezplot会无法运行ezplot(cmd,[-3*a,a],[-2*a,2*a])结果为心形线 第二题 A=[2 2 3;1 -1 0;-1 2 1];B=[5;6;8];X=A...
求Matlab大神帮忙做下下面2题
（1）.程序如下：A=[1 5 -3;3 1 2;0 1 5]b=[2 1 -42]'B=[A b];n=4;R_A=rank(A)R_B=rank(B)format rat if R_A==R_B&R_A==n %判断有唯一解 X=A\\b elseif R_A==R_B&R_A<n %判断有无穷解 X=A\\b %求特解 C=null(A,'r') %求AX=0的基础解...
怎么用matlab求这个问题,求编写程序
题主的问题是当扇形的圆心角θ为多大时使该漏斗的容积最大。这个问题我们可以这样来考虑：首先，根据初等数学几何公式，写出V与R、θ的关系式。其二，求V对θ变量的一阶导数，得到V关于θ的表达式。其三，求V对θ变量的二阶导数，得到dV关于dθ的表达式。其四，求d^2V\/dθ^=0的方程，求出θ的值...
关于matlab编程的问题
下面是我对这道题思路：1：对于第一题的微分方程组在matlab中输入[f,g]=dsolve('Df=exp(f*sin(t)+g)','Dg=exp(g*cos(t))+f','t')答案是 Explicit solution could not be found，即没有解析解，因此只能求其近似解，也就是数值解。求解微分方程的数值解需要知道“函数的初值”“自变量...
matlab 程序问题
不知道你的程序是要什么样的结果，单纯从代码来看，就是全角和半角的问题。我已经修改如下：建立M文件eq1．m function dy=eq1(x,y)dy=zeros(2,1);dy(1)=y(2);dy(2)=1\/5*sqrt(1+y(1)^2)\/(1-x);取x0=0，xf=0．9999，建立主程序如下：x0=0;xf=0.9999;[x,y]=ode15s('eq1...
matlab程序问题
不知道你的程序是要什么样的结果，单纯从代码来看，就是全角和半角的问题。我已经修改如下：建立M文件eq1．m function dy=eq1(x,y)dy=zeros(2,1);dy(1)=y(2);dy(2)=1\/5*sqrt(1+y(1)^2)\/(1-x);取x0=0，xf=0．9999，建立主程序如下：x0=0;xf=0.9999;[x,y]=ode15s('eq1...
matlab程序问题。需要用到蒙特卡洛方法
你给出的解是正确的 首先假设有编号为1~16的16个球，其中 编号1~8，8个球是红色，那么9~16，8个球是白色 n=1e6; %游戏100万次 A=0;B=0;C=0;D=0;E=0; %得奖统计清零 for i=1:n examp=randperm(16); %随机打乱1~16，16个自然数 num=sum(examp(1:8)<=8); %...
用matlab程序题:求解下列方程数值解并画出曲线
syms y(x)Dy = diff(y); D2y = diff(y,2);y = dsolve((1+x^2)*D2y == 2*x*Dy,y(0)==1,Dy(0)==3);x0 = -1:0.01:1;y0 = eval(subs(y,x,x0));plot(x0,y0)y-x图像
matlab函数程序设计问题
4、程序设计题 编写一个函数文件实现如下功能：输入一自然数x，若x为偶数，则分别计算出s1=1、s2=1+2、s3=1+2+3….sx=1+2+….+x，并把1到x按序存入向量r中，s1到sx存入相量s中，然后利用plot函数画出r与s的数学关系曲线；若x为奇数，则分别计算出s1=1\/1、s2=(1+2)\/2、s3=(1+...
matlab考试题求帮助!谢谢
第4题：可以用符号运算的dsolve直接求解，也可以用ode45的数值方法求解。y=dsolve('D3y+1.2*D2y+1.6*Dy+0.9*y-1=0','y(0)=0,Dy(0)=0,D2y(0)=0');ezplot(y,[0 20])axis auto 第5题：这里有个问题，100分到五分制的转换其实有不同的划分规则，程序只是提现转换思路。function ...
席澜宝龄： function y=zhongy(x) % 保存为 zhongy.m 文件 y=[2*x(1)-x(2)-exp(-x(1));-x(1)+2*x(2)-exp(-x(2))]; ------------------------------------------------------------ [x,fval,exitflag]=fsolve('zhongy',[0.1;0.1]) % x 为所求的解 %fval为解代入方程的验证结果 % exitflag=1表示函数收敛于解;反之,=0
拉孜县13843398694： matlab中用牛顿法求方程的根用matlab求解以下问题:用牛顿法求解方程x=e^ - x在x0=0.5附近的根(ε=10^ - 5) - ？
席澜宝龄：[答案] 代码如下:function rtn=newton1(fx,dfx,x0,tol,N)% Newton Method% The first parameter fx is a external function with respect to viable x.% The second parameter dfx is the first order diffential function ...
拉孜县13843398694： 如何用matlab计算下列方程 - ？
席澜宝龄： >> A=[5 -1 -4;-1 6 -3;-8 -10 23];b=[-1000;400;400];u=A\b u =-275.8621-23.7548-88.8889说明:上面U中的三个数分别是U1,U2,U3
拉孜县13843398694： 用MATLAB解下列联立方程 3x+4y - 7z - 12w=4 5x - 7y+4z+ 2w用MATLAB解下列联立方程 3x+4y - 7z - 12w=45x - 7y+4z+ 2w= - 3x +8z - 5w=9 - 6x+5y - 2z+10w= - ... - ？
席澜宝龄：[答案] A=[3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10]; b=[4;-3;9;-8]; r=rank(A); x=inv(A)*b
拉孜县13843398694： 求救!期末作业,MATLAB程序设计的题目~~求答案~1、求下列联立方程的解3x+4y - 7z - 12w=45x - 7y+4z+2w= - 3x+8z - 5w=9 - 6x+5y - 2z+10w= - 8 求系数矩阵的秩... - ？
席澜宝龄：[答案] A = [3,4,-7,-12;5,-7,4,2;1,0,8,-5;-6,5,-2,10]; %A为系数矩阵 b = [4;-3;9;-8]; Z = rank(A);%Z为秩 X = A\b;%解
拉孜县13843398694： matlab怎样求解一元二次方程 - ？
席澜宝龄： 一、工具:Matlab2012b 二、操作步骤: A.解一元方程【1】先举一例,解方程＂x^2+100*x+99=0＂在matlab ”Command Window＂中输入如下命令:x=solve('x^2+100*x+99=0','x')见下图【2】回车后,matlab就求出了这个一元二次方程...
拉孜县13843398694： 怎么用matlab解以下方程,写代码 - ？
席澜宝龄： 因为涉及三角函数,往往方程组无解.所以,为得到解,必须把所有常数用数值代入,不要用字母,如L1,L2,L3要用数值,避免解不出.
拉孜县13843398694： 求助:如何用MATLAB解出下面的方程. - ？
席澜宝龄： clc;clear f='2185.65*(1-x^3*(0.25+x^2)^(-3/2))+28.67*((1.1*x+1.21*atan(1.1/x)+x^2*atan(1.1/x))/(1.21+x^2))-1.83*x' ezplot(f,[-100 100]),grid x=fzero(f,9) 结果:x = 9.11523183213223
拉孜县13843398694： 用matlab,lindo,lingo,求解下面方程式, - ？
席澜宝龄： matlab里面有专门的工具箱解决这种问题.lingo中这种问题应该很简单,直接在命令窗口中输入:min=(((25+x^2)^(1/2))5.6)+(21.6542(25+(40/(20-x))^2))+(6((64+(20-x)^2)^(1/2))) x<=20 然后就可以编译得到解.
拉孜县13843398694： 怎么用matlab求解下列方程组 - ？
席澜宝龄： 上面是解的顺序出了点问题 下面应当可以.syms x1 x2 x3 x4 x5 x6 x7 x8 x9 xx0 xx1 [x1,x2,x3,x4,x5,x6,x7,x8,x9,xx0,xx1]=solve('x1*(1+0.03519)=xx1','x2*(1+0.03978*2)=xx1','x3*(1+0.0459*3)=xx1','x4*(1+0.0459*3)*(1+0.03519)=xx1','x5*(1+0.049725*...
