【求助】求高手帮忙 - MATLAB爱好者论坛-LabFans.com

fansni · 2008-03-21, 23:56

我编了一下程序，求方程组的解
先建立方程函数文件，并保存为myfun.m：
function F=myfun(x)
ne=1.55737476478271;
no=1.54812911343683;
a=0/180*pi;
t=0;
F=[x(1)-no*ne/(no^2*sin(x(2))^2+ne^2*cos(x(2))^2)^(1/2);sin(a)-x(1)*sin(pi/2-t-x(2))]
然后调用优化程序
x0=[1,1];
options=optimset('Display','iter');
[x,fval] = fsolve(@myfun,x0,options）
可以求出解，现在我想画出t与X（1）的关系，也就要改变t值求相应的x（1），不知道这个循环程序应该怎么编，谢谢

fansni · 2008-03-22, 20:38

好可怜啊
都没人理会阿:cry:

2008-03-21, 23:56	#1
fansni 初级会员注册日期: 2008-01-27 年龄: 40 帖子: 3 声望力: 0	【求助】求高手帮忙我编了一下程序，求方程组的解先建立方程函数文件，并保存为myfun.m： function F=myfun(x) ne=1.55737476478271; no=1.54812911343683; a=0/180pi; t=0; F=[x(1)-none/(no^2sin(x(2))^2+ne^2cos(x(2))^2)^(1/2);sin(a)-x(1)*sin(pi/2-t-x(2))] 然后调用优化程序 x0=[1,1]; options=optimset('Display','iter'); [x,fval] = fsolve(@myfun,x0,options）可以求出解，现在我想画出t与X（1）的关系，也就要改变t值求相应的x（1），不知道这个循环程序应该怎么编，谢谢

主题工具
显示可打印版本邮寄本页给好友
显示模式
平板模式切换到混合模式切换到树形模式

相似的主题
主题	主题作者	版面	回复	最后发表
[求助]是gads 方便呢还是gatbx或goat方便?	lizhanglin	MATLAB论坛	1	2008-12-21 17:13
扩频通信的载波恢复方法	winegdsky	MATLAB论坛	2	2008-12-20 01:19
超分辨率图像重建毕业论文	jiangjingjiu	MATLAB论坛	0	2008-04-30 01:05
【求助】plot怎样实现循环，急！急！	ncnc	MATLAB论坛	6	2008-01-23 15:31
【讨论】Matlab动态显示的问题讨论	tqr591513792	MATLAB论坛	0	2007-12-20 14:50

2008-03-22, 20:38	#2
fansni 初级会员注册日期: 2008-01-27 年龄: 40 帖子: 3 声望力: 0	好可怜啊都没人理会阿:cry: