matlab将求解sin隐式解,Matlab隐式符号方程求解和赋值

本文介绍了一种利用Matlab求解隐式方程的方法,并通过一个具体的例子展示了如何将符号解转换为数值解的过程。针对一个复杂的隐式方程,作者首先尝试使用Maple进行求解但结果过于复杂,最终选择Matlab完成目标方程y=f(theta)的求解。

近日处理了一个隐式方程的求解,由于方程含有较多的未知数,而且这些参数均是跟实验相关的一些参数,所以,必须得到需要求解的解与

这些参数之间的一个表达式。之前是考虑用的Maple推导求解了该隐私方程,求解结果非常繁杂,不便于编程计算,于是就考虑用Matlab的

符号运算,将解得到之后,再将符号表达是化为数值形式,以便编程计算。

需要求解的目标隐函数为:

0818b9ca8b590ca3270a3433284dd417.png

这里,我们需要求解y关于theta的一个表达式,即y= f( theta )

这里,我在Matlab中编写了一个函数来求解该方程,并且将a,b,c,H,theta作为参数传入该函数,将方程的根的数值解返回,

具体代码实现如下:

function [y]=find_y(a,b,c,H,theta)

ans=solve('sqrt((a*(cos(theta))^2+b*(sin(theta))^2)^2-y^2)-(a*(cos(theta))^2+b*(sin(theta))^2)=sqrt(c^2-(y*sin(theta))^2)-H-c','y');

F=@(a,b,c,H,theta)eval(ans);

y=F(a,b,c,H,theta);

这里用到了eval函数,关于这个函数的具体细节,大家可以在Matlab的帮助文档中进行更多的了解。

运行的主程序为:

clear all

close all

clc

a=100;

b=2;

c=3;

H=10;

theta=0;

y=find_y(a,b,c,H,theta)

计算结果为:

0818b9ca8b590ca3270a3433284dd417.png

评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符  | 博主筛选后可见
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值