matlab求解五元三角函数方程组哪有？
matlab求解五元三角函数方程组哪有？
08-09-25 &
数据点太少了，用非线性回归。x=[11.96,12.22,19.76,27.56,27.82]y=[1.7,7.3,4.5,1.0,8.0]fun=inline('(a(1)*sin(a(2)*x)+a(3)).*sin(a(4)*x)+a(5)','a','x')a=nlinfit(x,y,fun,[1 2 3 4 5])yy=(a(1)*sin(a(2)*x)+a(3)).*sin(a(4)*x)+a(5)plot(x,y,'o-',x,yy)结果：a =  1.0e+003 *   -0.0119    0.0016    1.3057    0.0040    1.1299yy =    2.7506    7.3000    2.4013    2.0523    7.9972【补充】我找到了另一个函数，把高潮时的数据回归了一下。A=[...0 7710.51.91.12.62.43.53.54.34.15.45.86.46.37.87.68.8.79.39.810.4 71510.11.11.12.12.13.475 70613.14.15.15.16.16.17.17.18.18.19.19.20.20.21.21.22.22.23.23.24.24.24.25.26.26.27.27.28.29.29.30.30.31.31.32.32.33.33.34.34.35.35.36.36.37.37.38.38.39.39.40.40.41.41.42.42.43.43.44.45.45.46.46.47.47.48.48.49.49.50.50.51.51.52.52.53.53.54.54.55.55.56.56.57.57.58.59.025 68459.];x=A(:,1);y=A(:,2);%  f(x) =  a0 + a1*cos(x*w) + b1*sin(x*w) + a2*cos(2*x*w) + b2*sin(2*x*w)fun=inline('a(1) + a(2)*cos(x*a(3)) + a(4)*sin(x*a(3)) + a(5)*cos(2*x*a(3)) + a(6)*sin(2*x*a(3))','a','x')%fun=inline('(a(1)*sin(a(2)*x)+a(3)).*sin(a(4)*x)+a(5)','a','x') a=nlinfit(x,y,fun,[600 40 0.2 12 50 -60]) %yy=(a(1)*sin(a(2)*x)+a(3)).*sin(a(4)*x)+a(5);yy=a(1) + a(2)*cos(x*a(3)) + a(4)*sin(x*a(3)) + a(5)*cos(2*x*a(3)) + a(6)*sin(2*x*a(3));plot(x,y,'o',x,yy)
请登录后再发表评论!查看: 20661|回复: 10|关注: 0
Matlab关于含有三角函数的方程的求解
<h1 style="color:# 麦片财富积分
新手, 积分 5, 距离下一级还需 45 积分
大家好，我是matlab初学者，有一个问题想向大家请教，如果一个方程中含有三角函数，怎么求解呢？比如方程组：x - (x^2 - 16) = 2*y ， 2*cos(x/2) - 1) * sin((x/2 + y)) = (y - x/2)
用matlab解出了x和y的值，但是我不知道是否可用，因为我试算了一下sin（x）=1，结果只有一个，但是我们都知道它有无穷多解，所以我就不知道我计算的结果到底可不可用了。。。
这个问题困扰我很久了，希望大家帮帮我，感激不尽！:handshake
MATLAB 数学、统计与优化版块优秀回答者
关注者： 18
回复 1# 郭林坪 的帖子
[x y]=solve('x-(x^2-16)=2*y','(2*cos(x/2)-1)*sin(x/2 + y)=(y-x/2)','x','y')
你的表达式都有问题，自己改改吧。
<h1 style="color:#4 麦片财富积分
关注者： 22
高中学过数形结合法吧？
到底有几组解，画个函数图看看几个交点不就行了？
另外通过自己扎实的数学推理也应该能够得到解集个数
隐函数图ezplot函数可以绘制
<h1 style="color:# 麦片财富积分
顶三楼:) :) :)
<h1 style="color:# 麦片财富积分
向三楼的Neptune_zx请教：像这样的方程：x - (x^2 - 16) = 2*y ， 2*cos(x/2) - 1) * sin((x/2 + y)) = (y - x/2)怎么用数行结合法呢？
<h1 style="color:# 麦片财富积分
谢谢二楼的zhanglu，matlab中我解的结果与你的式子解的结果是一样的，但是那不是我想要的结果，我想要matlab显示这个方程组的所有的解，而不是一个解。
MATLAB 数学、统计与优化版块优秀回答者
关注者： 18
回复 6# 郭林坪 的帖子
我说了你写的表达式有问题呢。。。自己看看自己写的表达式吧。。我可不知道你要表达什么式子。。
<h1 style="color:# 麦片财富积分
回复 7# zhanglu 的帖子
呵呵，明白你的意思o(∩_∩)o
我是把需要解的方程写出来了(*^__^*)
<h1 style="color:#4 麦片财富积分
关注者： 22
%% 数形结合
range=[-20 20];
set(gcf,'position',[200 200 900 600]);
curve_handle(1)=ezplot('x-(x^2-16)=2*y',range);
curve_handle(2)=ezplot('(2*cos(x/2)-1)*sin(x/2 + y)=(y-x/2)',range);
curve_handle(3)=ezplot('y=x/2',range);
set(curve_handle(1),'color','b');
set(curve_handle(2),'color','r');
set(curve_handle(3),'color','g','linestyle','-');
title('function');
legend('x-(x^2-16)=2*y','(2*cos(x/2)-1)*sin(x/2 + y)=(y-x/2)','y=x/2',...
& & 'location', 'NorthWest');
%通过函数图形可知有且仅有2组解（曲线共2个交点）
%% 利用数值函数fsolve求解
F=@(x)[x(1)-(x(1)^2-16)-2*x(2);
& && & (2*cos(x(1)/2)-1)*sin(x(1)/2 + x(2))-(x(2)-x(1)/2)];
%在交点附近取初始值供fsolve求解，这里2个初始值
S0_1=[-5 -5];
S0_2=[ 5&&5];
%求解结果放在S1和S2中
[S1 Fval1]= fsolve(F,S0_1,optimset('Display','iter')) %#ok&*NOPTS&
[S2 Fval2] = fsolve(F,S0_2,optimset('Display','iter'))
Solution={S1;S2};
%% 数据可视化
plot(S0_1(1),S0_1(2),'.','markersize',15,'color','g')
plot(S0_2(1),S0_2(2),'.','markersize',15,'color','g')
plot(S1(1),S1(2),'.','markersize',15,'color',[0 0.5 0])
plot(S2(1),S2(2),'.','markersize',15,'color',[0 0.5 0])
legend('x-(x^2-16)=2*y','(2*cos(x/2)-1)*sin(x/2 + y)=(y-x/2)','y=x/2',...
& &['S0_1 ','(',num2str(S0_1(1)),',',num2str(S0_1(2)),')'],...
& &['S0_2 ','(',num2str(S0_2(1)),',',num2str(S0_2(2)),')'],...
& &['S1& &','(',num2str(S1(1)),',',num2str(S1(2)),')'],...
& &['S2& &','(',num2str(S2(1)),',',num2str(S2(2)),')'], ...
& &'location', 'NorthWest');
%% 下面程序可加可不加
% annotation(gcf,'arrow',[0.709 0.813],...
%& &&&[0.223 0.896],...
%& &&&'Color',[0.61 0.897 0.904]);
% annotation(gcf,'arrow',[0.581 0.011],...
%& &&&[0.394 0.142],...
%& &&&'Color',[0.61 0.897 0.904]);
untitled.png (23.99 KB, 下载次数: 250)
20:40 上传
[ 本帖最后由 Neptune_zx 于
20:45 编辑 ]
<h1 style="color:# 麦片财富积分
高手啊！！！
太感谢你了！！！绝对的高手啊！！！我QQ号，方便的话加我吧(*^__^*)
我还有问题想向你请教o(∩_∩)o
站长推荐 /3
Powered by求大神帮忙！带有三角函数方程组化简后求解，结果很是诡异？_matlab吧_百度贴吧
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&签到排名：今日本吧第个签到，本吧因你更精彩，明天继续来努力！
本吧签到人数：0成为超级会员，使用一键签到本月漏签0次！成为超级会员，赠送8张补签卡连续签到：天&&累计签到：天超级会员单次开通12个月以上，赠送连续签到卡3张
关注：162,637贴子：
求大神帮忙！带有三角函数方程组化简后求解，结果很是诡异？收藏
1、问题描述cos(a)+sin(b)=1
①sin (a)-cos(b)=y
②已知y为常数，求解a和b化简后：①^2+②^2 得：1+2*sin（b-a）=y^2
③2、对化简前后matlab求解2.1 化简前equ1=&#39;cos(a)+sin(b)=1&#39;;equ2=&#39;sin(a)-cos(b)=y&#39;;s=solve(equ1,equ2,&#39;a,b&#39;);a=subs(s.a)b=subs(s.b)结果：当y=0时，a =
pi/3 -pi/3b =
pi/6 (5*pi)/6检验结果正确2.2 对化简后求解s=solve(&#39;1+2*sin(b-a)=y^2&#39;,&#39;a,b&#39;);s1=s.a;s2=s.b;a=subs(s1)b=subs(s2)结果：当y=0时，a =
pi/6 (5*pi)/6b = 0 0验证后结果显然错误当y=其他值时候，b总是等于03、疑问3.1 化简后为什么会出现这种错误？
低氮锅炉，RBI中国总代理,加拿大原装进口,技术安全可靠,超低氮排放
就没有人愿意解答么
登录百度帐号推荐应用
为兴趣而生，贴吧更懂你。或