1回答
我想要识别图像上车牌的一个白色像素。我不知道如何接受它，因为汽车的颜色比车牌的白色更“白”。我写了一个可运行的代码，但我必须用鼠标选择一个车牌点，然后我找到了车牌，我可以在它周围画一个矩形。你知道如何自动找到车牌上的一个白色像素吗？
代码如下：
clear all;
close all;
I= imread('volvo_frame_0001.bmp');
figure, imshow(I)
level=0.5;
BW=im2bw(I,level);
figure, imshow(BW);
BW2 = bwselect(BW,4);
figure, imshow(BW
浏览 1提问于2013-10-26得票数 1
5回答
用c程序可以求函数的导数吗？我使用matlab是因为它有一个内置的函数diff()，可以用来求一个函数的导数。
f(x)=x^2
用c可以求出上述函数的导数吗?算法是什么？
浏览 2提问于2012-10-16得票数 16
回答已采纳
1回答
求DJ(i)的第n项。
给定DJ(i)是i出现在DJ(i)中的次数
数组DJ(i)如下所示：
i 1 2 3 4 5 6 7 8 9 10 .....
DJ(i) 1 2 2 3 3 4 4 4 5 5 .....
我想要系列DJ(i)的第n个术语。我必须对这个项进行一些进一步的计算。
浏览 0提问于2012-10-14得票数 0
回答已采纳
1回答
矩阵中的回溯求最小元
、、、、
不确定我是否应该在这里张贴，如果不让我知道请！
问题在于在二级结构中的能量最小化的回溯过程，该过程采用动态规划。
所以，我有一个上三角矩阵W，它充满了这个公式：
W(i, j) = min([W(i,j-1),
W(i+1,j),
V(i,j),
min_k // I will use an example, where this argument is NOT used
]);
同时，我还使用以下逻辑填充回溯表B：
if(min was the first option)
浏览 5提问于2014-05-15得票数 0
回答已采纳
1回答
我希望找到包含给定列值的行。第一个()返回预期的3行号，但last语句只返回1，而不是我预期的12 4。
这个例子说明了。
> x <- matrix(c(1,0,NA,1,1,1,0,1,0,1,1,NA),nrow=4,ncol=3)
> x
[,1] [,2] [,3]
[1,] 1 1 0
[2,] 0 1 1
[3,] NA 0 1
[4,] 1 1 NA
> which(!is.na(x[,1]))
[1] 1 2 4
> which(!is.na(x[,1]) &&
浏览 1提问于2014-07-11得票数 0
回答已采纳
1回答
我需要帮助通过使用类获得使用Antoine的临时时间。
我的根本行动失败了，我不知道为什么。
我的代码：
from __future__ import division, print_function
%matplotlib inline
import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import root
class Chemical(object):
def __init__(self, name_1, balance_1, name_2, balance_2):
s
浏览 0提问于2019-02-18得票数 0
回答已采纳
1回答
求圆上点的角度
、、、
我得到了(x，y)点，我知道半径。我想找出圆上每个点的角度。将来我想用这个代码来测试角度是否在区间0,2Pi上均匀分布。
我尝试过的代码是
b <- matrix(c(d[,2],d[,5]),ncol=2)
S <- (b/sqrt((rowSums(b^2))))
#S <- matrix(S[!is.na(S)],ncol=2)
lim <- c(-1,1)
plot(S,xlim=lim,ylim=lim,xlab=expression(S1),ylab=expression(S2))
#S<-S[!is.na(S)]
U <- matrix(c(0
浏览 0提问于2020-05-14得票数 0
3回答
我试图找到无向图的度分布。我尝试了以下代码：
graph = { "a" : ["c"],
"b" : ["c", "e"],
"c" : ["a", "b", "d", "e"],
"d" : ["c"],
"e" : ["c", "b"],
"f&
浏览 1提问于2014-03-16得票数 2
回答已采纳
3回答
如何求对角值
、、
我想要找到数组值的和。数组可能是这样的：
const maxtrix = [
[1, 2, 3, 4],
[5, 10, 7, 5],
[8, 11, 15, 9],
[6, 19, 4, 14],
]
我想在数组中找到类似于这个[4, 7, 11, 6]的值。
我试过这样做，但我无法组织代码。
下面是我写的代码：
const maxtrix = [
[1, 2, 3, 4],
[5, 10, 7, 5],
[8, 11, 15, 9],
[6, 19, 4, 14],
]
// right - 4, 7, 11, 6
let left
浏览 1提问于2018-09-12得票数 0
回答已采纳
我似乎在使用合适的代码(或者一个我可以编辑和利用的基的例子)来使用Eratosthenes方法的筛子来找出质数的和时遇到了麻烦。我知道下面的代码更容易工作，但如果可能的话，我想要掌握更大的代码：
require 'prime'
Prime.each(1000).reduce(:+)
虽然有Eretosthenes筛子和sum的例子，但它们是在c++中的，因此对于Ruby来说是不可靠的。
浏览 0提问于2015-08-20得票数 0
2回答
我需要找到点(X,Y)到由四系数C0, C1, C2, C3定义的曲线的最小距离，比如y = C0 + C1X + C2X^2 + C3X^3。
我用np.linspace和np.polyval的数值方法为曲线生成离散的(X,Y)，然后用shapely的Point、MultiPoint和nearest_points找到最近的点，最后用np.linalg.norm求出距离。
这是一种离散曲线的数值方法。
我的问题是，如何通过分析方法找到距离并对其进行编码？
浏览 58提问于2022-10-26得票数 0
回答已采纳
2回答
我想我不像我想的那么了解which.min。
我试图在一个最小值为正的向量中找到这个现象。
TIME <- c(0.00000, 4.47104, 6.10598, 6.73993, 8.17467, 8.80862, 10.00980, 11.01080, 14.78110, 15.51520, 16.51620, 17.11680)
我想知道1到19的值z，这是上面向量时间的索引，它包含最接近但高于z的值。我尝试了以下代码：
vec <- sapply(seq(1,19,1), function(z) which.min((z-TIME > 0)))
vec
浏览 3提问于2015-05-26得票数 2
回答已采纳
1回答
我最近遇到了以下问题：
我们有一个M个连续整数的序列A，从A1 = 1: 1,2，...M开始(例如:m=8，A= 1,2,3,4,5,6,7,8 )
我们有一个集合T，由A的L_T个连续项组成的所有可能的子序列(例如L_T =3，子序列是{1,2,3}，{2,3,4}，{3,4,5}，...)。让我们称T的元素为"tiles“。
我们有一个集合S，由长度为L_S的A的所有可能的子序列组成(例如L_S = 4，子序列如{1,2,3,4}，{1,3,7,8}，...{4,5,7,8} )。
我们说，如果T中存在K个瓦片，则S的元素s可以被T的K个“瓦片”“覆盖”，使得它们的项集的并集包含
浏览 0提问于2014-04-02得票数 1
3回答
求正值和负值之和
、、、
我有一个DataFrame，如下所示。
d = {}
d['Model'] = ['M1', 'M2'] * 4 * 3 * 5
d['Support'] = ['S1', 'S1', 'S2', 'S2'] * 2 * 3 * 5
d['Group'] = ['G1', 'G2', 'G2'] * 2 * 4 * 5
d['Case'] = ['C11', 'C2
浏览 1提问于2021-01-13得票数 2
回答已采纳
1回答
求两点之间的欧氏距离C
、、、
我试图写一个函数来计算两个浮点之间的欧几里德距离。我对函数进行了编码，并调用了它，但我得到了一些错误。我怎样才能解决这个问题？
我的代码：
#include <stdio.h>
#include <math.h>
//Function Declarations
float CalculateDistance(float, float, float, float);
int main ()
{
/* Function Call */
float distance1 = CalculateDistance(float x1,float y1,float x2,floa
浏览 4提问于2022-04-20得票数 0
回答已采纳
5回答
求集合的交集
、、
我们得到了n组不同大小的整数。每个集合也可以包含重复项。我必须找到集合的交集。如果一个元素在所有集合中出现多次，则应将其添加到结果中。
例如，假设有三个集合{0,5,5,3,4} {5,2,3,5,6} {1,3,5,5,6}。给定集合的交集应为{3,5,5}
我的方法是：
1.对数组进行排序。
2.从最小的数组开始比较每个元素并更新计数。
有没有更有效的方法来找到交叉口？
浏览 3提问于2013-03-28得票数 3
回答已采纳
2回答
我想从(-\infty, a]中集成一个概率密度函数，因为cdf不是封闭形式的。但我不确定如何在C++中做到这一点。
这个任务在Mathematica中非常简单；我所需要做的就是定义函数，
f[x_, lambda_, alpha_, beta_, mu_] :=
Module[{gamma},
gamma = Sqrt[alpha^2 - beta^2];
(gamma^(2*lambda)/((2*alpha)^(lambda - 1/2)*Sqrt[Pi]*Gamma[lambda]))*
Abs[x - mu]^(lambda - 1/2)*
浏览 0提问于2012-05-28得票数 6
回答已采纳
3回答
这句话的意思是什么
// create arrays of 1M elements
const int num_elements = 1<<20;
在下面的代码中？它是特定于CUDA的，还是可以在标准C中使用？
当我printf 'ed num_elements时，我得到了num_elements==1048576
结果是2^20，那么<<运算符是不是C中幂运算的简写呢？
// This example demonstrates parallel floating point vector
// addition with a simple __global__
浏览 2提问于2011-11-05得票数 3
回答已采纳
2回答
求矩阵的特征频率
、、
我有一个大小为8*8的方阵，其中一些项是频率(ω)的函数。我想写一个在给定范围内搜索特征频率的函数，比如(0 -1 1kHz)。 我已经包含了下面的函数。这里术语'tx'，'ki1'，'ki2‘是ω的函数。为了求特征频率，矩阵的行列式应该是零。但是如果没有给出所有的值，我就找不到矩阵的行列式。 基本上，我不想给出一个频率，然后得到特征值。我想要解出矩阵，给出特征值，这些特征值就是特征频率。 你能为此推荐一些方法或函数吗？ 任何帮助或建议都将不胜感激。提前感谢！ import numpy as np
def mat(l1,l2,omega1):
浏览 10提问于2020-09-11得票数 0
1回答
我试图从Euler项目中解决一个问题: 13195的素因子是5，7，13和29。数字600851475143中最大的素因子是什么？以下是我的尝试：
def getPrimeNumbers(testNumber):
'''
For a given number it returns a list of all prime numbers from 2 to given number.
>>> getPrimeNumbers(15)
[2,3,5,7,11,13]
'''
浏览 0提问于2020-01-07得票数 -4
回答已采纳
3回答
我有一个向量v，如
v = c(0, 0, 1, 1, 1.5, 2, 2, 2, 5, 5, 5, 5, 5, 5, 10)
res = rep(0, length(unique(v)))
for(i in 1:(length(unique(v)))){
res[i] = sum(v == i) / sum(v > i)
}
> res
[1] 0.1818182 0.4285714 0.0000000 0.0000000 6.0000000
[6] 0.0000000
愿望的结果如下：
res = c(2/13
浏览 5提问于2014-02-04得票数 1
回答已采纳
1回答
关于数据: 3D数据(100x100x20)。
计算/公式和方法:我试图计算不同时差沿行、列和角度的平方差之和。
%% Grid and time paramters
%# Grid parameters
nRows=100;
nCol=100;
InitLag_Row=0;
InitLag_Col=0;
InitLag_Ang=0;
LessLags=20;
nLags_Row=nRows-LessLags;
nLags_Col=nCol-LessLags;
nLags_Ang=min(nRows,nCol)-LessLags;
%# Time parameters
T=1:1:20;
浏览 0提问于2012-07-17得票数 2
1回答
求直线与圆的交点
、、
我想弄清楚这个函数是干什么的。这是我老师给我的，我就是搞不懂公式背后的逻辑是什么，找到x和y坐标。从我的数学课上，我知道我找到拦截的公式，但它在代码中令人困惑。所以我有一些问题，他们如何定义a，b，c的公式，以及求坐标x和y的公式。
void Intersection::getIntersectionPoints(const Arc& arc, const Line& line) {
double a, b, c, mu, det;
std::pair<double, double> xPoints;
std::pair<double, doub
浏览 1提问于2016-06-10得票数 0
回答已采纳
2回答
寻找关于我在C++中解决的一个问题的反馈。这是一个leetcode问题，我用分而治之来解决这个问题。
给定整数数组nums的最大子数组，查找具有最大和的连续子数组(至少包含一个数字)并返回其和。例如:输入：-2，1，-3，4，-1，2，1，-5，4，输出: 6解释：4，-1，2，1有最大和=6。后续:如果您已经找到了\mathcal{O}(n)解决方案，尝试使用分而治之的方法编写另一个解决方案，这比较微妙。
我正在寻找在逻辑和实现方面的反馈。如果我在线搜索C++解决方案，我只看到使用数组解决问题，并且在使用迭代器时使用向量来解决问题。
任何反馈都会很好。
谢谢!
int FindMaxim
浏览 0提问于2020-01-09得票数 2
回答已采纳
1回答
在给定地图的情况下，我该如何在C中找到该地图的逆，因此f(25) =0，依此类推。有什么建议可以帮助我朝着正确的方向前进吗
i f(i)
0 25
1 11
2 24
3 0
4 19
5 2
6 7
7 15
8 27
9 14
10 17
11 22
12 13
13 26
14 20
15 16
16 8
17 10
18 18
19 4
20 12
22 23
23 3
24 5
25 21
26 6
27 1
浏览 0提问于2011-03-21得票数 1
示例如下：
我有以下矩阵：
4 0 3
5 2 6
9 4 8
现在，我想找出两个最小值，以及每一行的索引。所以结果是：
row1: 0 , position (1,2) and 3, position (1,3)
row2...
row3....
我使用了很多for循环，这是相当复杂的。那么有没有办法使用MATLAB函数来实现我的目标呢？
我试过了，但没有结果：
C=min(my_matrix,[],2)
[C(1),I] = MIN(my_matrix(1,:)) &find the position of the minimum value in row 1??
浏览 0提问于2013-01-10得票数 3
回答已采纳
2回答
如何返回此vba函数的值？如果我将debug.print添加到函数本身，它会打印正确的值。
Sub testIt()
Debug.Print getAccountNumber("keyB")
End Sub
Function getAccountNumber(myKey)
Dim c As Collection
Set c = New Collection
c.Add "4005", "keyA"
c.Add "4015", "keyB"
c.Add "4030", "keyC"
浏览 1提问于2015-03-05得票数 0
回答已采纳
我要查出车号！
看这张照片
使用此代码：
import numpy as np
import imutils
import cv2
# Read the image file
image = cv2.imread('Car_Image_1.jpg')
# Resize the image - change width to 500
image = imutils.resize(image, width=500)
# Display the original image
cv2.imshow("Original Image", image)
#
浏览 0提问于2019-03-01得票数 1
回答已采纳
15回答
这一次，您的目标是在不使用分支或循环的情况下，在二进制2's补码中找到3个整数的最大值(从-( 2^31 )到2^31-1)。
你只允许使用
不等式/等式(==、>、>=、<、<=、!=)都算作两个记号。
算术(+，-，*，/)
逻辑运算符(! not，&& and，|| or)
位运算符(~ not，& and，| or，^ xor，<<，>>，>>>算术和逻辑左、右移位)
常量。0代币
变量分配0代币
输入3个变量，如a、b和c。输出最大数目。
标准原子代码-适用高尔夫球规则。如果您有任何问题
浏览 0提问于2014-06-25得票数 17
回答已采纳
1回答
我用Matlab求出了雅可比迭代矩阵的谱半径，其中A=[4 2 1;1 3 1;1 1 4]。
在5次迭代后，我似乎无法输入正确的命令来获得错误的大小。有人能帮我吗？
以下是我到目前为止在Matlab中输入的命令列表：
A=[4 2 1;1 3 1;1 1 4]
A =
4 2 1
1 3 1
1 1 4
D=diagonal(diagonal(A));L=(A,-1);U=(A,1);
b=([3 -1 4])
x0j=zeros([0 0 0]);
x=D\(-(U+L)*x0j+b);r=b-A*x %Jac
浏览 2提问于2013-11-15得票数 2
3回答
我当时正在做K&R的练习2-1，目标是计算不同变量类型的范围，下面是我的函数，用来计算short int可以包含的最大值：
short int max_short(void) {
short int i = 1, j = 0, k = 0;
while (i > k) {
k = i;
if (((short int)2 * i) > (short int)0)
i *= 2;
else {
j = i;
while (i + j <
浏览 0提问于2017-08-25得票数 5
回答已采纳
1回答
我试着用ggplot2得到一个基于方程的图表。以下是我迄今所做的工作：
eq <- function(x,a) {
g<- (x^2+a+1)
return(g)
}
ggplot(data.frame(x = c(-10000, 10000)), aes(x)) +
stat_function(fun = eq, args = list(a=2), aes(colour = "Time"))
这给了我情节。我希望将范围保持在x轴上，即c(-10000,10000)，我想更改a，但在下面的代码中保留x以获得绘图，例如eq(x,3)。
所
浏览 2提问于2021-01-02得票数 1
回答已采纳
1回答
求凸壳指数
、、
我想计算一组点的凸包。我在网上发现的大多数算法都会返回一个点列表，但是我需要一个索引列表。
为此，我使用了一些计算点的现有代码，并尝试将其更改为返回点索引。
原始函数返回凸包点向量
#include <iostream>
#include <vector>
#include <tuple>
using namespace std;
typedef std::tuple<int, int> point;
// returns true if the three points make a counter-clockwise turn
bool
浏览 6提问于2022-02-12得票数 1
回答已采纳
1回答
求n维矩形的坐标
、、
在C++中，计算n维矩形(超矩形)坐标的最佳方法是什么？
我在一维向量中定义了矩形的维数，如：
[min1, max1, min2, max2, ...., minN, maxN]
例如，在2d中，向量是
[min1, max1, min2, max2]
我要找的角落坐标是
[min1, min2], [min1, max2], [max1, min2], [max1, max2]
对于n维的超矩形，我们如何做到这一点？
浏览 4提问于2019-07-16得票数 3
回答已采纳
1回答
我试图求解一个方程，并从它中找到x的值。我使用以下代码：
mu1 = 0
mu2 = 1
sigma1 = 0.5
sigma2 = 0.6
prior1 = 0.3
prior2 = 0.7
boundary = function(x) {
return(( 1 / sqrt(2 * pi * sigma1)) * exp(-0.5 * ((x - mu1) / sigma1)^2)*prior1) -
((1 / sqrt(2 * pi * sigma2)) * exp(-0.5 * ((x - mu2) / sigma2)^2)*prior2)
}
uniroot(bou
浏览 0提问于2014-02-06得票数 0
回答已采纳
2回答
我有以下示例数据集
Time <- c(1,2,3,4,5,6,7,8,9,10,11,12)
Value <- c(0,1,2,3,2,1,2,3,2,1,2,3)
Data <- data.frame(Time, Value)
我希望自动为Value列找到每个最大值，并创建一个只有值和相关时间的新数据框架。在本例中，最大值每四次出现一次。我想将数据分组到回收箱中，并找到相关的最大值。
不过，为了说明起见，我将我的示例保持简单，但请记住：
每个最大值在我的数据集中将是不同的每个最大值不能保证以相同的间隔发生，而是可以保证每个最大值都发生在一个时间值的范围内(即一个b
浏览 5提问于2022-05-27得票数 0
回答已采纳
2回答
我正在使用一个马尔可夫模型，用R语言编程的矩阵来解决一个健康经济问题。
我有一个矩阵，我想对它求幂200次。矩阵中的值随时间变化。例如，位置x12上的值在第一个矩阵0.1、第二个0.2和第三个0.3中。我不想写下200个矩阵，而是写下一个矩阵，并将其乘以200。对于x12:我是否可以使用一个具有200个值的向量，并且矩阵将使用该向量？
对于常量值，这是没有问题的：
############### Imatinib Base-line Strategy ##########################
trans_matrix_dasa_no2nd <- matrix (,nrow=3
浏览 0提问于2013-08-01得票数 0
1回答
用用户输入求内插
、、、、
我对算法做了一些修改，并在函数中得到了一个值(尽管不正确)，但是返回值拒绝将它发回主函数。此外，我也不能让他们
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <float.h>
/*Prototype for functions used*/
float getdistance(float[], float[], float, int);
float printvalue(float);
int main()
{
int maxtime = 18
浏览 0提问于2014-04-01得票数 0
回答已采纳
2回答
反求等距投影算法
、、、
我有这样的代码：
const a = 2; // always > 0 and known in advance
const b = 3; // always > 0 and known in advance
const c = 4; // always > 0 and known in advance
for (let x = 0; x <= a; x++) {
for (let y = 0; y <= b; y++) {
for (let z = 0; z <= c; z++) {
for (let p = 0; p <=
浏览 4提问于2022-07-13得票数 1
回答已采纳
1回答
我正在写一个程序：
例如，输入是5(可以不只是5)数字，并且我读取数组中的数据：1, 2, 3, 4, 5。我可以从这个数组中选择一些元素(不是第一个或最后一个)，例如3，然后在数组中删除这个数字，并将从左到左相乘后的元素与从右到右相乘的元素(在本例中意味着2*4 )添加到sum (最初是0)。得到的数组是1, 2, 4, 5，然后我一次又一次地做，直到元素的数量等于2(确切地说是1 and 5，因为我们不能删除这些数字)。
例如：(其中A、B、C、D是数字1和2、2和3等的对。)
A B C D
1 2 3 4 5
订单删除元素有6种可能的组合(并将左右乘法添加到sum中)：
A (B (
浏览 3提问于2012-11-28得票数 3
回答已采纳
1回答
我正在试图找到我在一个锯齿状的2d数组上的平均行，我已经打印出了这个数组，但是平均部分是打印错误的平均值。我的密码怎么了？
//print array
for(int k = 0; k < arr.length; k++) {
for(int j = 0; j < arr[k].length; j++) {
System.out.print(arr[k][j] + " ");
}
System.out.println();
}
//average
for (int r =
浏览 1提问于2022-02-01得票数 1
回答已采纳
1回答
我对安全问题很陌生，我对转位密码学的密钥也有一些疑问。
如果给我一个转位密码的编码密钥，我如何获得它的解码密钥？我试着在googling上搜索，但没有找到计算的任何步骤。例如，在几个与密码相关的网站上，有一些很棒的工具，当我进入4 3 1 6 2 5时，我会得到3 5 2 1 6 4…。但为什么和怎么做？
浏览 0提问于2016-03-04得票数 1
4回答
我必须得到一个素数序列。但是我的代码不能工作。怎么可能修复它呢？
var num1 = parseInt(prompt('Enter a number'));
var num2 = parseInt(prompt('Enter a number'));
var num3 = 0;
function primeSeq(num1, num2) {
var b = 1;
var c = '';
if (num1 > num2) {
num3 = num2;
num2 = num1;
num1 = num3;
}
for (va
浏览 0提问于2016-11-07得票数 0
1回答
我有一个数据集(x,y)，其中我用三项拟合高斯(gauss3)：
[gaussianFit,gof]=fit(x,y,'gauss3');
现在，我想找出这个函数的最大值。我知道对于gauss1，最大值只是b1参数，但是当它是多个高斯的和时，我不知道如何找到最大值。这些术语似乎是有权重的(否则，我的想法只是简单地取三个b-terms的平均值)。我的第二个想法是用所有参数来求解整个微分函数：
syms x
diff(a1.*exp(-((x-b1)./c1).^2) + a2.*exp(-((x-b2)./c2).^2) + a3.*exp(-((x-b3)./c3).^2)==
浏览 1提问于2015-04-29得票数 0
1回答
使用R，我想自动找到收益递减曲线的切线，以便进行最佳的搜寻图形练习。曲线如下所示：
ins <- 40
t <- 30
avg <- 30
curve(ins/2^(x/10+1-1),
0,2*t,xlim=c(-2*t,2*t),ylim=c(ins,0),
xlab="time",ylab="food",type="l",lty=1,col=4,lwd=3,axes=FALSE)
axis(1,pos=ins); axis(2,pos=0)
切线从点(-40,40)开始，理想情况下，我会在曲线上找到切线接触
浏览 0提问于2018-05-21得票数 0
1回答
我正在尝试下一段代码，看看predict是否可以帮助我找到2阶多项式的因变量的值，在这种情况下，它是明显的y=x^2：
x <- c(1, 2, 3, 4, 5 , 6)
y <- c(1, 4, 9, 16, 25, 36)
mypol <- lm(y ~ poly(x, 2, raw=TRUE))
> mypol
Call:
lm(formula = y ~ poly(x, 2, raw = TRUE))
Coefficients:
(Intercept) poly(x, 2, raw = TRUE)1 poly(x, 2, raw
浏览 0提问于2012-09-26得票数 6
回答已采纳
我将一段用fortran 77编写的代码转换为Matlab代码。此函数使用QL算法计算矩阵的特征值和特征向量。由于某些原因，我不能在matlab中使用eig函数的结果。用该方法得到的特征值与用eig函数得到的特征值不完全相同，有的相同，有的不同。我不知道问题出在哪里。感谢您的帮助和建议。如果运行和查看结果需要，我可以给出输入数组。
下面是fortran代码：
SUBROUTINE tqli(d,e,n,np,z)
INTEGER n,np
REAL d(np),e(np),z(np,np)
CU USES pythag
INTEGER i,
浏览 0提问于2012-07-29得票数 0
回答已采纳
1回答
求p值的简单方法
、、、
我不是数学，请温柔点..。
我正在寻找一个简单的方法，以求p-值对应于给定的x-平方(一自由度)。我想把它放到一个计算机程序(FileMaker)中；准确性并不重要。
从这次讨论中：
我有三个公式：
第一个公式：
p(x)=a/x^3 + b/x^2 + c/x +d
其中x=x-平方和a，b，c和d是常数。
范围从0,1到3，然后:A 0.00350853 /b -0.0623318 /c 0.357489 /d -0.0049905
范围为3至12，则为:a -3.1258 /b +2.81657 /c -0.387984 /d 0.0152614
第二个公式：
您可以使用与一个df的
浏览 0提问于2018-08-21得票数 0
1回答
让我们考虑一下这个Dataframe：
$> df
a b
0 6 50
1 2 20
2 9 60
3 4 40
4 5 20
我想计算列D的基础是：
index中B列的最大值：integer 0A片段
因此，我在我的数据中创建了一个列C(全部为零)，以便使用DataFrame.max(axis=1)。但是，除了在DataFrame上使用apply或循环之外，我不知道如何分割输入值。预期结果将是：
$> df
a b c d
0 6 50 0 60
1 2 2
浏览 3提问于2022-02-08得票数 0
回答已采纳
1回答
例如：a，b，c和n是三个常数，在特定的方程中需要用数据拟合的方法来计算。
如何计算参数的统计量(均值、标准差、方差、偏差值和学生t检验值)作为自定义方程(例如二次平台方程)？
示例：
x=[0,40,80,100,120,150,170,200],
y=[1865,2855,3608,4057,4343,4389,4415,4478]
y=a*(x+n)^2+b*(x+n)+c, x < xc(Ymax) ....(1) y=yp, x >= xc(Ymax) ....(2)
我已经通过给定的代码拟合了这个方程：
yf = @(b,x) b(1).*(x+n).^2+b(
浏览 9提问于2018-12-31得票数 0
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
即时通信 IM
对象存储活动推荐
腾讯云开发者
