| 查看: 2884 | 回复: 34 | ||||
| 本帖产生 1 个 数学EPI ,点击这里进行查看 | ||||
[求助]
求解二元二次方程组
|
||||
|
x-y=A(1-x)^2 y^2=2B(x-y)(1-y)基中A,B是已知的,且A,B相对应有100组数,所以我想在excel中计算,算一组就自动得出其他!!求高人指点!如果不能用excel应该怎么处理! |
» 猜你喜欢
26申博自荐
已经有6人回复
最近几年招的学生写论文不引自己组发的文章
已经有9人回复
带资进组求博导收留
已经有8人回复
A期刊撤稿
已经有4人回复
» 本主题相关价值贴推荐,对您同样有帮助:
三元二次的三角方程如何用MATLAB求解??
已经有10人回复
matlab解偏微分方程求助
已经有12人回复
求解二元原函数
已经有3人回复
求二元回归方程的做法
已经有7人回复
偏微分方程阶次
已经有3人回复
Matlab求解二阶偏微分方程组,希望能给出相应的求解代码
已经有10人回复
微分方程当做一元二次方程来解
已经有5人回复
什么语言都行,求一个二元二次方程,实在不会计算机语言。。
已经有6人回复
有限差分离散二阶偏微分方程后如何编译matlab程序解方程组?
已经有3人回复
求助:如何用MATLAB的偏微分方程工具箱画出不同区域内的温度分布图?
已经有5人回复
求助——求解常系数二次微分方程
已经有6人回复
matlab代码求解水动力学方程看不懂,请高手指点
已经有4人回复
【求助】求助一个二元二次回归方程的解法和答案
已经有11人回复
【求助】用MATLAB求解偏微分方程的数值解 悬赏20
已经有8人回复
【求助】求助非线性偏微分方程求解
已经有3人回复
【求助】求助用matlab解偏微分方程
已经有15人回复
【求助】利用matlab求解微分方程得出解函数中含0的意思
已经有3人回复
wurongjun
专家顾问 (职业作家)
-

专家经验: +831 - 数学EPI: 9
- 应助: 791 (博后)
- 贵宾: 0.308
- 金币: 24609
- 散金: 310
- 红花: 75
- 帖子: 3004
- 在线: 881.2小时
- 虫号: 1368482
- 注册: 2011-08-14
- 性别: GG
- 专业: 计算数学与科学工程计算
- 管辖: 数学

2楼2013-10-02 16:58:55
wurongjun
专家顾问 (职业作家)
-

专家经验: +831 - 数学EPI: 9
- 应助: 791 (博后)
- 贵宾: 0.308
- 金币: 24609
- 散金: 310
- 红花: 75
- 帖子: 3004
- 在线: 881.2小时
- 虫号: 1368482
- 注册: 2011-08-14
- 性别: GG
- 专业: 计算数学与科学工程计算
- 管辖: 数学

4楼2013-10-02 22:32:12
wurongjun
专家顾问 (职业作家)
-

