(window.slotbydup=window.slotbydup || []).push({
id: '2081942',
container: s,
size: '1000,60',
display: 'inlay-fix'后使用快捷导航没有帐号？
只需一步，快速开始
查看: 3547|回复: 11
最后登录贡献度36 金元3622 银行存款0 注册时间阅读权限50帖子主题精华0积分506UID2903066
高级玩家, 积分 506, 距离下一级还需 94 积分
贡献度36 金元3622 注册时间积分506精华0帖子
11:50 上传
11:50 上传
就这两个任务啦，给那个什么果扎带刀鞘，按照提示找到裂谷城家里，就看到一个提示标签，什么都没有了
怎么弄啊？？？？无解~~~
还有那个给戴尔芬找追随者，不知道怎么开始啊~~~也没有提示。上哪去找z追随者呢？？？？
最后登录贡献度36 金元3622 银行存款0 注册时间阅读权限50帖子主题精华0积分506UID2903066
高级玩家, 积分 506, 距离下一级还需 94 积分
贡献度36 金元3622 注册时间积分506精华0帖子
还有一个就是芬恩的鲁特琴也是，好像无解~~~
最后登录贡献度36 金元3622 银行存款0 注册时间阅读权限50帖子主题精华0积分506UID2903066
高级玩家, 积分 506, 距离下一级还需 94 积分
贡献度36 金元3622 注册时间积分506精华0帖子
求做个任务的达人解答的~~~
最后登录贡献度2043 金元88929 银行存款0 注册时间阅读权限90帖子主题精华0积分17065UID389973
资深玩家, 积分 17065, 距离下一级还需 935 积分
贡献度2043 金元88929 注册时间积分17065精华0帖子
那个追随者我也是一头雾水啊~~求任务帝啊~~
最后登录贡献度48 金元8624 银行存款0 注册时间阅读权限90帖子主题精华0积分1054UID3174293
游戏狂人, 积分 1054, 距离下一级还需 946 积分
贡献度48 金元8624 注册时间积分1054精华0帖子
芬恩的鲁特琴 应该是 你接任务之前进过那个地洞 开了箱子 拿掉了任务物品& &接了任务后去
箱子里面没有 完成不了任务&&
追随者的任务 你随便送个跟班就行了 神马 酒馆的雇佣兵 买房子送的 女仆之类的都行
最后登录贡献度36 金元3622 银行存款0 注册时间阅读权限50帖子主题精华0积分506UID2903066
高级玩家, 积分 506, 距离下一级还需 94 积分
贡献度36 金元3622 注册时间积分506精华0帖子
shu3129230
那任务就是卡住了~~~~
还有那个果扎的，好像有也是卡住了~~~
最后登录贡献度0 金元3575 银行存款0 注册时间阅读权限50帖子主题精华0积分358UID206798
高级玩家, 积分 358, 距离下一级还需 242 积分
贡献度0 金元3575 注册时间积分358精华0帖子
我也有一些失败的任务。等开2周目再研究吧
最后登录贡献度18 金元6369 银行存款0 注册时间阅读权限70帖子主题精华0积分709UID2769561
超级玩家, 积分 709, 距离下一级还需 291 积分
贡献度18 金元6369 注册时间积分709精华0帖子
带追随者的话 就是带个跟班去交任务就行了
那个剑鞘的话 可能是卡住了 本来应该是指一个地方你去找的&&你大概是以前就拿了 然后放家里了吧&&这种以前拿了 卡住的任务有好几个的说&&我也卡了几个
最后登录贡献度52 金元6049 银行存款0 注册时间阅读权限70帖子主题精华0积分813UID1259969
超级玩家, 积分 813, 距离下一级还需 187 积分
贡献度52 金元6049 注册时间积分813精华0帖子
本帖最后由 zmx2005war 于
16:08 编辑
还有一个就是芬恩的鲁特琴也是，好像无解~~~
<font color="#76852 发表于
分恩的鲁特琴。。应该是在外面要找的。。我曾经在一个洞洞里看到过- -但是纯当路过了没有拿。。现在后悔了=0=。。话说不管去山洞或地下城之前有没有任务，任务物品大部分都在那里放着的，拿了之后，就是接了任务就直接多出来一个选项：我已经拿到xxx了，直接就可以交的。一般不会bug，琴如果找不到了。。我怀疑你有任务之前去过拿到了，然后丢给商人了- -
那个xxx的剑鞘。。是本书，不是剑鞘。。给任务的时候有说的，因为果扎的徒弟太笨了，她怀疑是没有好好看教程课本。。所以让你去找这本书
最后登录贡献度36 金元3622 银行存款0 注册时间阅读权限50帖子主题精华0积分506UID2903066
高级玩家, 积分 506, 距离下一级还需 94 积分
贡献度36 金元3622 注册时间积分506精华0帖子
zmx2005war
& & 琴在身上，任务物品卖不掉的~~~回去找NPC对话确没有对话选择，估计是卡了
最后登录贡献度26 金元1356 银行存款0 注册时间阅读权限0帖子主题精华0积分240UID1967463
头像被屏蔽
贡献度26 金元1356 注册时间积分240精华0帖子
提示: 作者被禁止或删除 内容自动屏蔽
最后登录贡献度18 金元3698 银行存款0 注册时间阅读权限50帖子主题精华0积分442UID2842808
高级玩家, 积分 442, 距离下一级还需 158 积分
贡献度18 金元3698 注册时间积分442精华0帖子
coc qasmoke你们懂得
Powered by查看: 2585|回复: 9|关注: 0
matlab求解微分方程组，解是1*1
sym量，怎么解决？
[gix,giy,pix,piy]=dsolve('Dgix=-a*zix+m*pix','Dgiy=-a*ziy+m*piy',...
& & 'Dpix=-c*pix+k*fzix','Dpiy=-c*piy+k*fziy','t');
其中zix、ziy、fzix、fziy 都是关于(gix,giy)的函数。
a、m、c、k都是常数。
上述为四个常微分方程组成的方程组，想解出gix和giy（应该也能解出pix,piy但是这个数据不是很重要），但是运行之后得到的是gix (1*1 sym)
giy(1*1 sym) pix(1*1 sym)&&piy(1*1 sym)，
而且 我想画出gix 和giy的图像，使用plot(gix,giy)会报错：Conversion to double from sym is not possible.
关注者： 3
double(gix)
double(gix)
谢谢回复！不过问题还是没有解决，1*1sym是不是没有解出来？现将源程序贴出如下：
k=52.2443;c=77.8411;b=39.8006;a=8.5688;m=4.7916;
gix=rand(10,1);giy=rand(10,1);
plot(gix,giy,'.');
zix=4.*(gix.^2+giy(ii).^2-1).*
ziy=4.*(gix.^2+giy.^2-1).*
plot(zix,ziy,'+');
fzix=(1-exp(-zix))/(1+exp(-zix));
fziy=(1-exp(-ziy))/(1+exp(-ziy));
plot(fzix,fziy,'*');
[gix,giy,pix,piy]=dsolve('Dgix=-a*zix+m*pix','Dgiy=-a*ziy+m*piy',...
& & 'Dpix=-c*pix+k*fzix','Dpiy=-c*piy+k*fziy','t');
% x=double(gix);
% y=double(giy);
% plot(gix,giy,'o');
% [x,y]=dsolve('Dx=4*x-2*y','Dy=x+y','t');
% ezplot(y);
关注者： 54
谢谢回复！不过问题还是没有解决，1*1sym是不是没有解出来？现将源程序贴出如下：
不知道图对否，但是我运行你的脚本是可以出图的(MATLAB2014B，别的版本没试)。
捕获.PNG (42.94 KB, 下载次数: 0)
09:57 上传
不知道图对否，但是我运行你的脚本是可以出图的(MATLAB2014B，别的版本没试)。
谢谢！出图肯定是没有问题的，我那样画图就是为了跟踪数据而已，最后我想得到的是gix和giy的方程组的解，你可以看到上面画的图是我前面三个画图函数的图，那个肯定是没问题的，只是后来加上的那句：
[gix,giy,pix,piy]=dsolve('Dgix=-a*zix+m*pix','Dgiy=-a*ziy+m*piy',...
& & 'Dpix=-c*pix+k*fzix','Dpiy=-c*piy+k*fziy','t');
没有办法解出gix,giy。你试试这个脚本：
k=52.2443;c=77.8411;b=39.8006;a=8.5688;m=4.7916;
gix=rand(10,1);giy=rand(10,1);
plot(gix(ii),giy(ii),'.');
zix=4.*(gix(ii).^2+giy(ii).^2-1).*gix(ii);
ziy=4.*(gix(ii).^2+giy(ii).^2-1).*giy(ii);
plot(zix,ziy,'+k');
fzix=(1-exp(-zix))/(1+exp(-zix));
fziy=(1-exp(-ziy))/(1+exp(-ziy));
plot(fzix,fziy,'*r');
[gix,giy,pix,piy]=dsolve('Dgix=-a*zix+m*pix','Dgiy=-a*ziy+m*piy',...
& & 'Dpix=-c*pix+k*fzix','Dpiy=-c*piy+k*fziy','t');
% x=double(gix);
% y=double(giy);
plot(gix,giy,'o');
% [x,y]=dsolve('Dx=4*x-2*y','Dy=x+y','t');
% ezplot(y);
就会发现matlab报错，画图只是前面那三个，方程组结出来的结果画不出来！
大神您怎么看？
关注者： 54
|此回复为最佳答案
谢谢！出图肯定是没有问题的，我那样画图就是为了跟踪数据而已，最后我想得到的是gix和giy的方程组的解， ...
k=52.2443;c=77.8411;b=39.8006;a=8.5688;m=4.7916;
gix=rand(10,1);giy=rand(10,1);
plot(gix(ii),giy(ii),'.');
zix=4.*(gix(ii).^2+giy(ii).^2-1).*gix(ii);
ziy=4.*(gix(ii).^2+giy(ii).^2-1).*giy(ii);
plot(zix,ziy,'+k');
fzix=(1-exp(-zix))/(1+exp(-zix));
fziy=(1-exp(-ziy))/(1+exp(-ziy));
plot(fzix,fziy,'*r');
[gix,giy,pix,piy]=dsolve('Dgix=-a*zix+m*pix','Dgiy=-a*ziy+m*piy',...
& & 'Dpix=-c*pix+k*fzix','Dpiy=-c*piy+k*fziy','t');
因为你方程解出来是这个东西：
-(exp(-c*t)*(m*(C1 + (fzix*k*exp(c*t))/c) - c*exp(c*t)*(C3 - (t*(a*c*zix - fzix*k*m))/c)))/c
======================
要是数值才能plot
%=====gix=============
symvar(gix) %看它的变量
k=52.2443;c=77.8411;b=39.8006;a=8.5688;m=4.7916;
fzix = 0.8;
zix = 0.3;
t = 0:0.1:1;
gix = subs(gix); %代入数值
gix = vpa(gix,2) %化简
plot(gix,'o');
大神！您的回复太给力了！我使用了您给的代码，画出了gix,图中显示gix为一个定值（理论上不应该是一样的）。不过我还是有一些不明白。仿照您的过程，我添加了giy的处理，试图绘制横坐标为gix,纵坐标为giy的图像，代码如下：
symvar(gix);
symvar(giy);%看它的变量
k=52.2443;c=77.8411;b=39.8006;a=8.5688;m=4.7916;
fzix = 0.8;
fziy = 0.8;
zix = 0.3;
t = 0:0.11;
gix = subs(gix); %代入数值
gix = vpa(gix,2) %化简
giy=subs(giy);
giy=vpa(giy,2);
plot(gix,giy,'o');
代码又出现了同样的问题。
我使用了symvar(giy），查看了它的变量，跟gix只有ziy,fziy,C2 不一样(zix,fzix,C1,C3),我给ziy,fziy,C2 随机给了三个值（不知道您当初是怎么给定这几个值的），然后下面都仿照您对gix的处理。但是结果还是不对。
大神，我现在问题很简单，就是要画出给定10个随机坐标点（gix,giy）在经历上述微分方程组之后，这10个随机点会收敛到新的10个点。并将这新的10个点绘制出来，这样对比可以看到他们的轨迹就行。
大神！您再支支招！不胜感激！
关注者： 54
大神！您的回复太给力了！我使用了您给的代码，画出了gix,图中显示gix为一个定值（理论上不应该是一样的 ...
你前面的代码出错的直接原因是giy有别的变量你没赋值，比如C4.
但是，我感觉问题倒不在这里。不好意思，之前只关注你提的问题是作图，所以只看了gix和giy。
clear
clc
%
k=52.2443;c=77.8411;b=39.8006;a=8.5688;m=4.7916;
ii=1:10;
gix=rand(10,1);giy=rand(10,1);
plot(gix(ii),giy(ii),'.');
hold on
zix=4.*(gix(ii).^2+giy(ii).^2-1).*gix(ii);
ziy=4.*(gix(ii).^2+giy(ii).^2-1).*giy(ii);
plot(zix,ziy,'+k');
hold on
fzix=(1-exp(-zix))/(1+exp(-zix));
fziy=(1-exp(-ziy))/(1+exp(-ziy));
plot(fzix,fziy,'*r');
hold on
[gix,giy,pix,piy]=dsolve('Dgix=-a*zix+m*pix','Dgiy=-a*ziy+m*piy',...
& & 'Dpix=-c*pix+k*fzix','Dpiy=-c*piy+k*fziy','t');
%%
gix
giy复制代码gix =
-(exp(-c*t)*(m*(C1 + (fzix*k*exp(c*t))/c) - c*exp(c*t)*(C3 - (t*(a*c*zix - fzix*k*m))/c)))/c
-(exp(-c*t)*(m*(C2 + (fziy*k*exp(c*t))/c) - c*exp(c*t)*(C4 - (t*(a*c*ziy - fziy*k*m))/c)))/c
理论上到这一步，只要给出积分常数C1，C2，C3，C4的值，就可以求gix,giy的值作图了。方法跟之前差不多，所以我之前没有看前面的代码，是随意给变量赋值的。
但实际上，即使给C1234赋值了，gix是求不出来的。
比如，gix，
gix = -(exp(-c*t)*(m*(C1 + (fzix*k*exp(c*t))/c) - c*exp(c*t)*(C3 - (t*(a*c*zix - fzix*k*m))/c)))/c
红色这部分，根据之前代码求出来的a,c,zix,fzix,k,m的值，
&& a*c*zix
&&629.1709
&&454.8579
&&416.5801
&&255.5145
&&-78.6813
& &47.9276
& fzix*k*m 的结果是个10*10的矩阵，
它俩相减就会出现维数不对的错误。
所以，
C1 = 0.6;
C2=0.4;
C3 = 0.9;
C4 = 0.2;
%%
t = 0:0.11;
gix_data = subs(gix); %代入数值复制代码Error using mupadmex
Error in MuPAD command: NEW arrays must be scalar or the same shape.
这个我现在也觉得比较奇怪，为什么求解出来的是这样。
可能dsolve对矩阵类型的参数有限制？还得去查查。
你前面的代码出错的直接原因是giy有别的变量你没赋值，比如C4.
但是，我感觉问题倒不在这里。不好意思， ...
您做事真是认真！谢谢！这个我自己想了一下，细节问题可能也存在。但是因为我知道这个代码的背景，你没有这方面信息，所以可能我在一条错误的路上请您指点呢！ 所以，我自己也重新认识了一下常微分方程组求解的问题。很多微分方程是不可求解析解的，dsolve是求解析解的matlab函数。我改为数值解（近似解）的函数ode45问题已经解决了。能够画出gix和giy的图像了。
同时再一次感谢您的帮助！
你前面的代码出错的直接原因是giy有别的变量你没赋值，比如C4.
但是，我感觉问题倒不在这里。不好意思， ...
柚耔大神！再追问您一个问题：比如画出 y=2x+1 (1&x&10)的曲线，我很关心它的起始点和终点，如何设置使得它的起始点为&*&终点为&o&,而它的中间轨迹是实线。大神请您不吝赐教！
站长推荐 /1
Powered by}