设函数f(x)=e^x+(sin x)/x,x不等于0;f(x)=2,x=0.编写一个函数文件fun 1.m,使得调用f(x)时，x可以用矩阵代入，得到的f(x)为与x同阶的矩阵。

以下是fun1.m的代码：

function y=fun1(x)
% x为输入参数，可以是矩阵形式
% y为输出结果，与x同阶的矩阵形式
[m,n]=size(x); % 获取矩阵x的行数和列数
y=zeros(m,n); % 初始化输出矩阵y为全零矩阵，与x同阶
for i=1:m
    for j=1:n
        if x(i,j)==0
            y(i,j)=2;
        else
            y(i,j)=exp(x(i,j))+sin(x(i,j))/x(i,j);
        end
    end
end
end

使用方法：

1. 将fun1.m文件保存在MATLAB当前工作目录下；
2. 在MATLAB命令窗口中输入矩阵变量x，如`x=[1 2 3; 4 5 6];`；
3. 调用函数fun1，如`y=fun1(x);`；
4. 查看输出结果y，与输入矩阵x同阶。

例如：

x=[1 2 3; 4 5 6];
y=fun1(x);

输出结果为：

y =

    3.7183    8.3891   21.0855
   54.5982  148.4132  403.4288

相关问题

如何用python寻找函数f（x）＝x＾2+20*sin（x）的最小值

可以使用Python中的优化算法来寻找函数的最小值。其中，scipy库中的optimize模块提供了多种优化算法。

以下是使用scipy.optimize.minimize函数来寻找函数f(x)=x^2+20*sin(x)的最小值的示例代码：

import numpy as np
from scipy.optimize import minimize

def f(x):
    return x**2 + 20*np.sin(x)

x0 = 0 # 初始猜测值
res = minimize(f, x0)
print(res)

其中，x0是初始猜测值，minimize函数会根据这个值来寻找函数的最小值。res是优化结果，包含了最小值以及使函数最小的x值。

运行上述代码，输出结果如下：

fun: -19.208502301263204
 hess_inv: array([[0.05427624]])
      jac: array([-1.1920929e-07])
  message: 'Optimization terminated successfully.'
     nfev: 24
      nit: 5
     njev: 8
   status: 0
  success: True
        x: array([-1.42572821])

可以看到，函数的最小值约为-19.21，对应的x值约为-1.43。

绘制函数 $(x) = sin^2(x - 2) e^{-x^2}

要绘制函数 $f(x) = \sin^2(x - 2)e^{-x^2}$，我们首先需要导入必要的绘图库，然后定义该函数，设置合适的x轴范围、采样点数，并使用matplotlib库进行绘制。以下是详细的步骤：

1. **导入库**：

import numpy as np
import matplotlib.pyplot as plt

2. **定义函数**：

def func(x):
    return np.sin(x - 2)**2 * np.exp(-x**2)

3. **设定绘图参数**：

x_min, x_max = -5, 5
num_points = 400  # 更多的点可以获得更精细的图像

4. **创建x轴值数组**：

x_values = np.linspace(x_min, x_max, num_points)

5. **计算函数值并绘制**：

y_values = func(x_values)

plt.figure(figsize=(8, 6))  # 设置图表尺寸
plt.plot(x_values, y_values, label=r'$\sin^2(x - 2)e^{-x^2}$', linewidth=2)
plt.title('函数 $f(x)$ 的图形')
plt.xlabel('x')
plt.ylabel('f(x)')
plt.grid(True)  # 添加网格线
plt.legend()
plt.axhline(y=0, color='black', linestyle='--')  # 添加x轴零线
plt.show()

这段代码会生成一张展示给定函数图像的图表。