专家经验: +831 - 数学EPI: 9
- 应助: 791 (博后)
- 贵宾: 0.308
- 金币: 24609
- 散金: 310
- 红花: 75
- 帖子: 3004
- 在线: 881.2小时
- 虫号: 1368482
- 注册: 2011-08-14
- 性别: GG
- 专业: 计算数学与科学工程计算
- 管辖: 数学
|
程序里面的AB.txt文件是你的数据! function fcqiujie AB=textread('AB.txt'); N=size(AB,1); A=AB(:,1);B=AB(:,2); tx=zeros(N,4); for i=1:N tx(i, =roots([(A(i)^2-2*B(i)*A(i)^2),(-2*A(i)+2*B(i)*A(i)),(-2*A(i)+1),2,1])';end x=tx+1; for i=1:N y(i, =x(i, -A(i)*(1-x(i, ).^2;end for i=1:N idx=find(imag(x(i, )==0&real(x(i, )<1);idy=find(imag(y(i, )==0&real(y(i, )<1);id=intersect(idx,idy); if length(id)==1 disp(['x(',num2str(i),')=',num2str(x(i,id),7),' ','y(',num2str(i),')=',num2str(y(i,id),7),]); else disp(['x(',num2str(i),')=',num2str(x(i,id),7),' ','y(',num2str(i),')=',num2str(y(i,id),7),]); end end |

28楼2013-10-05 13:49:46
wurongjun
专家顾问 (职业作家)
-

专家经验: +831 - 数学EPI: 9
- 应助: 791 (博后)
- 贵宾: 0.308
- 金币: 24609
- 散金: 310
- 红花: 75
- 帖子: 3004
- 在线: 881.2小时
- 虫号: 1368482
- 注册: 2011-08-14
- 性别: GG
- 专业: 计算数学与科学工程计算
- 管辖: 数学
【答案】应助回帖
|
function fcqiujie AB=textread('AB.txt'); N=size(AB,1); A=AB(:,1);B=AB(:,2); tx=zeros(N,4); for i=1:N tx(i, =roots([(A(i)^2-2*B(i)*A(i)^2),(-2*A(i)+2*B(i)*A(i)),(-2*A(i)+1),2,1])';end x=tx+1; for i=1:N y(i, =x(i, -A(i)*(1-x(i, ).^2;end for i=1:N idx=find(imag(x(i, )==0&real(x(i, )<1);idy=find(imag(y(i, )==0&real(y(i, )<1);id=intersect(idx,idy); if length(id)==1 disp(['x(',num2str(i),')=',num2str(x(i,id),7),' ','y(',num2str(i),')=',num2str(y(i,id),7),]); else disp(['x(',num2str(i),')=',num2str(x(i,id),7),' ','y(',num2str(i),')=',num2str(y(i,id),7),]); end end |

29楼2013-10-05 13:52:58
wurongjun
专家顾问 (职业作家)
-

专家经验: +831 - 数学EPI: 9
- 应助: 791 (博后)
- 贵宾: 0.308
- 金币: 24609
- 散金: 310
- 红花: 75
- 帖子: 3004
- 在线: 881.2小时
- 虫号: 1368482
- 注册: 2011-08-14
- 性别: GG
- 专业: 计算数学与科学工程计算
- 管辖: 数学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
376981713: 金币+50, ★★★很有帮助 2013-10-08 17:02:27
376981713: 金币+50, ★★★很有帮助 2013-10-08 17:02:27
|
x(1)=0.9997559 0.9998253 y(1)=-0.8649408 0.04513399 x(2)=0.9997509 0.9998288 y(2)=-0.941755 0.08290534 x(3)=0.9997463 0.9998319 y(3)=-1.013453 0.1155603 x(4)=0.9997421 0.9998346 y(4)=-1.081342 0.1444069 x(5)=0.9997381 0.9998372 y(5)=-1.146264 0.1702852 x(6)=0.9997343 0.9998395 y(6)=-1.208794 0.193772 x(7)=0.9997307 0.9998416 y(7)=-1.269348 0.215283 x(8)=0.9997272 0.9998436 y(8)=-1.328237 0.2351294 x(9)=0.9997239 0.9998455 y(9)=-1.385702 0.2535512 x(10)=0.9997206 0.9998473 y(10)=-1.441932 0.2707376 x(11)=0.9997175 0.999849 y(11)=-1.497078 0.2868413 x(12)=0.9997144 0.9998506 y(12)=-1.551267 0.3019872 x(13)=0.9997115 0.9998522 y(13)=-1.604602 0.3162789 x(14)=0.9997086 0.9998536 y(14)=-1.65717 0.3298037 x(15)=0.9997057 0.9998551 y(15)=-1.709046 0.3426356 x(16)=0.999703 0.9998564 y(16)=-1.760292 0.3548384 x(17)=0.9997003 0.9998577 y(17)=-1.810963 0.3664669 x(18)=0.9996976 0.999859 y(18)=-1.861109 0.3775692 x(19)=0.999695 0.9998602 y(19)=-1.91077 0.3881872 x(20)=0.9996924 0.9998613 y(20)=-1.959985 0.3983582 x(21)=0.9996802 0.9998666 y(21)=-2.20039 0.4435457 x(22)=0.9996687 0.9998712 y(22)=-2.43336 0.4812982 x(23)=0.9996579 0.9998753 y(23)=-2.660799 0.5135192 x(24)=0.9996477 0.9998789 y(24)=-2.883966 0.5414669 x(25)=0.9996378 0.9998822 y(25)=-3.103737 0.5660189 x(26)=0.9996284 0.9998852 y(26)=-3.320749 0.5878119 x(27)=0.9996192 0.999888 y(27)=-3.535481 0.6073231 x(28)=0.9996104 0.9998905 y(28)=-3.748297 0.6249193 x(29)=0.9996018 0.9998929 y(29)=-3.959487 0.6408881 x(30)=0.9995935 0.9998951 y(30)=-4.169279 0.6554595 x(31)=0.9995854 0.9998971 y(31)=-4.377861 0.6688198 x(32)=0.9995775 0.9998991 y(32)=-4.585385 0.681122 x(33)=0.9995697 0.9999009 y(33)=-4.791979 0.6924936 x(34)=0.9995621 0.9999026 y(34)=-4.997748 0.7030411 x(35)=0.9995547 0.9999042 y(35)=-5.202785 0.712855 x(36)=0.9995474 0.9999057 y(36)=-5.407165 0.7220124 |

31楼2013-10-08 16:21:16
3楼2013-10-02 17:27:44
5楼2013-10-03 17:10:40
6楼2013-10-03 17:33:41
feixiaolin
荣誉版主 (文坛精英)
-

专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3430
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
7楼2013-10-03 19:46:38
feixiaolin
荣誉版主 (文坛精英)
-

专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3430
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
8楼2013-10-03 19:48:42
9楼2013-10-04 09:38:48
feixiaolin
荣誉版主 (文坛精英)
-

专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3430
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
376981713: 金币+20, ★有帮助 2013-10-05 13:49:06
感谢参与,应助指数 +1
376981713: 金币+20, ★有帮助 2013-10-05 13:49:06
|
fid = fopen('路径及文件名','r'); [AB,COUNT] = fscanf(fid,'%5f',[2,50]); %% eg. 2列50行 参数A A=AB(1, ;参数B xB=AB(2, ;for k=1:50 %% 1+2*(x-1)+(-2*A+1)*(x-1)^2+(-2*A+2*B*A)*(x-1)^3+(A^2-2*B*A^2)*(x-1)^4=0 C4=A^2-2*B*A^2; C3=-2*A+2*B*A; C2=-2*A+1; C1=2; C0=1; C=[C4,C3,C2,C1,C0]; Root=(roots(C)); %% 打印或显示: x=1+Root(4); x=1+Root(3); x=1+Root(2); x=1+Root(1); end; |
10楼2013-10-04 15:58:31













回复此楼
=roots([(A(i)^2-2*B(i)*A(i)^2),(-2*A(i)+2*B(i)*A(i)),(-2*A(i)+1),2,1])';