高分求用MATLAB编译的高数或大学物理题(40~50行左右)
%% 画陆地图
subplot(1,2,1);
% 等高线图
x=-1000:5:1000;
y=x;
[x y]=meshgrid(x,y);
z=30*(1-(x.^2+y.^2)/10^6);
[c,h]=contour(x,y,z,10);%画10条等高线
set(h,'ShowText','on');%显示等高线的值
axis('equal');
title('小岛等高线');
%三维图
subplot(1,2,2);
mesh(x,y,z);
title('小岛三维图');
%% 求高潮与低潮时小岛露出水面的面积比。
%由于等高线是圆,故只需求出高低潮时小岛等高线圆面积即可
%涨潮时:z=-2,既有30*(1-(X^2+Y^2)/10^6)=-2
%此等线围成面积 area_h=pi*r^2=pi*(X^2+Y^2)
%故小岛高潮时面积:
area_h=pi*(1+2/30)*10^6;
%同理,低潮时:
area_l=pi*(1-2/30)*10^6;
%故面积比为
area_ratio=area_h/area_l
http://zhidao.baidu.com/question/159115790.html?si=3
clc;clear
a=0;b=1;
fa=1-a-sin(a);
fb=1-b-sin(b);
c=(a+b)/2;
fc=1-c-sin(c);
if fa*fb>0,break,end
while abs(fc)>0.5*10^(-4)
c=(a+b)/2;
fc=1-c-sin(c);
if fb*fc>0
b=c;
fb=fc;
else
a=c;
fa=fc;
end
end
format long
fx=fc,x=c
结果:
fx =
-2.414986223420179e-005
x =
0.510986328125000
精确解:
>> x=solve('1-x-sin(x)')
x =
.51097342938856910952001397114508
clear;
clf;
q=2e-6;
k=9e9;
a=1.5;
b=-1.5;
x=-6:0.6:6;
y=x;
[X,Y]=meshgrid(x,y);
rp=sqrt((X-a).^2+(Y-b).^2);
rm=sqrt((X+a).^2+(Y+b).^2);
V=q*k*(1./rp-1./rm);
[Ex,Ey]=gradient(-V);
AE=sqrt(Ex.^2+Ey.^2);
Ex=Ex./AE;Ey=Ey./AE;
cv=linspace(min(min(V)),max(max(V)),49);
contourf(X,Y,V,cv,'k-')
%axis('square')
title('\fontname{隶书}\fontsize{22}偶极子的场'),hold on
quiver(X,Y,Ex,Ey,0.7)
plot(a,b,'wo',a,b,'w+')
plot(-a,-b,'wo',-a,-b,'w-')
xlabel('x');ylabel('y'),hold off
求助利用MATLAB进行复杂函数的定积分方法! 详情见补充,多谢了!_百度...
如简单的梯形法、辛普生(Simpson)•法、牛顿-柯特斯(Newton-Cotes)法等都是经常采用的方法。它们的基本思想都是将整个积分区间[a,b]分成n个子区间[xi,xi+1],i=1,2,…,n,其中x1=a,xn+1=b。这样求定积分问题就分解为求和问题。2.数值积分的实现方法基于变步长辛普生法,MATLAB给出...
求用matlab计算一个六重积分的具体程序,绘制仿真图
问题 首先,说一下条件给的不清楚或可能有错的地方:1、按照公式中使用的符号,下面的常数d、h应为D、H,a1、a2应为R1、R2;2、公式中,第三重积分的下限确定是2*sqrt(R1^2-x)?【后面有讨论】3、求积分得到的结果就是一个具体的数,标题要求的【绘制仿真图】是什么意思?参考代码 下面给出...
用matlab求解二阶微分方程
用matlab的dsolve函数计算该二阶微分方程的解析值是有一定的难度。只能用数值分析的方法来求解其数值解。当然求其数值解,必须还得已知a、θ的具体数值。该问题可以用matlab的ode45函数求得。求解方法如下:1、建立微分方程的自定义函数,odefun(x,y)2、当a=10、θ=π\/6时,执行下列代码 theta=pi...
matlab如何实现用二分法求代数方程在区间内的解?
1、在MATLAB中,求解符号微分方程通解的指令格式为:y=dsolve('equation','x')%equation指符号微分方程,x为符号变量。2、如:>> syms a bfun='Dy=a*x+b';y=dsolve(fun,'x')。3、符号微分方程的特解y=dsolve('equation','codition','x')%equation为符号微分方程condition为微分方程的定解...
用matlab求差分方程的解?
如何用matlab解差分方程?差分方程是指未知函数差分的方程,其一般形式为 F(t,yt,Δyt,Δy^2t,···,Δy^nt)=0 差分方程的通解是指含有相互独立的任意常数,且任意常数的个数等于方程的价数的解;差分方程的特解是指不含有任意常数的解;为了说明如何用matlab解差分方程问题,特举例:例1...
如何用matlab求解微分方程并画图
如何用matlab求解微分方程并画图,可以先用dsolve()或ode()求出其微分方程(组)的解析解或数值解,然后用plot()绘制其图形。例如:解微分方程 y'=y-2t\/y,y(0)=1,0<t<4 1、用dsolve()求解,代码及结果如下 2、用ode45()求解,代码及结果如下 3、当然喽,使用dsolve()或ode(...
怎么用matlab算含字母的定积分? 比如:
含参量的定积分用matlab是不能直接求得。首先,设定r等于某一定值 然后,求定积分 如:r=100;int(r*cos(x)+sqrt(400-r^2*(sin(x))^2)-2*sqrt(100-r^2*(sin(x))^2),x,0,pi)ans =40*elliptice(25)- 40*elliptice(100)elliptice——椭圆积分 或直接求不定积分 syms x r int(r...
如何用matlab求解二阶微分方程,以及程序实例
如何用matlab求解二阶微分方程,对于一般的微分方程(组)可以用dsolve()函数求得其解析解或数值解,对于较复杂的微分方程(组)可以用ode45()函数求得其数值解。例如:微分方程,用dsolve和ode45计算t=0.1,0.2时y的值。y = dsolve('D2y+0.5*Dy+2*y =0','Dy(0)=0,y(0)=1')t=0....
如何用matlab求解二阶微分方程,以及程序实例
1、对于求数值解的微分方程,你可以用ode45()函数求解。如求下列微分方程 func。m %自定义微分方程的函数 function z = func(t,y)z =[y(2);(1-y(1)^2)*y(2)-y(1)];main。m %主程序 clear all;close all;clc y0 = [0.25;0];h = 0.1;a = 0;b = 20;[t1 y1] = ...
如何利用MatlAB软件用矩阵求逆、矩阵除法、矩阵分解求下面线性方程组的...
A=[5,1,-1,1;1,3,-1,-1;-1,-3,-1,5;1,1,-1,4]C=inv(A) \/*这是求矩阵 A的逆矩阵 B=[1 2 3 4]'C*B得解 D=[A B]rref(D) \/*这是求矩阵D的行简化阶梯矩阵,得出未知数的解 矩阵的除法:A\\B命令
仲长裕血康: 】研究偶极子(Dipole)的电势(Electric potential)和电场强度(Electric field density) clear; clf; q=2e-6; k=9e9; a=1.5; b=-1.5; x=-6:0.6:6; y=x; [X,Y]=meshgrid(x,y); rp=sqrt((X-a).^2+(Y-b).^2); rm=sqrt((X+a).^2+(Y+b).^2); V=q*k*(1./rp-1./rm); [Ex,Ey...
黄石港区17052367737: 急!!!!高分求函数最大化问题.希望得到matlab或者lingo程序,请高手指点. 问题见附件. - ?
仲长裕血康: 用lingo求解的,程序如下:model: max=(2-0.9*x1)*(200-10*x2)*0.4*(1-(3-x1)/0.9*(x2-x1))-0.9*x1*(200-10*x2)*0.04*(1-(3-x1)/0.9*(x2-x1)^2); (200-20*x2+10*x1)*0.04*(0.9*(x2-x1)-3+x1)+0.08*(3-x1)=0; 20/9>x1; x2 3 x1+9*x2>30; end 结果为y的最大值为:513.0279 此时x1为:1.981881x2为:11.00527 希望回答对你有帮助····
黄石港区17052367737: 高分求~数学matlab编程题,求过程?
仲长裕血康: 尝试一下如下的matlab代码: syms x; xlist=[10:10:1000,1500,inf]; % value list of X fun = int(sin(x)/x); % define the integrate, reference : http://www.mathworks.com/help/toolbox/symbolic/int.html ylist=sin(xlist)./xlist; % compute the function value jflist ...
黄石港区17052367737: 求用Matlab命令编写代码计算下列数学算式的值;答案对的话,给高分.在线等,谢谢 - ?
仲长裕血康: clear all;clc; y=((sin(log(5))+log10(108)+64^(1/5))/asin(0.5))*exp(abs(sind(175))) y = 11.1067
黄石港区17052367737: 高分求解一个MATLAB编程计算程序 - ?
仲长裕血康: 利用fsolve求解此方程组:12 f = @(x)[x(1)^2-(1:5).*x(2:end)-1 sum(x(2:end))-10]; fsolve(f, zeros(1, 6));
黄石港区17052367737: 求MATLAB验证高数公式? - ?
仲长裕血康: >> clear>> syms x>> limit(sin(x),x,0),limit((1+1/x)^x,x,inf)ans =0ans =exp(1)验证了两个重要极限,即 limit(sin(x),x,0)=0 和 limit((1+1/x)^x,x,inf)=e 的正确.其中exp(1)即为常数e .
黄石港区17052367737: 超高分数求用极大似然法估计信号频率,实现matlab仿真 - ?
仲长裕血康: 题目给得太模糊了,就算大家想帮你也帮不了啊.首先,SNR=8dB,但是滤波器带宽呢?理论上,噪声的功率是无限的,如果你不给出适合的带宽,仿真无法进行.而且,你的噪声是高斯白噪?还是别的?其次,最大似然只是一种方法,是抽象的.你可以在时域进行最大似然,也可以在频域进行最大似然,这都不定的.在我眼中,估计一个信号的频率,一个FFT就够了,如果要最大似然,也没什么好似然的.你找到FFT结果的最大值,也就是信号的频率,其实也是一种最大似然,不过是频域上的而已.
黄石港区17052367737: matlab中输入20个数,求最大值和最小值,用循环语句和调用MAX函数和MIN函数来实现 - ?
仲长裕血康: 将这20个数用数组A表示,则最大值为max(A),最小值为min(A).具体在matlab命令行中输入doc max和doc min会有更详细的使用说明.另外用for求max(A),记为maxA: maxA=A(1)for i=1:length(A)if(A(i)>maxA)maxA=A(i)end end
黄石港区17052367737: matlab如何求函数的最大值 - ?
仲长裕血康: matlab如何求函数的最大值 (1)先在命令窗口给s,A赋值 (2)>> y=@(x)0.65*(s+x)*x/(A+x) (3)[x,fval]= fminbnd(y,0,inf)%用于求函数y在x在[0,inf]内的最小值 则在取负即为求最大值.不过本题我试验了,本题好像有问题.手算都没有最大值的,因为对y求导后y'在x>0内恒大于零,故y为增函数,所以应为在x取无穷大时,y最大,为无穷大.
黄石港区17052367737: 求一段matlab代码,高分送上 - ?
仲长裕血康: mathematica的解法要不要……在mathematica下很容易做到,说白了就是判你那式子的最大值和最小值:Minimize[{n0/n1 - n0/n3 + n1/n2 - 1,0.1 < n0 < 0.8 && 0.05 < n1 < n0 && 0.05 < n2 < n1 &&0.05 < n3 < n2}, {n1, n2, n0, n3}]解出:{-15., {...