求解一个函数三个变量如何运行

求救,各位大佬

怎么运行3000*g-f-F.*cosd(j)这一函数

里面的g f F三个函数都是由arrayfun带入数组得到的,三个数组大小不一

代码如下

本题所用参数如下
空气密度p(km/m^3)
风速v(m/s)
空气阻力系数C
空气密度p1(g/l)
物体迎风面积S(标准小型汽车尺寸m^2)
物体速度v1(km/h)
摩擦力系数M
汽车质量m(kg)
重力加速度g(m/s)
风吹向汽车时的角度(°)
p=1.29
v=(0:0.1:56)
C=0.4
p1=1.293
S=1.6*1.8
v1=(0:1:100)
M=0.018
m=(1000:3000)
g=9.8
j=(0:360)

求不同风速下风力大小

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
可以使用 Matlab 中的 ode45 函数求解三个变量的微分方程。 假设三个变量分别为 x, y, z,微分方程可以写成如下形式: dx/dt = f(x,y,z,t) dy/dt = g(x,y,z,t) dz/dt = h(x,y,z,t) 其中 f(x,y,z,t), g(x,y,z,t), h(x,y,z,t) 分别是 x, y, z 的函数。 以下是求解微分方程的步骤: 1. 定义函数 首先需要定义一个函数,输入参数为 t 和一个三元向量 [x; y; z],输出值为一个三元向量,分别为 dx/dt, dy/dt, dz/dt。函数的定义如下: ```matlab function dxdt = myODE(t, xyz) x = xyz(1); y = xyz(2); z = xyz(3); dxdt = zeros(3,1); dxdt(1) = f(x,y,z,t); dxdt(2) = g(x,y,z,t); dxdt(3) = h(x,y,z,t); end ``` 其中,f(x,y,z,t), g(x,y,z,t), h(x,y,z,t) 是你定义的三个函数。 2. 求解微分方程 使用 ode45 函数求解微分方程,输入参数为定义的函数 myODE,以及一个时间区间,初始状态 [x0; y0; z0]。代码如下: ```matlab tspan = [t0, tf]; xyz0 = [x0; y0; z0]; [t, xyz] = ode45(@myODE, tspan, xyz0); ``` 其中,t0 和 tf 分别是时间区间的起点和终点,x0, y0, z0 分别是初始状态下的 x, y, z 值。 3. 绘制图形 使用 plot 函数绘制 x, y, z 随时间变化的图形,代码如下: ```matlab figure; plot(t, xyz(:,1), 'r-', t, xyz(:,2), 'g-', t, xyz(:,3), 'b-'); xlabel('Time'); ylabel('Value'); legend('x', 'y', 'z'); ``` 其中,xyz(:,1) 表示 xyz 矩阵中第一列的所有元素,即 x 随时间变化的值;xyz(:,2) 表示 y 随时间变化的值,xyz(:,3) 表示 z 随时间变化的值。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值