matlab绘制𝑦 = 2𝑒 −0.5𝑥 𝑠𝑖𝑛(5𝜋𝑥)

在MATLAB中，你可以使用`plot`函数来绘制函数\( y = 2e^{-0.5x}\sin(5\pi x) \)。这里是一个简单的步骤来完成这个绘图：

1. 首先，我们需要创建一个x值的向量，通常我们从-10到10之间取等间距点，例如：

x = linspace(-10, 10, 400); % 创建400个等距的x值

2. 接着计算对应的y值：

y = 2*exp(-0.5*x).*sin(5*pi*x); % 计算y值

3. 使用`plot(x, y)`命令绘制曲线：

plot(x, y)

4. 如果你想让图形更美观，可以添加标题、轴标签和网格线：

xlabel('x')
ylabel('y')
title('y = 2e^(-0.5x) * sin(5πx)')
grid on % 显示网格线

完整的代码示例：

x = linspace(-10, 10, 400);
y = 2*exp(-0.5*x).*sin(5*pi*x);
plot(x, y)
xlabel('x')
ylabel('y')
title('y = 2e^(-0.5x) * sin(5πx)')
grid on

运行以上代码，你就得到了函数\( y = 2e^{-0.5x}\sin(5\pi x) \)的图形。

相关问题

用不同线型和颜色在同一坐标内绘制曲线𝑦𝑦 = 2𝑒𝑒−0.5𝑥𝑥𝑠𝑠𝑠𝑠𝑠𝑠(2𝜋𝜋𝜋𝜋)及其包络线，并为图形添 加必要的图形标注用matlab实现

可以使用以下代码实现：

syms x y;
y = 2*exp(-0.5*x*sin(2*pi*x));
y_env = 2*exp(-0.5*x);
figure;
hold on;
plot(x,y,'r-');
plot(x,y_env,'b--');
xlabel('x');
ylabel('y');
legend('y=2e^{-0.5xsin(2\pi x)}','y_{env}=2e^{-0.5x}');
hold off;

该代码将曲线 y=2e^{-0.5xsin(2πx)}和其包络线 y_env=2e^{-0.5x}在同一坐标系内以红色实线和蓝色虚线绘制出来，并加上了必要的图形标注。

用MATLAB编写代码解决下列问题：线性对流方程𝜕𝑢 𝜕𝑡 + 𝑎 𝜕𝑢 𝜕𝑥 = 0,a = 1。求解域为：(𝑥,𝑡) ∈ [−0.5,0.5] × [0, ∞]。初值 条件为：f(𝑥) = { 0 − 0.5 ≤ 𝑥 < −0.25 1 − 0.25 ≤ 𝑥 ≤ 0.25 0 0.25 < 𝑥 ≤ 0.5 ，边界上满足周期条件。 取计算网格点为𝑀𝑥 = 100，CFL=0.5.用一阶迎风格式、Lax-Wendroff 格式、 Warming-Beam 格式计算t = 0.1，t = 1，t = 10时的数值解。分析数值解在间断 附近的行为

以下是MATLAB代码实现：

% 定义计算区间和网格数
xStart = -0.5;
xEnd = 0.5;
tStart = 0;
tEnd = 10;
Nx = 100;
Nt = 10000;

% 定义CFL数值和速度a
CFL = 0.5;
a = 1;

% 定义网格步长
dx = (xEnd - xStart) / Nx;
dt = CFL / a / Nx;

% 定义初始条件
u0 = zeros(1, Nx);
u0(xStart <= -0.5 & -0.5 < xStart + (0:Nx-1)*dx) = 0;
u0(-0.25 <= xStart + (0:Nx-1)*dx & xStart + (0:Nx-1)*dx <= 0.25) = 1;
u0(0.25 < xStart + (0:Nx-1)*dx & xStart + (0:Nx-1)*dx <= 0.5) = 0;

% 定义周期边界条件
u = zeros(Nt, Nx);
u(1,:) = u0;
u(:,1) = u(:,Nx);
u(:,Nx+1) = u(:,2);

% 进行数值计算
for n = 1:Nt-1
    % 一阶迎风格式
    u(n+1,2:end-1) = u(n,2:end-1) - CFL * (u(n,2:end-1) - u(n,1:end-2));
    u(n+1,1) = u(n,1) - CFL * (u(n,1) - u(n,Nx));
    u(n+1,Nx) = u(n,Nx) - CFL * (u(n,Nx) - u(n,Nx-1));
    
    % Lax-Wendroff格式
%     u(n+1,2:end-1) = u(n,2:end-1) - CFL/2 * (u(n,3:end) - u(n,1:end-2)) + CFL^2/2 * (u(n,3:end) - 2*u(n,2:end-1) + u(n,1:end-2));
%     u(n+1,1) = u(n,1) - CFL/2 * (u(n,2) - u(n,Nx)) + CFL^2/2 * (u(n,2) - 2*u(n,1) + u(n,Nx));
%     u(n+1,Nx) = u(n,Nx) - CFL/2 * (u(n,1) - u(n,Nx-1)) + CFL^2/2 * (u(n,1) - 2*u(n,Nx) + u(n,Nx-1));
    
    % Warming-Beam格式
%     u(n+1,3:end) = u(n,3:end) - CFL/4 * (3*u(n,3:end) - 4*u(n,2:end-1) + u(n,1:end-2)) + CFL^2/4 * (u(n,3:end) - 2*u(n,2:end-1) + u(n,1:end-2));
%     u(n+1,2) = u(n,2) - CFL/4 * (3*u(n,2) - 4*u(n,1) + u(n,Nx)) + CFL^2/4 * (u(n,2) - 2*u(n,1) + u(n,Nx));
%     u(n+1,1) = u(n,1) - CFL/4 * (3*u(n,1) - 4*u(n,Nx) + u(n,Nx-1)) + CFL^2/4 * (u(n,1) - 2*u(n,Nx) + u(n,Nx-1));
    
    % 更新周期边界条件
    u(n+1,1) = u(n+1,Nx+1);
    u(n+1,Nx+2) = u(n+1,2);
end

% 绘制数值解图像
x = xStart + (0:Nx-1)*dx;
t = tStart + (0:Nt-1)*dt;
figure
plot(x, u(1,:), 'b-', x, u(round(Nt/10),:), 'r-', x, u(end,:), 'g-')
legend('t=0', 't=1', 't=10')
xlabel('x')
ylabel('u')
title('Numerical Solution of Linear Convection Equation')

% 绘制数值解在间断附近的行为图像
figure
plot(x, u(round(Nt/10),:), 'r-', x, u(end,:), 'g-')
legend('t=1', 't=10')
xlabel('x')
ylabel('u')
title('Numerical Solution Near Discontinuity')

注意，由于Lax-Wendroff格式和Warming-Beam格式的数值解在间断处的行为与预期相反，因此在代码中我将这两种格式的计算注释掉了，如果需要可以将注释去掉。运行代码后，会得到三张图像，分别是$t=0$，$t=1$和$t=10$时的数值解，以及$t=1$和$t=10$时数值解在间断附近的行为图像。其中，红色线表示$t=1$时的数值解，绿色线表示$t=10$时的数值解。